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Abstract 



In this work we give a theoretical description of the elastic vector meson production in diffrac- 
tive DIS developed within the /cj-factorization formalism. Since the /c^-factorization scheme 
does not require large values of Q"^ +my, we conduct an analysis that is applicable to all values 
of from photo- up to highly virtual production of vector mesons. The basic quantity in 
this approach — the unintegrated gluon structure function — was for the first time extracted 
from the experimental data on F2p, thoroughly investigated, and consistently used in the vec- 
tor meson production calculation. Moreover, by limiting ourselves to the lowest Fock state 
of the vector meson, we were able to construct in a closed form the theory of spin-angular 
coupling in the vector meson. This allowed us for the first time to address the production of a 
vector meson in a given spin-angular state. We performed an extensive analytical and numer- 
ical investigation of the properties of IS, 2S, and D-wave vector meson production reactions. 
Treating the physical ground state vector mesons as purely IS states, we observed a good 
overall agreement with all available experimental data on vector meson production. For the 
excited states, our analysis predicts a picture which is remarkably different from IS'-state, so 
that such reactions can be regarded as potential sources of new information on the structure 
of excited states in vector mesons. 
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Chapter 1 
Introduction 



In the past 30 years the particle physics theory has proved numerous times to provide a 
good, consistent, unified description of the great variety of nuclear, low and high energy 
particle physics experiments. Being based on the ideas of QFT applicability, gauge approach 
to fundamental interactions, symmetry and naturalness considerations, the Standard Model 
managed to explain virtually all phenomena in electromagnetic, weak and strong interactions, 
to predict new particles and effects. Though the questions of fundamental origin lie beyond 
the scope of the Standard Model, its precision in description, for instance, QED phenomena 
reaches the magnitude of 10~^°. 

However, the current situation is not that optimistic in the domain of strong interactions. 
The gauge-based formulation — the Quantum Chromodynamics (QCD) — seems to offer 
reasonably good description only of the energetic enough processes (more accurately: only 
when every vertex involves at least one highly virtual particle) thanks to the famous asymptotic 
freedom. The major difficulty lies in the behavior of the QCD coupling constant: as{Q'^) 
exhibits infrared growth and becomes comparable to unity at ~ 1 GeV^. The net result 
is that the perturbation theory — the only prolific universal treatment of various high-energy 
processes — fails to give even qualitative description of low-energy, essentially non-perturbative 
phenomena. Additional difficulties arise from the non-abelian nature of QCD, chiral symmetry 
breaking, non-trivial QCD vacuum, instantons etc. 

On the other hand, many separate concepts have been developed, which do not cling to 
the perturbative QCD (pQCD) and provide reasonably good description of phenomena in 
their applicability regions. The fundamental problem of the theory of strong interactions is 
that these heterogeneous approaches do not match^. They do not comprise a unified picture 
of strong interactions. Given such a lack of universal, rigorously derived results, one must 
admit that the subject of our investigation belongs to the realm of phenomenology rather 
than rigorous theory. 

^ Just a few examples of poor accordance among various approaches: the quark generated ladder diagrams 
do not appear to correspond uniquely to any of c;xperimentally observed Regge trajectories. Another example 
is the vague status of = const BFKL results in true QCD. 
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1.1 Diffractive processes and Pomeron 



In the light of these problems, the careful examination of regions where two or more approaches 
overlap (or visa versa, where neither of the concepts exhausts the interaction) are of great 
interest. Diffractive Deep Inelastic Scattering (DIS) is exactly one of these fields. 




a) P P B) 



Figure 1.1: Examples of deep inelastic scattering process: (a) hard DIS and (b) diffractive 
DIS. In the latter case Mj^^j:^ <S s and the process proceeds via pomeron ^-channel exchange. 

A typical hard DIS process (Fig. 1.1a) occurs when a virtual photon^ strikes a proton to 
produce a hard system X with large invariant mass'^ s and large enough multiplicity, final 
state hadrons being distributed over whole rapidity range approximately smoothly. 

However, as it was noted long ago, sometimes the proton survives, being only slightly 
deflected, and a virtual photon turns into a so-called diffractive system Xdiffr with invariant 
mass Afjjjj^ <^ s. In this process the proton and the diffractive system are naturally separated 
by a large rapidity gap and a condition which appears necessary for the rapidity gap formation 
is -C s, or in terms of Bjorken x 

x=^«l. (1.1) 
s 

This is one of the most common cases of diffractive DIS (DDIS) processes. In fact, the 
class of diffractive processes is not confined within DIS; it is much broader. There are many 
other reactions which possess the generic features — the rapidity gap and smallness of Mj^^y^. 
— and therefore can be classified as diffractive processes (for a recent review see [1]). 

How can a typical diffractive process occur? Certainly, it must be kind of a peripheral in- 
teraction: if the photon struck directly one of the valence quarks, the proton would 'explode', 
providing no way for the large rapidity gap formation. What remains is the possibility of 
the t-channel exchange by not-too-energetic 'particle' (Fig. 1.1b), which would be a natural 
mechanism of the experimentally observed weak proton deflection and small Mj^jj^. Further 
experimental features suggest that this 'particle' should be chargeless and colorless, its inter- 
action AA'itli otlier ])aitieles sliould be of strong (not EM or weak) nature, its 'propagation' 

^We will always imply that the virtual photon is emitted by an electron, which means the photon is always 
space-like: if q is photon momentum, then = —q^ > 0. 

^In hard DIS phenomenology this quantity is usually labeled as W"^. However, for simplicity we will use 
notation s. 
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should be independent of the specific process (7^, 77, pp, pp, etc), and it should be of spin 1 
(due to approximately s-constant pp cross section). In the early 60s this 'particle' was dubbed 
Pomeron (symbol IP). 

Further properties come from combining the Regge picture and BFKL results with ex- 
perimental observations (for a detailed review of Regge theory see [2]). They include, first 
of all, the asymptotic equality of total pp and pp cross sections (the Pomeranchuk theorem). 
Formulated long ago, it was experimentally verified only recently. Then, the Regge theory 
predicts the power-like s-dependence of the total pp cross section a (x ^^(aip-i)^ which has also 
been experimentally observed, with intercept Sjp = ajp — 1 fa 0.08. On the other hand, the 
BFKL equation [3, 4] succeeded in reproducing such power- like dependence in QCD, but in 
a simplified case as =const. In this approach the hard pomeron is treated as two reggeized 
gluons — an ansatz used currently in diffraction phenomenology with great success. However, 
the predictive power of the BFKL approach for the numerical value of the pomeron intercept 
is still limited and not all issues with sensitivity of the result to the infrared region have been 
understood. For further reading on pomerons, a topic very intriguing by itself, we refer to [5] . 

1.2 Vector meson production in diffractive DIS 

There arc several possible final states X in a typical diffractive DIS (DDIS) process 7*j9 Xp: 
system X can be a real photon, a qq continuum pair forming two jets or qq bound state, for 
example, a vector meson. Let us now focus specifically on exclusive vector meson production 
in diffractive DIS. This reaction has been studied extensively at fixed target DIS experiments 
at CERN and FNAL and more recently by the HI and ZEUS collaborations at HERA. 

Despite the great deal of theoretical work on vector meson production in diffractive DIS 
[6, 7, 8, 9], there is a number of issues that have not yet been seriously analyzed and still need 
a closer investigation. 

One of them concerns the vector meson production in a seeming soft region, namely, 
at small values of and m\r (or, to put it short, at small values of + ^y)- Indeed, 
the majority of early calculations treated the vector meson productions within the DGLAP- 
inspired approach, the production amplitudes being expressed in terms of the intergrated 
gluon density G^x^Q"^). Certainly, this line of calculation is not applicable at small enough 
values of + my, say at -\- my <j 1 GeV^. However it is necessary to understand that 
the DGLAP-based approach not only can be avoided but also should be avoided when one 
studies diffractive scattering. Indeed, at high energies and small to moderate values of the 
dynamics of the amplitudes is governed by large logarithms of log(l/a;) rather than log((5^), 
and the correct and the most natural method to treat processes in this kinematical region is 
/c^-factorization approach. This approach does not place any restriction on the value of as 
long as one works at large enough energies. 

Although the strategy of the evaluation of the vector meson production amplitudes within 
the /cj-factorization approach is essentially clear, performing reliable numerical prediction is not 
a straghtforward task. The impediment consists in presence of purely soft, non-perturbative 
quantities in the calculation, namely, the gluon content of the proton and the wave function 
of the vector meson. 

The presence of the former quantity is a specific feature of the particular final state we 
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investigate, however, the unintegrated gluon density 

is the basic quantity in all /ct-factorization calculations. Unfortunately, no reliable Ansatz or 
parametrization has been developed, and this gap needs to be filled. 

Another issue that has never been brought under scrutiny is the spin-angular coupling 
inside the vector meson. In an off-forward scattering 7^^ — > Vxy the s-channel helicity flip 
amplitudes can be non-vanishing. Because of the well known quark helicity conservation in 
high energy QCD scattering, such a helicity flip is possible only due to the internal motion 
and spin-angular momentum coupling of quarks in a vector meson. This issue was accurately 
analyzed only in very recent papers [6, 7], where it was shown that helicity non- conserving 
amplitudes are not negligible, as they had been thought before. Thus, as such, the helicity 
flip amplitudes would offer a great deal of unique information of internal constituent motion 
and spin-angular momentum structure of vector mesons, unaccessible in other experiments. 
In addition, the vector meson decays are self-analyzing and the full set of helicity amplitudes 
can be measured experimentally. For unpolarized incident leptons, the angular distribution 
of decay products is parameterized in terms of 15 spin-density matrix elements, which can 
be calculated via five — two helicity conserving plus three helicity violating — basic helicity 
amplitudes [10]. 

Certainly, the helicity structure of the vector meson production amplitudes must be ana- 
lyzed only along with a careful treatment of the spinorial structure of the qq V transition. 
It is thus rather surprising that the above issue of sensitivity of the production amplitudes 
to the spin-angular momentum coupling has not been addressed before. Namely, in a typical 
vector meson production calculation, a vector meson has been implicitly taken as IS state 
and at the same time an unjustified ansatz was used for qq ^ V transition spinorial structure, 
namely, of u'^fj,u ■ type. Being a mere analogy of qq'-j vertex, this ansatz in fact corresponds 
neither to pure S nor to pure D wave state but rather to their certain mixture. Only in [11] 
the cases of 15* and 2S vector mesons were compared and the necessity of similar calculation 
for D wave states was stressed. Such calculations however have been missing in literature 
until now. 

In addition to purely theoretical needs, there arc more issues that call upon a thorough 
analysis of the D-wave effects. For instance, different spin properties of the S- and D-wave 
production may resolve the long standing problem of the D-wave vs. 25'-wave assignment for 
the p'(1480) and p'(1700) mesons (as well as the lo' and (f)' mesons). Furthermore, the deuteron 
which is a spin-1 ground state in the pn system is known to have a substantial D wave admix- 
ture, which mostly derives from the tensor forces induced by pion exchange between nucleons. 
Recently, there has been much discussion [12] of the nonperturbative long-range pion exchange 
between light quarks and antiquarks in a vector meson, which is a natural source of the S-D 
mixing in the ground state p and cu mesons. 



In the present work we addressed both issues. We performed an accurate determination of 
the unintegrated gluon density from the experimental data on the structure function and 
gave its convenient and ready-to-use parametrizations. Besides, we constructed the consistent 
description of the vector mesons with spin-angular coupling taken into account, which enabled 
us to calculate diffractive production amplitudes for pure S-wave and D-wave states as well 
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as for an arbitrary S/ D wave mixture. This resulted in a complete theory of the vector meson 
production in diffractive DIS within the fcj-fectorization approach. 

1.3 The strategy of the thesis 

The guideline of the thesis is the following. The main text is comprised of three Parts. Part 
I is an introduction to the /crfactorization approach. Here we calculate some basic scattering 
processes, such as the virtual Compton scattering, and introduce the concept of the differential 
densities of partons. The discussion on the similarities and distinctions between the DGLAP- 
motivated description and /ct-factorization description of diffractive processes can also be found 
here. 

In Part II we turn to the vector meson production amplitudes. These are preceeded by 
the theory of vector meson structure within the truncated Fock space, that is, when the 
vector meson is assumed to be a bound state of qq pair only. Upon obtaining the closed 
analytical expessions for vector meson production amplitudes, we perform the twist expansion 
and illustrated some of the most salient properties of the 5'-wave and D-wscve vector meson 
amplitudes. 

Part III contains the numerical analysis of the expressions obtained and the concrete pre- 
diction to various experimentally observed quantities. At first we perform an extraction of 
the differential gluon density of the proton and thoroughly investigate its properties. Having 
brought the differential glue under control, we turn to the vector meson production ampli- 
tudes and give a large number of predictions for IS", 25", and D-wme states. Whenever the 
experimental results are available, we confront our predictions with the data. This Part is 
concluded with a detailed consideration of the effect of the Coulomb tail of the vector meson 
wave function and of the S/ D-waye mixing. 

Finally, we summarize our main findings in Conclusions. Some lengthy calculations can be 
found in Appendices. 

All the results presented in this thesis were derived by the author. The text is based on 
publications [13] and [14], and on the works in progress [15]. Preliminary results have been 
presented as talks at [16, 17]. 
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Part I 

Basics of A;t-factorization scheme 
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Chapter 2 

The virtual Compton scattering 



We start our acquaintance with the A;t-factorization scheme with calculation of the imaginary 
part of the forward virtual Compton scattering amplitude. By means of optical theorem it 
is related with the total photoabsorption cross section and with the structre functions of the 
proton. 

The purpose to get started with this quantity is twofold. First, during this calculation we 
will follow all steps and discTiss all major feature of the /cj-factorization scheme of calculations. 
Being rather simple, the Compton scattering amplitude will keep us from being distracted by 
inessential technical complications, which would arise in other diffractive reactions. 

The second purpose is to derive the well-known expression for the structure function in 
terms of unintegrated gluon density of the proton: the basic quantity in any /cffactorization 
calculation. These expressions will be used later, when we discuss the determination of the 
unintegrated glue from the experimental data. 

The third aim is to use the simplicity of this amplitude to gain as much insight into the 
dynamics of photon-proton peripheral interactions. This information will be used later in 
deriving the vector meson production amplitudes thanks to a remarkable similarity between 
the virtual Compton scattering and the vector meson electropro duct ion processes. Indeed, in 
the proton rest frame, both can be viewed as follows: a photon dissociates into a qq pair, which 
interacts with gluon content of the proton and then is projected onto the final state. The hard 
dynamics in both cases arc the same, only difference lurking in the final state projection. 

2.1 Modeling virtual photoabsorption in QCD 

The quantity which is measured in deep inelastic leptoproduction is the total cross section of 

photoabsorption 7*77 X summed over all hadronic final states X, where yU, z/ = ±1,0 are 
helicities of (T) transverse and (L) longitudinal virtual photons. One usually starts with the 
imaginary part of the amplitude A^^ of forward Compton scattering 7*]? — > 7*p', which by 
optical theorem gives the total cross cross section of photoabsorption of virtual photons 
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where W is the total energy in the 7*p c.m.s., m.p is the proton mass. is the virtuahty of 
the photon and = Q"^ / (Q^ + — nip) is the Bjorken variable. Hereafter we will suppress 
the subscript bj and use x = x^j. 




P 



Figure 2.1: The pQCD modeling of DIS in terms of multiproduction of parton final states. 

In pcrturbative QCD (pQCD) one models virtual photoabsorption in terms of the multiple 
production of gluons, quarks and antiquarks (fig. 2.1). The experimental integration over the 
full phase space of hadronic states X is substituted in the pQCD calculation by integration 
over the whole phase space of QCD partons 

/ iM^l^dr^ ^ Y: \Mn\' n /' —d'i^i , (2-3) 

where the integration over the transverse momenta of partons goes over the whole allowed 
region 

0</?2< 1^.2^ Q!(l^_ (2.4) 
~ ' ~ 4 Ax ^ ^ 

The core of the so-called DGLAP approximation [18] is an observation that at finite x the 
dominant contribution to the multiparton production cross sections comes from a tiny part of 
the phase space 

1 > > X2--- > Xn-l > Xn > X , 

0<«?<«i...<«2_i^P<g2 ^ (2.5) 

in which the upper limit of integration over transverse momenta of partons is much smaller 
than the kinematical limit (2.4). At very small x this limitation of the transverse phase space 
becomes much too restrictive and the DGLAP approximation is doomed to failure. 

Hereafter we focus on how lifting the restrictions on the transverse phase space changes 
our understanding of the gluon structure function of the nucleon at very small x, that is, very 
large - . In this kinematical region the gluon density g{x, Q"^) is much higher than the density 
of charged partons q{x,Q'^),q{x,Q'^). As Fadin, Kuraev and Lipatov [19] have shown, to the 
leading log - (LL-) approximation the dominant contribution to photoabsorption comes in 



16 



this regime from multigluon final states of fig. 2.1; alternatively, to the LL- splitting of gluons 
into gluons dominates the splitting of gluons into qq pairs. As a matter of fact, for the purposes 
of the present analysis we do not need the full BFKL dynamics, in the /Cffactorization only the 
qq loop is treated explicitly to the LL- approximation. In this regime the Compton scattering 
can be viewed as an interaction of the nuclcon with the lightcone qq Fock states of the photon 
via the exchange by gluons, fig. 2.2, and the Compton scattering amplitude takes the form 

^.M = *:,AA®^55-®*M,AA (2.6) 

Here \E'^ is the and q, q hclicity A, A dependent lightcone wave function of the photon 
and the QCD pomeron exchange gg-proton scattering kernel Agg does not depend on, and 
conserves exactly, the q,q helicities, summation over which is understood in (2.6). 




WWC3 




3ww wwC 




Figure 2.2: The kt-factorization representation for DIS at small x. 



The resummation of diagrams of fig. 2.1 defines the unintcgratcd gluon structure function 
of the target, which is represented in diagrams of fig. 2.2 as the dashed blob. 



2.2 Details of calculation 

Suppose that there were no interaction between the gluons exchanged in i-channel, so that 
the full calculation of the Compton scattering amplitude amounts only to picking up the Born 
diagrams. Consider one of such diagrams, i.e. Diagr.a in Fig. 2.2, but without the dashed blob. 
A virtual photon turns into qq pair which interacts with a proton via two-gluon exchange. The 
general expression for this amplitude reads: 

f d^k f d^K ... ,/ R/, . f) — k + m , . „/ 4/, 

lA = / — — / — — u\-ig^^ t^ )t ^ ^^(-^g^^^ )u 

J [2n)* J (27r)* [{p — ny — + ie\ 

■y—i)^ 2 , . ■ \—'>')^ 2 , • ■ ^/ 

— + le K'^ — jji'^ -\- le •' 
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Sp |(— ie)e i{k — q + m) {—ie)e* i{k + m) {—igYt^) i{k + k + m) {—ig'-)H^) i{k + m)| 



(2.7) 



[k'^ — w? + ie] [{k + — + ie] [{k — qY — rn? + ie] 

Let's first calculate the numerator. 
Color factor 

If we consider photon scattering off a single quark, we have 

^Mt-t-) . = = \^(Nl -^) = \C.= 1 (2.8) 

However, we should take into account that quarks are sitting inside a colorless proton, whose 
color structure is 

Scalar = -^f^^^^q' (2.9) 

In this case there are two ways a pair of gluons can couple 3 quark lines. In the first way 
both gluons couple to the same quark. Since the quark momentum does not change after 
these two interactions, the nucleon stays in the same state: {N\N) = 1. In the second case, 
gluon legs are attached to different quark lines, so that extra momentum k circulates between 
quarks, which gives rise to the factor (A^| exp(i«;ri — iKr2)\N) , i.e. to the two-body formfactor. 
Therefore, for the lower line instead of 

^Mt^t^} = (2.10) 

one has 

^e^'^ {35a.'Swti.tf.^, + 66aa'CaN\ exp(mri - i«r2)|7V)) e"'^'^' 
= Sp{t^t^} - Sp{t^t^}{N\ exp(iKri - irera) | A^) 

= ^6ABil-{N\exp{^f^ri-inr2)\N)). (2.11) 

Note also that a similar calculation for Nc number of colors would yeild the same result. Thus, 
the overall color factor is 

2V{k) = 2(1 - {N\ exp(mri - iKr2)\N)). (2.12) 

As known, the highest power s contribution comes from so-called nonsense components of 
gluon propagator (density matrix) decomposition: 

2p'q', 2p'q', , 2p'q', 
9m' -—^^—^ + 9^,^l' —■ l^.idj 

Therefore, the lower trace is calculated trivially 

\sp{p'c^p'q'}^s\ (2.14) 



18 



So, combining all factors, one has for numeractor of Eq.(2.7) 

4 . . . . . 

{AnasY 47rQ;em ■ 2— ■ Sp |e {k — q + m) e* {k + m) q' {k + k + m) q' (k + m)j 

= {47ra,f47ra,me}-8-2s^-l(''\j*^j*). (2.15) 

Note that we factored out because it will appear later in all trace calculations. So, the 
resulting expression for amplitude (2.7) looks like 

A = (Ana,)Una,^ 16e} ■ [ -^J ■ ^ ^—1^-16) 

J [ZttP J (zTTj^ [{p — K)'^ — + ie\ [K^ — fl^ + l€\ 

hr ^ 7*) 



[k'^ — w? + ie] [{k + kY — + ie] [{k — qY — w? + ie] 

One can now immediately write similar expressions for three other diagrams (Fig. 2.2 
b,c,d). Indeed, they will differ from Eq.(2.16) only by the last line. Aside from different 
expressions for traces, the quark line propagator structures will read: 

(6) [k"^ - + ie] [{k - qf - + ie] [{k + nf - + ie] [{k + k - qf - + ie] 

(c) [k"^ — + ie] [{k — qY — rn^ + ie] [{k — kY — + ie] [{k — k — qY — rn^ + i^] 

(d) [k^ - + ie][{k - qY - + ie][{k - kY - + ie][{k - qY - + ie] (2.17) 



2.3 Denominator and trace evaluation 

Now we turn to calculation of denominators. As usual, we implement Sudakov's decomposition 
and 

k = yp + zq' + k; 
K = ap' + Pq' + K 

q = q — xp , (2-18) 

and make use of relation 

d!^k ^ -sdy dz (fk . (2.19) 

The complete analysis of denominator hierarchy and their integrals is performed in Appendix 
A. We show there that, for example, for diagram A, the imaginary part of the desired integral 
is equal to 

^ f f dy dz da dB 1 An"^ r'^ , 1 - z 1 1 

Iml — ^l = — / dz ^ r^. 2.20 

[J propagators j Jo z [k^ + rn'' + z{l- z)Q^]'' [^^ + //^J^ 

The answers for the other three diagrams differ only by replacements k ^ k + k in quark prop- 
agators whenever appropriate. The whole expression for the imaginary part of the amplitude 
is then 



32 2 f J j2r f d'^^iV{K)a 



ImA = s ej ■ Uem ■ / dzd'^k j 



(27r)2^/ J 'V («2 + ;,2)2 
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X 



+ 



J (a) ^ 
Jib) + /(c) 



(2.21) 



[{k + kY + m? + z{l - z)Q^] [A;2 + + z{l - z)Q^] 
The short-hand notation should be in fact understood as 

Q;s(lower) ■ Q;s(upper) = as{K^) ■ OLs{q^) ■ (2.22) 

with 

= max[P + + z{l - z)Q^, tf] . (2.23) 

Now we calculate the integrands / which enter Eq.(2.21) Sp{...} = 2s^ ■ I. We will do 
this via light cone helicity amplitude technique. In the subsequent discussion we will use the 
following convention: 

q — q+rij^ 

i.e. the light cone direction + is taken along photon propagation. 

A crutial point that justifies the usage of the helicity amplitude technique for all quarks 
lines inside the loop is that in the trace calculation all fermions can he treated as on-mass shell 
thanks to the presence of n_ vertices. This property comes from the following arguments. 

Note that every intermediate quark line in any diagram couples at least to one of the 
t-channel gluons. Algebraically, it means that every k + m stands near the factor 7t,_. Let us 
apply the Sudakov decomposition to the 7 matrix: 

Y = 7+^+ + i-n- + T ; 

7+ = n_ = ;^(7o + 73) , 7- = ^+ = ;^(7o - 73) ■ 
Now decompose the propagator numerator of the constituent, to which this n_ leg couples: 

A — * 

k -\-m — A;+7_ + A;_7+ — k^-\-m (2.24) 
and rewrite it using notation of (4.9), (4.11) as 

k + m^ k+-f^ + k*_-f+ -k^ + m + (k^- k*_)-f+ ^k* + m-{ 7+ . (2.25) 

2k+ 

In other words, we expressed the virtual quark propagator as the sum of on shell quark 
propagator and an additional "instantaneous interaction" term. However, since n_ is inserted 
between two {k + m) factors, this item does not work due to identity 7+7+ = 0. The net 
result is that wherever n_ appears, both constituents can be treated on mass shell in the trace 
calculation, which completes the proof. 

Having established that the fermion lines in the trace calculation can indeed be taken as 
if the quarks were real, we can now decompose the numerator of each ofthe quark lines as 

k + m ^ k* + m — '^u\Ux , (2.26) 
where spinors ux are for an on-mass shell fermion. 
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In the case of antiquark line the above derived property is vahd as welL The only thing to 
remember here is that antiquark propagates upstream the fermion arrow, so that 



X 

i.e. each antiquark propagator gives rise to factor —1. 

The derivation is given in Appendix B in full detail and yields 
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X < m 



Yl 4 • "em ■ J dz (fk J 



2f f d'^KV{K)ag 



1 



{k + + + z{l - z)Q'^ k"^ + m? + z{l - z)Q'^ 

k + K 



x + [z^ + {i- zy] 



(2.27) 



(2.28) 

2, 



{k + ay + m2 + z{l - z)Q^ k^ + m^ + z{l - z)Q^' 
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(2vr)2 , 
x4z'(1-^)2q2 



- f d^KV{K)a 



+ /i2)2 



1 



1 



(A; + k)2 + m2 + z{l - z)Q^ k^ + + z{l - z)Q^ 



.(2.29) 



2.4 Gluon density 

It is obvious that Eqs.(2.28) and (2.29) are not directly related to the real experimental 
situation, for up to now we assumed that the two exchanged gluons do not interact. Such 
interaction will definitely change the properties of the entire t-channel exchange, and in fact, 
as predicted by the BFKL equation, the resultant Pomeron has rather little in common with 
the initial two perturbative gluons. 

Since the BFKL evolution necessarily involves soft gluons (see below more on soft-to-hard 
diffusion), it does not allow for accurate perturbative calculations. It must be understood 
however that although we do not know what happens "inside the Pomeron" on the way from 
proton to the quark-antiquark pair, we nevertheless know — and the knowlegde is based on 
the leading order BFKL analysis — that eventually the qq pair will interact with nothing else 
but two gluons. We underline that this conclusion does not require the gluons to be hard, but 
rather it relies on the fact that higher Fock states of the t-channel can be, to the leading log ^ 
approximation, absorbed in the two gluon state [20]. 

Thus, the only thing we need to know is the momentum distribution of the uppermost 
gluons, or, to put it exactly, the probability distribution to find a gluon with given lightcone 
momentum fraction Xg and the transverse momentum k. 

dng = J^{K,Xg)^^. (2.30) 

This distribution is called the unintegrated (or differential) gluon structure function, DGSF, 
or simply the unintegrated gluon density. 
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Since the differential gluon density is uncalculable with pQCD, a reasonable way to proceed 
in our computation of the compton scattering amplitude further consists in finding out the 
correspondence V{k) <-> J-'{k,, Xg). Namely, we will calculate the unintegrated gluon density at 
the Born level J- Bom in terms of V{k), and then postulate that the BFKL dynamics amounts 
to replacement Tsom ^- This procedure will give us a unique prescription how to correctly 
incorporate the unintegrated gluon density into the /ct-factorization calculations. 

In order to provide a gentle introduction into the concept of the unintegarated parton 
densities, we start with famous Fermi- Weizsacker- Williams approximation in QED. We will 
find the expression for the unintegrated photon densities in the case of a single charged particle 
and charge neutral positronium and then translate the results to the case of color forces. 



2.4.1 Differential density of gauge bosons: the QED primer 




Figure 2.3: The Fermi-Weizsdcker- Williams diagram for calculation of the flux of equivalent 
photons 

For the pedagogical introduction we recall the celebrated Fermi- Weizsacker- Williams ap- 
proximation in QED, which is the well known precursor of the parton model (for the review 
see [21]). Here high energy reactions in the Coulomb field of a charged particle are treated as 
collisions with equivalent transversely polarized photons — partons of the charged particle. 
Fig. 2. 3. The familiar flux of comoving equivalent transverse soft photons carrying a lightcone 
fraction <^ 1 of the momentum of a relativistic particle, let it be the electron, reads 

7 ^ ^em ^ dx^ ^em dXry . n 

an — . — , l^-3-l-j 

TT (K^ -|- K-^^Y X-y TT K'^ X-y 

Here k, is photon transverse momentum and the photon longitudinal momentum 

in the electron Breit frame. The origin of in the numerator is in the current conservation, 
i.e. gauge invariance. Then the unintegrated photon structure function of the electron is by 
definition 

77 f ->2\ ^^e <^em ( '^'^ \ (r, onN 

^ ^ ai^ ^ ""-^ dx.diogK' ^ ~ir [w^j ■ ^ ^ ^ 

2.4.2 Differential density of photons in a positronium 

If the relativistic particle is a positronium. Fig. 2.4, destructive interference of electromagnetic 
fields of the electron and positron must be taken into account. Specifically, for soft photons 
with the wavelength A = - ^ op, where op is the positronium Bohr radius, the electromagnetic 
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fields of an electron and positron cancel each other and the flux of photons vanishes, whereas 
for X <^ ap the flux of photons will be twice that for a single electron. The above properties 
are quantified by the formula 

:F^{x,,n^)=N,^ y^Jvin), (2.33) 

where factor A^c = 2 is a number of charged particles in the positronium and corresponds to 
the Feynman diagrams of Fig. 2.4a, 2.4b. The vertex function V{k) is expressed in terms of 
the two-body formfactor of the positronium, 

V{k) = 1 - F2{k, -k) = 1 - (P| exp(iK(r_ - r+)) |P) , (2.34) 

where r_ — r_|_ is the spatial separation of e"*" and e~ in the positronium. The two-body 
formfactor F2{k,—k) describes the destructive interference of electromagnetic fields of the 
electron and positron and corresponds to the Feynman diagrams of Fig. 2.4c, 2.4d. It vanishes 
for large enough k ^ Op^, leaving us with V{k) — 1, whereas for soft gluons one has 

V{k) oc K^al (2.35) 

One can say that the law (2.35) is driven by electromagnetic gauge invariance, which guarantees 
that long wave photons decouple from the charge neutral system. 









7 K 










a) 


b) 




Figure 2.4: The Fermi- Weizs acker- Williams diagrams for calculation of the flux of equivalent 
photons in positronium. 

Finally, recall that the derivation of the differential fiux of transverse polarized photons 
would equally hold if the photons were massive vector bosons interacting with the conserved 
current, the only change being in the propagator. For instance, for the charge neutral source 
one finds 

T^{xy,K') = i^^^\v{K) . (2.36) 

TT my J 
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Recall that the massive vector fields are Yukawa- Debye screened with the screening radius 



Rc^ . (2.37) 

my 

To the lowest order in QED perturbation theory the two exchanged photons in figs. 2. 3, 2.4 do 
not interact and we shall often refer to (2.36) as the Born approximation for the differential 
vector boson structure function. One can regard (2.36) as a minimal model for soft k behavior 
of differential structure function for Yukawa-Debye screened vector bosons. 

2.4.3 Differential gluon density in a proton 

The expression for the Born level unintegrated gluon density in color neutral proton can be 
obtained immediately by generalization of (2.33) and (2.36). The only thing one as to do is 
an accurate color algebra, which leads to 

j,(Born) ^ c!fN,^^V{k) . (2.38) 
71 

Therefore, a prescription how to include unintegrated gluon density is as follows: 

CfNc^^^V{h) = ^ T . (2.39) 

TT 

Note that in this prescription one of the two strong coupling constants in (2.21) (the lower 
one) is absorbed into the definition of J^. 

2.5 Final expressions 

With all pieces brought together, the answer for the total photoabsorption cross section reads 

n ^ Jo J J 
xj^ixg, k") [m^$g + [z' + (1 - zf]$l\ ; (2.40) 

aUx,Q') = ^Y^^U'dz [d'k f'^a,iq').Tix,,K')4Qh'il-zr^l, (2.41) 

TT , Jo J J 



where 



1 1 ^ k k — K 
$0 = ^ =^ ; ^1 = ^ ^ • 2.42 

Here 

= z{l - z)Q^ + m) , (2.43) 

and the density of gluons enters at 
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Here is the transverse mass of the produced qq pair in the photon-gluon fusion — > qq: 

„ ml + P ml + (k- k)'^ , , 

= — — + ^- . (2.45) 



z z 



No restrictions on the transverse momentum in the qq loop, k, and gluon momentum, 
K, are imposed in the above representations. The above used BFKL scheme defines DGSF 
uniquely in terms of physical observables. 

We note that equations obtained are for forward diagonal Compton scattering, but sim- 
ilar representation in terms of the unintcgratcd gluons structure function holds also for the 
off-forward Compton scattering at finite momentum transfer A, for off-diagonal Compton 
scattering when the virtualities of the initial and final state photons are different, Q'} ^ Q}-, 
including the timelike photons and vector mesons, = —i^v-i in the final state. 

The photoabsorption cross sections define the dimensionless structure functions 

Ft,lM') - ^—<^T,L (2.46) 
and F2 = Ft + which admit the famihar pQCD parton model interpretation 

Ft{x, Q') = Yl 4fe/(^' Q') + ^/(^' Q')] ' (2-47) 

f=u,d,s,c,b,.. 

where qf{x,Q'^),qf{x,Q^) are the integrated densities of quarks and antiquarks carrying the 
fraction x of the lightcone momentum of the target and transverse momenta < Q. 

2.6 The virtual Compton amplitude in the impact pa- 
rameter space 

A deep further insight into the Compton amplitude — and into the diffractive processes in 
general — can be gained by switching to the impact parameter (the transverse coordinate) 
representation. We will see that the answer will allow for a simple probabilistic form 

A = Jdzj ci¥(7(r)|*^(^,f)|^ (2.48) 

The quantity cr(r) has the meaning of the dipole cross section, while the 'if^{z,r) represents 
the photon wave function. 

In order to switch to the impact parameter space, we perform the 2-dimensional Fourier 
transform. We use the following equalities: 

^ = i-/d¥e-7o(r); Mr) ^ J d^k'^ ^ Mer) ; 

I d^re'^^y^ir) ; => h{r) = ( -i-^) [ d'kC^ 



-i-^] Ko{er) = -le-^K^ier) . (2.49) 
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This leads to representation 



(ffe'''^Ko{er) (l 



ie)-Ki{er) (l 
r ^ 



(2.50) 



Let us now take an approximation that there is no other k dependence in the photoab- 
sorption cross section. Then one will immediately have 



(frKl{er) 2[1 - cos(Kr)] . 



(2.51) 



Substituting these expressions into photoproduction cross sections, one gets 



TT 



dz d f 



mlKl^er) + [z" + (1 - zf]e^Kl{eT) 



—asJ^{xg, ff) 2[1 - cos(Kf)] ; 
^ V ef C dz I d'^f- AQh'^il - zf ■ Klier) 

71 "~ Jo J 



X 



-z^asT{xg^ K^) 2[1 — cos(/?r)] 



(2.52) 



2.6.1 Dipole cross section 



The above results can put into the form (2.48) by breaking Eqs. (2.52) into some positively de- 
fined cross section and the square of the photo wave function. This is done in an unambiguous 
way by defining the dipole cross section. 

Let us first consider the total quark-proton cross section 



'qp 



27r 



d^K 



(2.53) 



Note that this expression does not depend on the quark transverse momentum k. This means 
that precisely this cross section corresponds not only to the plane wave, but also to any 
transverse wave packet. In particular, a localized state in the impact parameter space (that is, 
a quark with a fixed separation from the proton) would be described by the same formula. 

Given this cross section, we can now ask for interaction of color dipole with the proton. In 
this each extra gluon that is attached to the antiquark rather then to the quark gives rise of 
extra phase factor exp(±irer) as well as extra minus sign. These factors for the four diagrams 
will then sum up to produce 

1 - e''^'' - e"^'^" + 1 = 2[1 - cos{Rr)] . 

The result for the dipole cross section is then 



I 



2n 



max 



di^ ( 
-^OLs I max 

AC \ 



AC 



A 



A 



T{x, k2)2[1 - cos(Kr)] 
J^(x,/?2)[1- Jo(Acr)]. 



(2.54) 
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Note that in contrast to (2.53) the argument of contains now the effect of a possible screening 
from the complementary charge. Indeed, even when is small, the strong coupling constant 
does not boost up, for such a boost requires the presence of soft gluon vertex correction loops, 
which are strongly suppressed by the color anti-charge of the antiquark. 

Having defined the dipole color cross section, we can now return to the photoabsorption 
cross section and cast it into form 

ar = j^dz j d^r adrp{r) ■ \^T^\^{z, f) , aL = dz J d^r ad^p{r) ■ \^^\^{z, r) ; 

r-)r = ^E^.' [rrAKlier) + [z^ + (1 - zfyK^^ier)] ; 

i 

= ^E^.' 4QV(1 - zfRUer) . (2.55) 

i 

Note that this representation literally represents the probabilistic form (2.6) of the forward 
scattering amplitude. 
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Chapter 3 

DGLAP vs. kt factorization 



The calculation of the forward Compton scattering amplitude conducted in the previous chap- 
ter within the framework of the /c^-factoriation approach can be used now to investigate the 
major similarities and the gross differences in comparison with the widely used DGLAP ap- 
proach to the computation of high energy reactions. 

It turns out that in the double logarithmic regime — that is, when both log | and log 
are large, we might expect that both approaches are applicable and their predictions should 
asymptotically converge. We are now going to demostrate that it is indeed the case, and 
during this analysis we will also show what sort of phase space restrictions DGLAP approach 
contains and what it leads to. 



3.1 How DGLAP and kt factorization approaches meet 
at high 

Recall the famihar DGLAP equation [18] for scaling violations at small x, 

dF2{x,Q^' 



dlogQ 



jjy[y^H^-ynG(^-,Q^^"-^G(2.,Q^)Y„ (3.1) 



where for the sake of simplicity we consider only light flavours. Upon integration we find 



FM') « EeJ /f ^^G(2.,C') , (3.2) 

f ' JiJ.^ Q 37r 



with //^ being the proper cut-off. In order to see the correspondence between the /cj-fa'Ctorization 
and DGLAP factorization it is instructive to follow the derivation of (3.2) from the k±_- 
representation. 

First, separate the K^-integration into the DGLAP part of the gluon phase space < 
Q = + k"^ and beyond-DGLAP region ■ One readily finds 



_o I3.dj 
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2 

Consider first the contribution from tlie DGLAP part of tlie pliase space <Q . Notice 
tliat because of tfie factor in (3.3), tlie straiglitforward integration of tlie DGLAP 

2 2 

component yields G{xg, Q ) and Q is precisely the pQCD hard scale for the gluonic transverse 
momentum scale: 

'^^ dK^ , -.0. I k k~K 



Jo K 



= y -qS ^ + =4 j W g ) (3.4) 

The contribution from the beyond-DGLAP region of the phase space can be evaluated in terms 
of J^{xg, Q ) and the rescaling factor C2: 

rdK^ / 2NW -.2N / 1 z{l-z)Q^ \ (1 z{l-z)Q^ \ ^2^_. ^2.^, ^2, 

)-^(2;fl,« ) (^=2 =4 j " ^ ) ^'^^ )^i^9,Q )HX9,Q ) 

= =8 ^-=6-^ — + =4 J ):^(2:s,<5 )logC2(a;s,Q ,z){3 

The latter form of (3.5) allows for convenient combination (3.4) and (3.5) rescaling the hard 
scale in the GSF 

G{xg,Q^) +J^(xg,Q'')\ogC2(xg,Q\z) ^ G(xg,C2(xg,Q\z)Q^) . (3.6) 

Here the exact value of I{xg, Q^) depends on the rate of the K^-rise of T[xg^ k^). At small Xg 

2 2 2 

and small to moderate Q one finds I{xg, Q ) substantially larger than 1 and C2{xg, Q , 2;) ^ 1, 

see more discussion below in section 9. 

— * 2 

Now change from dk"^ integration to dQ and again split the z,Q^ integration into the 
DGLAP part of the phase space <^ \Q^, where either zK^orl — z< and the 
beyond-DGLAP region ^jQ^, where < z < 1. As a result one finds 



30 Q2 



^2f^ ^^2,n4 o^n -.\n2 

+ -4 

Q Q Q 

if g'<ig2 



2A,^ - 2A^f + Ao^y if g'^ig^ 



(3.7) 



where ^ 

[z^ + (1 - zflz'^il - z)"" (3.8) 



/ dz[z 
Jo 



Let C2 be C2{xg, Q ,z) at a. mean point. Notice also that ~ g^, so that ~ 2x. Then 

2 

the contribution from the DGLAP phase space of Q can be cast in precisely the form (3.2) 
^^(^•.«^i...P -^ft"' f ^G(2..Q^) . (3.9) 
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The beyond-DGLAP region of the phase space gives the extra contribution of the form 

^ ) non-nOLAP ^ ^ 'f^^ ^ ^ ^ 

-E4^4^^(2^'^')- (3-10) 

Eqs.(3.9) and (3.10) immediately reveal the phenomenological consequences of lifting the 
DGLAP restrictions in the transverse momenta integration. Indeed, the DGLAP approach 
respects the following strict inequalities 

K'^<^k'^ and k"^ <^ . (3.11) 

2 

As we just saw, removing the first limitation effectively shifted the upper limit in the Q 
integral to ^Q"^ 7^ Q^, while lifting the second constraint led to an additional, purely non- 
DGLAP contribution. Although both of these corrections lack one leading log-Q^ factor they 
are numerically substantial. 

The above analysis suggests that the DGLAP and /cf-factorization schemes converge loga- 
rithmically at large Q^. However, in order to reproduce the result (3.9) and (3.10) for the full 
phase space by the conventional DGLAP contribution (3.2) from the restricted phase space 
(2.5) one has to ask for DGLAP gluon density Gpt{x, Q"^) larger than the integrated GSF in 
the /cj-factorization scheme and the difference may be quite substantial in the domain of strong 
scaling violations. 

3.2 The different evolution paths: soft-to-hard diffusion 
and vice versa 

The above discussion of the contributions to the total cross section from the DGLAP and 
non-DGLAP parts of the phase space can conveniently be cast in the form of the Huygens 
principle. To the standard DGLAP leading logQ^ (LLQ^) approximation one only considers 
the contribution from the restricted part of the available transverse phase space (2.5). The 
familiar Huygens principle for the homogeneous DGLAP LLQ^ evaluation of parton densities 
in the x-Q^ plane is illustrated in Fig. 3.1a: one starts with the boundary condition p{x, Ql) as 
a function of x at fixed Ql, the evolution paths {z, Q"^) for the calculation of p{x, Q"^) shown in 
Fig. 3.1a are confined to a rectangle x < z <1, Q\<Q^ <Q^-, the evolution is unidirectional 
in the sense that there is no feedback on the x-dependence of p(x, Q^) from the ,x-dependence 
of p[x^Q\) at Q\ > Q\. In Fig. 3.1a we show some examples of evolution paths which arc 
kinematically allowed but neglected in the DGLAP approximation. Starting with about flat 
or slowly rising G{x, Ql), one finds that the larger Q^, the steeper the small-x rise of G{x, Q'^). 

At X <S 1 the DGLAP contribution from the restricted transverse phase space (2.5) no 
longer dominates the multiparton production cross sections, the restriction (2.5) must be 
lifted and the contribution to the cross section from small and large > can no longer 
be neglected. The Huygens principle for the homogeneous BFKL evolution is illustrated in 
Fig. 3.1b: one starts with the boundary condition T{xq, Q'^) as a function of Q"^ at fixed Xq <S 1, 
the evolution paths (z, Q^) for the calculation of p{x, Q^) are confined to a stripe x < z < Xq, 



30 




-1 V -3 -2 -1 

10 1 ^10 10 10 1 



a) b) 

Figure 3.1: The Huygens principle for Q'^.x evolution of DIS structure functions with (a) 
DGLAP restricted transverse phase space and (b) for the BFKL x evolution without restric- 
tions on the transverse phase space and hard-to-soft & soft-to-hard diffusion. 

in contrast the the unidirectional DGLAP evolution one can say that under BFKL evolution 
the small-x behaviour of p{x, Q"^) at large Q"^ is fed partly by the x-dependence of soft p{x, Q^) 
at larger x and vice versa. The most dramatic consequence of this soft-to-hard and hard-to- 
soft diffusion, which can not be eliminated, is that at very small x the ,x-dependence of the 
gluon structure in the soft and hard regions will eventually be the same. The rate of such a 
hard-to-soft diffusion is evidently sensitive to the infrared regularization of pQCD, the model 
estimates show that in the HERA range of x it is very slow [22, 23]. 
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Part II 

Derivation of vector meson production 

amplitudes 
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Chapter 4 

Description of a vector meson 



In this chapter we first introduce the vector meson hght cone wave function (LCWF) and show 
how it emerges in diagrammatic calculations. Then, describing S and D wave type vector 
particles, we give at once expressions for S and D wave vector meson spinorial structures, 
which we then prove by computing the normalization condition for LCWF. Finally, we also 
calculate V — > e'^e" decay constants to be used afterwards. 

4.1 Bound states in QFT 

While describing particle motion in non-relativistic Quantum Mechanics, one usually deals 
with a configuration space particle wave function, which is a good description because the 
number of particles is conserved. So, when one has a system of particles and shows that the 
wave function corresponding to their relative motion descreases at large relative distances at 
least exponentially, one can speak of a bound state. 

In Quantum Field Theory (QFT) this approach needs an update, since the field function 
becomes an operator in Fock space. Besides, since a bound state always implies the presence 
of interaction, the projection of a physical bound state onto the Fock space of free, non- 
interacting, plane-wave state vectors has a rather complicated structure: 

\Vphys) = Co\qq) + Ci\qqg) + C2\qqgg) + c^lqqqq) + ... (4.1) 

We emphasize that in this decomposition quarks and gluons are assumed free, i.e. on mass 
shell. Coefficients c, can be called 'wave functions' of the given projection of a physical vector 
meson, with |cjp being the probability of finding a vector meson in a given state. 

The exact treatment of any reaction involving the vector meson must account for all terms 
in the above expansion. Demostration of a method that would account for all these terms is 
however still an unresolved task, and currently one is bound to the term-by-term analysis of 
the vector meson reactions. 

Given a large number of papers devoted to the high-energy reactions involving vector 
mesons, and in particular, the process of difractive vector meson production in DIS, one might 
expect that the lowest Fock state in the above decomposition has been already thoroughly 
studied. It turns out however that it is not so, for in all early calculations the importance of 
the spin-angular coupling inside the vector meson and dramatic effects it entails was heavily 
overlooked. 
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In this Chapter we close this gap. We construct a full and exact theory of the vector 
meson stucture, provided the vector meson Fock space is saturated only to the lowest qq state. 
Being only an approximation, this approach still is of vital importance to the whole field, 
for it results in a complete, self- consistent and self-contained spin-angular description of the 
vector meson. To our best knowlegde, our work represent the only satisfactory theory of spin 
structure of the vector meson. 



4.2 LCWF and vertex factor 

Let us now outline how a wave function of a bound state appears in the diagrammatic language. 

In the non-relativistic quantum mechanics, the two-particle bound state problem can be 
immediately reformulated as a problem for one particle of reduced mass /x, moving in the 
external potential. This reformulation allows one to split the wave function into two factors: 
the wave function of the motion of the composite particle as whole and the wave function 
corresponding to the internal motion of constituents. The former part factors out trivially, 
while the latter wave function obeys the following Schodinger equation 



ijj{r) = E'tPir) . (4.2) 



Since the wave function ^(r) and the interaction operator V{r) exhibit good infinity behavior, 
one can rewrite this equation in the momentum representation 

l^^(p) + J d^kV{k)ij{p -k) = Ei^ip) ; 

(p1 _ ^l;(p) = / d^kV{k)i;ip - k) . (4.3) 

In this notation, this equation can be viewed as a homogeneous non-relativistic Bcthc-Salpcter 
equation for the wave function '^(p) that describes the relative motion of constituents inside 
a composite particle. 
Let us now introduce 

r(p) ij{p) (4.4) 

Then Eq.(4.3) can be rewritten as 



E 



m 



This equation has an absolutely straightforward diagrammatic interpretation (Fig.4.1). One 
sees that T{p) stands for bound-state constituents transition vertex, with p being the relative 
momentum of the constituents, factor — E\ describes propagation of qq pair and V{k) 

stands for the interaction between constituents. Of course, the kinetic energy ^ E, the 

total energy, which is in fact negative, so no pole arises in the propagator. 

In the relativistic case, i.e. in QFT, it is not clear a priori whether the whole picture 
that involves somehow defined wave function and representation of the vector meson as free 
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Figure 4.1: The diagrammatic interpretation of the integral equation for vertex function r(p) 
(p is the relative constituents momentum). 

non-interacting constituents would work at all. So, in our approach we will be aiming at 
introducing an appropriately defined wave function and demonstrating that hard processes 
involving vector mesons can be expressed in terms of expectation values of qq amplitudes 
between wave functions, i.e. we intend to treat a hard process in a probabilistic, quantum 
mechanics-like manner. 

In the following we will show that this program succeeds. Namely, we will introduce the 
radial wave function of the qq state of a vector meson as 

^(g) = ^ J^^^ , (4.6) 

(the angular dependence of the wave function will be treated separately, see Sect.4.4) Here 
r(g) is the vertex factor, is the eigenvalue of the relativistic kinetic operator of the on 
mass shell qq state and niy is eigenvalue of the total relativisitic Hamiltonian, which is of 
course equal to the mass of the vector meson squared. Then, during an accurate and honest 
analysis of a hard process Feynman diagrams, we will always make sure that wave function 
(4.6) automatically appears in calculations and the rest looks the same as if both fermions 
were on mass shell. If we see that fermion virtualities modify the results, or if different Fock 
states mix during hard interactions of the vector meson, it would signal the invalidity of the 
free particle parton language and consequenly the breakdown of the whole approach. This 
restriction must always be taken into account when obtaining and interpreting the parton 
model-based results. 



4.3 Light cone formalism 

The term " light cone approach" to high-energy process calculations can have different mean- 
ings. Some prefer to re-formulate the whole QFT within light cone dynamics, introduce 
light cone quantization and derive light cone Feynman rules (on Light Cone Field Theory see 
[24, 25]). However, one should keep in mind that even within the usual QFT the light cone 
formalism can be freely used as a means to greatly simplify intermediate calculations. This is 
exactly the way we will use it. 

It was noted long ago [26] that the calculation of a high energy collision is simplied if one 
decomposes all momenta in terms of light cone nt. and transversal components, which we 
will always mark with the vector sign over a letter (so called Sudakov's decomposition): 

n+ = ;^(1, 0, 1) ; nf^. = -^(1, 0, -1) ; (n+n_) = 1 , (n+n+) = (n_n_) = 

if = p+< + p-n'^ +p^; = 2p+p_ - f . (4.7) 
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Indeed, imagine two high energy particles colhding with momenta and respectively 
and equal masses m. Then one can choose such a frame of reference that in the Sudakov's 
decomposition 

— q+n'^ + q-u't + q^ 

quantities p_ and q+ are large p_, q+ m while p+ — (p^ + w?)/{2p_) and q^ — -\- 
w?) I (2q'+) are small. The total energy squared of these two particles will be defined as 

s = 2g+p_ . (4.8) 

Note that our definition of s is somewhat different from the more familiar one (p + Q')^ by terms 
oc w?,p^. However, it is not of any imporatnce for us, since in the course of all calculations 
we will keep track only of the highest power s contributions, i.e. everything will be calculated 
in the leading power s approximation. 



Figure 4.2: Kinematics oiV qq vertex on the light cone. Vector meson momentum qy is 
taken incoming, constituents momenta are outgoing. 

Let us now go further and examine the kinematics of a typical qqV vertex (Fig.4.2). The 
Sudakov's decomposition of all momenta reads: 

q'^ = qv+n'+ + qv-nl ; 

k'i = ki+n'l + ki_nt + k^ ^ zqy+n'^ + yqv-u't + k" ; 

k^ = k2+n^ + k2-n^- - = (1 - z)qv+n^+ + (1 - y)qv-nt - k^ (4.9) 




so that 



ql = 2qy+qy_=ml- k'^ + k!^ = q>^ ; kf ^ m\ (4.10) 



i.e. quarks can be off mass shell. Now let us introduce momenta A;*, which would correspond 
to on mass shell fermions. The only component in ki subject to modification is k^^ component, 
or absolutely equivalently, the energy. Large ki^ components are insensitive to (reasonable) 
quark virtuality variations. So, to obtain the on mass shell momenta, one has to replace 

Then the 4- vector 

qi' = + kT (4.12) 
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squared is equal to 



= g2 = 2g+g_ = 2qy^ (kl_ + kl_) = . (4.13) 

And again we emphasize that the Feynman invariant mass (i.e. the total 4-momentum 
squared) of the virtual quark-antiquark pair is rriy. The quantity is the invariant mass 
of the free, non-interacting qq state (see 4.6). However, it is precisely M, not my that will 
govern the hard interaction of qq pair with gluons. 

Finally, it is useful to introduce the relative momentum of free qq system: 

= {K - k*), . (4.14) 

Then, trivial algebra leads to 

M2 = 4m^ + 4p2; ^ -p^ ; (pg) = , (4.15) 

where p is the 3-dimensional relative momentum in the qq pair rest frame of reference. Its 
components are 

P = {p,Pz); P = k; p, = ^{2z-l)M. (4.16) 



4.4 Spin structure of a vector particle 

Let us start with a well known example of a deuteron, which is a non-relativistic analogy of a 
vector meson: they arc both vector particles built up of two fermions. To have the correct P- 
parity, proton and neutron must sit in the spin-triplet state, thus leaving us with two possible 
values of their angular momenta: L = and 2. 

In the conventional non-relativistic language one describes the spin-angular coupling by 
the Clebsh-Gordan technique. The non-relativistic Feynman diagram calculations can be best 
performed in an alternative approach. Here a deuteron, being a vector particle, is described by 
a 3 dimensional polarization vector V. So, while calculating high energy processes involving 
d ^ pn transitions, one can use the following spin structure of deuteron-nucleon-nucleon 
vertex: 

0:r0,-V (4.17) 

Since both nucleons can be treated on mass shell, only two terms enter Fj, which can be 
written as: 

0n [u{p)cr' + w{p){3py - (5V)c^^] (f>p ■ V' (4.18) 

Here cr' are Pauli matrices and p is the relative proton-neutron momentum. One immediately 
recognizes here spin structures corresponding to pn pair sitting in S and D waves respectively. 
In particular, squaring the above expression gives 

(VV*) for 

3(pV)(pV*) - (VV*)p2 for SD interference 

3p^(pV)(pV*) + (VV*)p^ for |L>p (4.19) 
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Now, let us go relativistic and turn to vector mesons. The polarization state of a vector 
particle is described by a four- vector V^. Therefore, a general form of qqV vertex has the form 

uT^u ■ ■ r(p) , 

where T{p) is the familiar vertex factor. Up to now, it has been customary in literature to 
choose the simplest form of the spinorial structure F^: 

u'^^u ■ ■ r(p) . (4.20) 

However, one must admit that (4.20) is simply an analogy of qq^ vertex and does not reflect 
the true internal structure of a vector meson. It is known [27] that the correct spinorial 
structure corresponding to pure the S wave qq state reads 

^^ = ^^- MTk = - m(M + 2m) ) " ^^^^^ " ^^'^^^ 

It is implied here that spinorial structures are inserted between on mass shell spinors in 
accordance with our principal guideline (see discussion in Sect. 4. 2). 

Once S wave spinorial structure is established, the expression for D wave can be obtained 
by contracting 5* wave with the symmetric traceless tensor of rank two SpiPj — Sijp^, rewritten 
in the Lorenz notation. To do so, one should replace 

■ A - M Mil 

Pi ~^ P/i ) 'Jij ~^ Qliv + 

(in the qq pair rest frame of reference q^ = (M, 0, 0, 0)). However, since inserted between 
on mass shell spinors gives zero due to the Ward identity, one obtains the required tensor in 
the form 3pi^p„ + g^uP^- Its contraction with yields 

n, ^ 2 f 2 M + m \ 
= {^Pi^Pu + Qi^uP ) ■ o^p7p = P 7m + + "^)Pm = [P 9hu^ Pi^Pu j Iv = T)^^-/,, . 

(4.22) 

We will prove below that structures (4.21), (4.22) after being squared indeed perfectly repro- 
duce (4.19), i.e. they indeed correspond to pure 5" and D waves. 



The quantities Sy^y and T>^y used in (4.21), (4.22) have the meaning of SjD wave projectors, 
which will be used in all subsequent calculations. Namely, all calculations will be at first 
performed for the naive qqV vertex (4.20) and then we will apply the projector technique to 
obtain expressions for S and D wave states. 



4.5 Vector meson LCWF normalization 

Before tackling the diffractivc vector meson production process, we first should have a pre- 
scription of normalization of the vector meson wave function. 
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1v TT 1v 



Figure 4.3: Diagram used for normalizing the vector meson LCWF. The amphtude of this 
diagrams is set equal to 2q^i. 

4.5.1 Naive qqV vertex 

A natural way to normalize the wave function of a composite system is to put the amplitude 
given by this diagram in Fig.4.3 equal to 2gf+i. Here extra leg carries zero momentum but 
couples to the fermion line as 7^n^!. Note that for a charged composite particle (a deuteron) 
this is precisely setting the electric formf actor equal to unity in the soft photon limit. 

As described above, we first treat qqV vertex as u''^^u ■ T{p). In this case the general 
expression for this amplitude is 

^ _ (-1) f ^ Sp{iViT ■ i{k - qv + m) ■ iV^'^* ■ i{k + m) ■ ih_ ■ i{k + m)] 



I 



(27r)^ J [fc^ — + ie] ■ [k"^ — m? + ie\ ■ [{k — qy^ — m? + ie] 



(27r)4 J [k'^ - ief[{k - qvY - -m? + ie\ 

where A^c = 3 is a trivial color factor originating from the quark loop. We deliberately 
recognized Vi and V2 as distinct entities just to make sure later that such a loop is indeed 
diagonal in polarization states. 

The first step is to rewrite this expression in terms of Sudakov's variables. As usual, one 
implements decomposition 

k" = zqv+n"^ + yqv-n''_ + k^ ; q^ ^ 2qv+qv- = rn^ 

and transforms ^ 

d'^k = -mycfkdydz . 

Now we note that vertex functions F do not depend on y (and neither does the trace, as will 
be shown later), so we can immediately perform the integrations over y by means of Cauchy 
theorem. Indeed, since the integral 

dy =^ ^ , (4.24) 

-00 [yzruy — {k"^ + m?) + ie]^ [(1 — y){l — z)my — {k^ + m?) + ie\ 

is convergent and has good infinity behavior, one can close the integration contour in the most 
convenient way. To do so, one should analyze the position of all poles on the complex y plane: 



P + ie P + 

yi,2 = ——2 —2-; ys = 1 - — 3—2- + 



zruy zruy ' {i — zjrriy [i — zjuiy 



{\-z)m\ (l-z)7 
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One sees that if 2; < or 2; > 1, a// poles lie on the same side of the real axis in the complex 
y plane, which leads to zero contribution. The contribution that survives comes from region 
< 2; < 1, which has in fact a simple physical meaning: all constituents must move in the 
same direction. In this region, we close the integration contour through the upper half-plane 
and take a residue at y = yz- Physically, it corresponds to putting the antiquark on mass 
shell. After this procedure, one gets for (4.24): 



2-Ki 1 



2Tii (1 



(1 - z)ml [yzml - {k^ + m?) + leY y=yz (1 - ^Wv [k'^ + m?- z{l - z)mlf 
_ 2m I 1 
~ ~ml z^{l-z) [M'^-mlf ' 

One immediately recognizes here the same two-particle propagator as in (4.6). Therefore, the 
equation for the amplitude reads 

We now turn to the trace calculations. As we found during the calculation of the forward 
compton scattering amplitude, all the fermions can be treated in the trace calculations as if 
they were on mass shell. We will use this property in all subsequent calculations. 

The easiest way to calculate the traces in our case is to do it covariantly, without involving 
further the Sudakov technique. Since quarks in numerator can be treated on mass shell, we 
first note that 

{k + m)n^{k + m) ^ 2(k*n^)(k* + m) 

so that 

-\sp{...] = -zq+Sp{ViCk* -q + m)V*Ck* + m)} = -2zq+ [M^ViV^) + 4{Vip){V*p) 

where p is the relative quark- antiquark momentum [see (4.15)]. Note that in the antiquark 
propagator we replaced k — qy — > k* — q, since the antiquark is now put on mass shell. Besides, 
we explicitly used here gauge condition (qV) = 0, which means that polarization vectors must 
be written for on mass shell qq pair, not the vector meson, — another important consequence 
of our approach. Substituting this into (4.25), one gets 

^ = (2# y ^'^1(1^'^'' [-M\V,V*) - A{V,p){V*p)] . (4.26) 

A prominent feature of this equation is the orthogonality of Vl and Vr polarization states — 
the necessary condition for any normalization prescription. 

The next step is to realize that the integral can be cast in the form of d^p integration by 
means of 

d k ^ -Tjdp^d^p = —d^-p . 



z[l-z) M " M 

Tus, the final expression for normalization condition is 

Nr. /■ „ 4 



1= " 



(27r)3 / d'^M\'l^?[-M\V,V*) - A{V,p){V*p)\ . (4.27) 

We see that the expression being integrated is explicitly spherically non-symmetric, which is 
a manifestation of a certain D wave admixture. Thus we now apply projector technique to 
obtain results for S and D wave states. 
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4.5.2 Normalization for S wave vector meson 

The correct expressions for pure S/D type vertices can be readily obtained with the aid of 
projector technique. Namely, to obtain an expression for S wave, replace 



Such a replacement for Vi leads to 



-M'^{Viv;) - A{Vip){v;p) ^ -m^{Viv;) - 

Then, one applies the same replacement to V2 to obtain 

AM 



AM 



M + 2m 



{Vip){V*p) 



-M'{ViV^ 



M + 2m 

Therefore, the answer for S wave states reads 



(V^p){V*p)^-M'{V,V*) 




which is manifestly spherically symmetric. 



(4.28) 



(4.29) 



4.5.3 Normalization for D wave vector meson 

Results for D wave states are derived in the same way. The replacements Vfj, 

-M^p\ViV;) + 3MV(Kp)(^*p) . 

After angular averaging 

one gets the normalization formula for D wave state: 




V,yV^^ lead to 



(4.30) 



(4.31) 

Several remarks are in order. First, S wave — > D wave transitions are forbidden. Indeed, 
such an amplitude will be proportional to 



M'[p\V,V*) + 3iV,p){V*p)]. 



(4.32) 



which vanishes after angular integration. Then, we emphasize that the structure of results 
(4.29), (4.30), (4.32) is absolutely identical to Eq.(4.19). This fact can be viewed as the proof 
that spinorial structures (4.21), (4.22) indeed correspond to pure S and D wave states. 
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a) 



S) 



Figure 4.4: Normalizing LCWF to V{V — > e+e ) decay width: (a) the diagram for F — > 7* 
transition, (b) the diagram for V — > e^e" decay. 



4.6 Decay constant 

An additional normalization condition consists in relating the vector meson wave function to 
the experimentally measurable physical quantity — F — > e^e~ decay width (Fig.4.4). The 
loop at Fig. 4. 4a describes transition y — > 7* and enters the amplitude of the decay V e^e~ 
(Fig.4.4b). Let us define the decay constant via relation 



(4.33) 



So defined fy has dimension dim[/y] = m^. The quantity Cy refiects the fiavor content of a 
vector meson (in the previous calculations it simply gave unity) and is equal to 



1 



1 2 



Cv = 



V2'3V2' 3' 3 



(4.34) 



for p, cu, (f), J/i/j mesons correspondingly. 

Knowing that such a loop does not mix polarization states, we can multiply both sides of 
eq.(4.33) by V* and get the expression 



V 



Sp{iV* ■ i{k + m) ■ iVT ■ i{k - qv + m)} 
{k"^ — m? + ie) ■ {{k — qvY ~ + i^) 



(4.35) 



Calculations similar to the above normalization condition derivation yield (for the naive 
type of vertex) 



/i 



Nr. 



V 



(27r) 



dz 



z 1 - z 



-d^k tlJv[-M^{VV*) - 4:{Vp){V*p)] 



(4.36) 



Applying now the projector technique, one gets in the case of S states (after proper angular 
averaging) 



(4.37) 



and in the case of D wave states 




(4.38) 
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Finally, one can write down the expression for the decay width in terms of fv- 



r{V e+e") 



1 



my 



4n\A\'^ 




(4.39) 



327r2m?; 



2 



3m?, 



This formula can be used to extract the numerical value of jy from experimental data. 



Later on, we will be presenting numerical analyses of vector meson production cross sections, 
for which we will need some wave function Ansatz. Here we describe two forms of the wave 
function that will be exploited there. By no means should they be expected to accurately 
represent the true radial wave functions in a vector meson. Our Ansatze are pure guesses, 
based on non-relativistic quantum mechanical experience, of how the wave function might look 
like. Undoubtedly, such an approach involved a certain degree of ambiguity to the numerical 
results, and in our subsequent analysis we will study this ambiguity in detail. 

4.7.1 Suppressed Coulomb wave functions 

A first guess for the vector meson wave function, especially in the case of heavy non-relativistic 
meson, would be a Coulomb-like form, similar to wave function of a positronium: 



where a is a typical size of the meson. 

However, such a hard wave function will not fit our course of calculations, since the expres- 
sion for the decay constant will be ultra-violet divergent. Furthermore, as wc will sec later, 
this hard wave function will lead to the vector meson production amplitudes saturated not at 
the scale < 1/a^, but will extend to 1/a^ <^ ^ Q^- 

Thus, it appears that the hard wave function Ansatz leads to complications, which do not 
seem to be resolvable within the lowest Fock state only. Therefore, starting from now, we will 
limit ourselves to the soft wave function Ansatz only, "soft wave function" meaning that all 
integrals of physically relevant amplitudes involving the wave function will be saturated by 



If we still prefer to have a Coulomb-like wave function, we can consider its slightly regu- 
larized form, which we will call the "suppressed Coulomb" wave function. Besides, in order 
to be able to conduct simple estimates, we will take as simple form as possible. So, in this 
ansatz the normalized wave functions read 



4.7 Ansatz for LCWF 



1 



(4.40) 



(l + p2a2)2 ' 




1 



V'25(P^) 




C2 i^node - ajP^) 

/M (l + aip2)3 ' 




1 



(4.41) 
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with normalization constants to be determined from Eqs.(4.29) and (4.31). Here parameters 
are connected to the size of a bound system: in the coordinate representation ip^ oc exp(— r/oj). 
For strict Coulomb functions one would have = 802/2 = 3ai = Rsohr, where Rsohr is the 
Bohr radius. However, this relation should be treated with care in the case of qq quarkonia, 
where the quark- antiquark potential is quite complicated and therefore a, should rather be 
considered as free parameters. Value of ^node pinpoints the position of the node in the 25* 
radial wave function. For the pure Coulombic system ^node = 1; but in our case the exact 
value of inode should be obtained from the requirement of the orthogognality between IS and 
25" states with ai and 02 fixed from other requirements. 



with a again being a typical size of the wave function. This wave function Ansatz corresponds 
to the case of a strong confinement. Although the approximately quadratic potential that 
leads to such an abrupt dcscrcasc at > is not exactly what is suspected about the 

color-singlet static quark- antiquark potential (in the quenched approximation), these wave 
function still possess the main confinement-like properties. 
In this ansatz one has 



Note again that for purely oscillator potential one also has relation an = 0,2 = «i, which might 
not hold in our case, since the oscillator type potential is also a crude approximation of the 
true quark-antiquark interaction. The position of node ^node would be equal to 3/2 for pure 
oscillator model, but again in our case its value can turn out different. 



4.7.2 Oscillator type LCWF 

By oscillator-type wave function we mean 



■0(p) OC exp 





(4.42) 
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Chapter 5 

Vector meson production amplitudes 

5.1 Preliminary notes 




p p p p 

a) b) 



Figure 5.1: The QCD-inspired diagrams for 7*p Vp process with two gluon i-channel. 
Only Diagr.(a) does contribute to the imaginary part of amphtudes. 

Having set up the notation and defined and described a vector meson by itself, we are now 
ready to evaluate the full set of amplitudes of its off-forward virtual diffractive photoproduc- 
tion. 

In the pQCD motivated approach to this process the pomeron exchange is viewed as a 
two-gluon exchange as it is shown in Fig. 5. la. Using the scalarization procedure, we will split 
the diagram into 2 pieces and will treat each of them separately. The upper blob describes 
the pomeron-assisted transition of the virtual photon into a vector meson. In the perturbative 

2 

QCD approach, which is legitimate here due to the presence of the relevant hard scale Q = 
+ z{l — z)Q^ , the qq fluctuation of the virtual photon interacts with two hard gluons and 
then fuses to produce a vector meson. This interaction is described by four diagrams given in 
Fig. 5. 2, with all possible two-gluons attachments to qq pair taken into account. All of them 
are equally important and needed for maintaining gauge-invariance and color transparency. 
The latter property means that in the case of very soft gluons the upper blob must yield zero, 
for the qq pair is colorless. 

The lower part of the general diagram Fig. 5. la is of course not computable in pQCD. The 
physically meaningful procedure is to relate it to the experimentally measurable gluon density. 
To do so, we will first calculate this lower blob in the Born approximation and then give a 
prescription how to introduce the unintergrated gluon density. In the course of this procedure. 
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we will neglect in the intermediate calculations proton off- forwardness and take it into account 
only at the very end, as a certain factor to the unintegrated gluon density. 



I I 

f ^ 

8\ ^ 



a) 
V 



c) 



t 



b) 



V 



d) 



Figure 5.2: The content of the upper blob in Fig.S.la in the pQCD approach. The true vector 
meson internal structure is approximated by qq Fock state. 




Figure 5.3: (a): A particularly useful convention of the loop momenta (only transverse com- 
ponents of quark momenta are shown), (b): the propagator notation used while calculating 
denominators. The crosses denote on mass shell particles. 



5.2 Notation and helicity amplitudes 

In our calculation we will use the following Sudakov's decomposition (see also Fig. 5. 3a) 

kp, = yPfl + zq^j! + /c^ ; 

= ap^ + (3q^ + ; 
A;. = + aq^ + (5.1) 

Here k and k are momenta that circulate in the quark and gluon loops respectively; A is the 
momentum transfer. Vectors and q^ denote the light-cone momenta: they are related to 
the proton and virtual photon momenta as 

Pm=p/ + — % = <lj-xpj-, q'^=p'^ = 0; = s = 2{p'q') . (5.2) 
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As it was mentioned in the Introduction, the condition x ^ 1 is necessary to speak about 
diffractive processes. The longitudinal momentum transfer can be readily found from kine- 
matics (see Fig. 5.1). To the higher power s terms it reads 

ml = p^ = {p- Af = ml - 2{pA) + A^; cr = — — ; 

m'y^{q + Af^-Q' + 2{qA) + A' ^ 5 = x + ^^^^^ . (5.3) 

s 

The final vector meson momentum reads: 

QVf. = (Itl + — P/x' + . (5.4) 

s 

— * 

Finally, throughout the text transverse momenta will be marked by the vector sign as k and 
3D vectors will be written in bold. 

There are several possible helicity amplitudes in the transition — > V\y. First of all, 
both photon and vector meson can be transversely polarized. The polarization vectors are 

eTM^e-;; = + (p' - g% . (5.5) 

s 

Note that we took into account the fact that the vector meson momentum has finite transverse 
component A. Then, the virtual photon can have the scalar polarization (which is often called 
longitudinal; we will use both terms) with polarization vector 

eo^ = ^(?' + V)/.- (5.6) 

Finally, the longitudinal polarization state of a vector meson is described by 

1 / A^ - 

Vl, = ^ U' + p; + A, I . (5.7) 

Note that, as we already mentioned, in the self-consistent approach we must take the running 
polarization vector for the longitudinal polarization state. It depends on M, not my, which 
reflects the fact that in our approach we first calculate the production of an on-shell qq pair 
with (whose dynamics is governed by M) and then projects it onto the physical vector meson. 
We stress that this projection will automatically arise in the course of usual Feynman diagram 
evaluation. 

Thus, there are 5 different amphtudes: 

L 

T^T {X^^ Xv) 
L 

L^T 

T^T {X^^ -Xv) (5.8) 
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The first two are helicity conserving amplitudes. They are dominant and almost insensitive 
to the momentum transfer A. The next two are single helicity flipping amplitudes. They are 
unavoidably proportional to | A| in the combination (eA) or (V^*A) and would be vanishing for 
the strictly forward scattering. Finally, the last amplitude corresponds to the double helicity 
flip and will be proportional to (eA) {V*A) . 



5.3 General amplitude 

We will take diagr.(c) at Fig. 5. 2 (it is shown in Fig.5.3) as a generic diagram and perform 
a thorough analysis for it. It turns out that the other diagrams are calculated in the same 
fashion. 

The general expression for the amplitude given by diagr.(c) reads: 

. , f d^k f df^K . , . „/ ,,/, p — + mi, , „/ 4/, 



(27r)4y (27r)4 " ' ^ [(^ _ _ ^2 + 

S-p \iee i(^4 + m) {—igYt^) «(^3 + m)iV* i(^2 + rn) {—ig^H^) i{ki + m)| 
[ki —w?-\- ie\ [k2 —w?-\- ie\ [/cg — w? -\- ie\ [k1 — w? -\- ie\ 

(5.9) 

Here cy is the same as in (4.34) and F is the famihar qq vertex function. Note that we 
introduced 'gluon mass' // in gluon propagators to account for confinement at a phenomeno- 
logical level. 

Let's first calculate the numerator. 



5.4 Color factor 

If we consider strictly forward gluon scattering off a single quark, we have 

1 1111 N'^ — 1 1 9 

j^Sp{t^'t^'} ■ 5AA'5BB'Sp{t^t^] = ^Jab-5ab = 2^]^ = 2^^ ^ 3 ^^'^^^ 

However, we should take into account that quarks are sitting inside a colorless proton, whose 
color structure is 

i^coiar = ^e^'^q^'q^q" (5.11) 

In this case there are two ways a pair of gluons can couple 3 quark lines (see Fig. 5. 4). In the 
first way both gluons couple to the same quark. Since the quark momentum does not change 
after these two interactions, the nucleon stays in the same state: {N\N) — 1. In the second 
case gluon legs are attached to different quark lines, so that extra momentum k circulates 
between quarks, which gives rise to the factor (iV| exp(mri — ii<ir2)\N), i.e. to the two-body 
formf actor. Therefore, for the lower line instead of 

^Mt-t^) = (5.12) 
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one has 

= Sp{t^t^} - Sp{t^t^}{N\ exp(mri - iKr2)\N) 

= ^6AB{l-{N\exp{iKn-iKr2)\N)). (5.13) 

Note also that a similar calculation for Nc number of colors would yield the same result. Thus, 
the overall color factor is 

^CfN,V{k) = 2V{k) = ^CfN,{1 - {N\ exp(mri - iKr2)\N)). (5.14) 



K 



K ' i K 



N 



N 



N 



N 



X Nc (Nc-1) 



Figure 5.4: The ways two glons can couple a colorless nucleon. 



5.5 Scalarization of upper and lower parts 

As known, the highest power s contribution comes from so-called nonsense components of 
gluon propagator (density matrix) decomposition: 

2p'g', 2p'q', I 2p'q', 
9nti' - — \ Hf^^' —■ I^-Loj 

The lower (proton) line gives then 

u'{p) ■ q'{p — ki + mp)q' ■ u{p — A) (5.16) 

To the highest power s order, it can be rewritten as 

= \sp{p'm}^s'. (5.17) 

forward ^ 



u'^q^p'q^Up 



As we mentioned, the effect of off-forwardness (skewedness) will be taken into account later. 
So, combining all factors, one has for numerator of Eq.(5.9) 

(47ra,)2 V47raem Cy ■ -CfN^V{k)—s^ ■ Sp le {k^ + m) q' (/cg + m) V* (fcg + m) q' {h + m)| 
= {inasf V^^cv2CfNcV{k)-2s^ -^X-f* ^V). (5.18) 
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Note that wc factored out 2,s^ from the trace because it will appear later in all trace calcula- 
tions. So, the resulting expression for amplitude (5.9) looks like 



(47ra,)V(«) 



[ki — m? + ie] [/cl — ra^ + ie] [k"^ — m? + ie] [kl — m? + ie 



{p — Ki)'^ — nip + ie [ki — //^ + ie] [k^ — //^ + ie] 

(5.19) 



One can now immediately write similar expressions for the other three diagrams (Fig. 5. 2 
a,b,d). Indeed, they will differ from Eq.(5.19) only by the last line: they will have different 
expressions for traces and propagator structures (i.e. the exact values for ki.). 



5.6 Denominator evaluation 

Now we turn to the calculation of denominators. As usual, we implement Sudakov's decom- 
position (5.1) and make use of relation 

d k — -sdydzd k . 

The physical picture of the way we will do the resulting integrals is the following. We are 
interested only in the imaginary part of these diagrams. In fact, it can be shown that at the 
level of accuracy used here the diagram in Fig. 5. la gives rise only to the imaginary part of the 
amplitude. The real part is given by Fig. 5. lb and can be readily found from analiticity (so 
that here is no need for additional calculations), however it turns out small due to smallness 
of pomeron intercept, so we will neglect it in our subsequent calculations. 

The imaginary part is computed by setting three particles in the s-channel cut on mass 
shell (this is illustrated in Fig. 5. 3b). One way to do so is to apply Cutkosky rule to modify 
our expression. Another, more straightforward and 'honest' way is to calculate three of the 
integrals (namely, over a,f3,y) via residues. That's what we are going to do. 

The details of this calculations are given in Appendix A. Here we cite the result: 



Im j dy dz da d(3 



r 



[all propagators] 
7ri\ / 27ri\^ r dz 



■ —, TilJviz.k') ■ — — ^(5.20) 

S)\ s) J Zil-zf"^' ' [;tf+m2 + ^(l-z)Q2](^2 + ^2)2^ > 

Here ki is the transverse momentum flowing through photon vertex along the fermion line. 

— * — * — * 

Particularly, for diagr.(c) it is equal to ki — k — {1 — z)/S. — K2 (with the specific quark loop 
momentum choice given at Fig.(5.3a)). 

Thus, the amphtude for diagr.(c) has the form 

11 f Tii\ f 27ri\2 1 



A = V^Ko^ ACfN^s^ Cv -s- (-y) ( — j^) 



(27r)8 

/ -r^^d'kMz^k') [ (4vraJ' • — • (5-21) 

J Z(l-Z) ^"^^ ' ^y(«2 + ;,2)2^ [^2+^2 + ^(l_^)Q2] ^ ^ 
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After bringing all coefficients together, one gets 



/(c) 



A;? + m2 + z(l - ^)g2 



(5.22) 



The other diagrams are calculated in the same way. The most important difference is that 
for each diagram we will have a propagator + z{l — z)Q^] with its own definition 

— * 

of ki, the transverse momentum in photon vertex: 



diagr.a kia = k — {1 — z)A = f — -A 

— * — * — * X ~* 

diagr.b k\\, — k — {\ — 2;) A + + -A = r + k 

diagr.c A^ic = — (1 — ^)A — + -A = r — k 

— * — * — * X — * 

diagr.d k]_^ — k^ ziS. — r + -A 



(5.23) 



Here r = ^— (1 — 2 2;) A/2. Thus, the whole expression for the imaginary part of the amplitude 
is 



A — IS— 



647r5 



X 



1-z 



/(«) 



Jib) 



kl + m2 + z{l - z)Q'^ kl^ + m2 + z{l - z)Q^ 



+-. 



+ 



lid) 



kl + m2 + z{l - z)Q'^ 1- zkla + m'^ + z{l - z)Q^ 



(5.24) 



5.7 Off- forward gluon density 

In Sect. 2.4 we discussed the forward unintegrated gluon density J-{xg, k) and developed a 
prescription (see Eq.(2.39)) how to introduce it into the fc^-factorization calculations. Being 
devised for forward scattering processes only, this gluon density bears a clear probabilistic 
sense, which is in fact reflected in the word 'density'. 

In the present case of vector meson production, the initial and final states are kinematically 
distinct, therefore the forward unintegrated gluon density in its pure form is not the relevant 
quantity; instead, off-forward (or skewed) gluon structure function [28] 

J^{xi,X2,K,A) (5.25) 

should be used. It depends on the lightcone momenta xi and X2 carried by the first and 
the second gluon, on the transverse momentum k inside the gluon loop, and on the total 
transverse momentum transfer A. At a; 1 = 0:2 and A = 0, the forward gluon structure 
function is recovered, which means that even the strictly forward vector meson production 
should be described by off-forward gluon structure function. 
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The experimental determination of the off-forward gluon density might be in principle 
possible from accurate measurements of the off-forward virtual Compton scattering process, 
but the lack of such measurements makes this analysis not feasible in the nearest future. Thus, 
the original idea to involve in the vector meson calculations as little ambiguity as possible by 
determining the gluon content of the proton from other reactions does not work. One might 
of course expect that replacement 

^(xi, A) ^ ^ f d] (5.26) 



should be quite legitimate in kinematical regimes when xi ~ 2:2, A <^ However, such a 
regime takes place only in the photoproduction of light mesons and is severely violated as we 
go to higher or higher my, therefore such a replacement would be a poor option for the 
most of the case we study. 

There is however way around, which allows us to reduce the unknown off-forward gluon 
distributions to the forward ones. As shown in [29] , if the energy behavior of the gluon densities 
is describable by a simple Regge-type behavior 

^oc(i)' (5.27) 

then in the case Xi ^ X2 the off-forward unintegrated gluon structure functions can be related 
to the forward unintegrated gluon density according to 

T{xu a;2 ^ 0, A ^ 0) = i?, • ^(xi, ti)- Rg= p ■ (5-28) 

Bearing this exact result in mind, one can hope that a similar relation will hold for gluon 
densities with a somewhat more complicated energy behavior, if the effective exponent A is 
calculated in the vicinity of the kinematical point 

Note that for the purposes of approximate numerical calculation correspondence (5.28) can 
be further simplified. Introduce an agrument shift c(A) such that 

22A+3 rcA + ^) / 1 \ ^ / 1 \^ 

J^(xi,0,«,0) oc — = (5.29) 

^ ' ' ' ^ r(A + 4) \xiJ \c{X)xiJ ^ ^ 

holds for all Xi. Simple arithmetics shows that c(A) changes from ?a 0.435 at A = to 0.4 at 
A = 1. Given this very fiat dependence, we can approximate c(A) by a constant value 0.41 so 
that 

T{xg, 0, K, 0) ^ J^{OAlXg, k) . (5.30) 

This form will be used in our numerical calculations. 

The effect of non-zero A comes both from the Pomeron-exchange and from the proton 
impact factor. Since Pomeron singularity moves in complex j plane as |t| = A^ changes, the 
value of the effective Pomeron intercept will be affected by A^, a customary representation of 
this effect (in the linear Rogge trajectory approximation) being 

aip(A2) =aip(0) -6]pA2. (5.31) 
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Although in our case the effective intercept of the ghion density is not an input number, but is 
generated dynamically, we still account for the above effect by multiplying the gluon density 
by factor 

r — * Q 1 

exp — 6ipA log(xo/x) , (5.32) 

with xo = 0.03 and trajectory slope 6ip being different for hard {x^^'^ = 0.07) and soft 
{x^^^ = 0.15) parts of the gluon density. Since the resulting contribution to the slope increases 
with energy growth, this effect is called the diffractive cone shrinkage. 

The second — and the most significant effect of non-zero A — comes from the proton 
impact factor. Effectively, it amounts to introduction of a proton formfactor F(A^), which 
would be equal to 1 at A^ = and would start significantly decreasing when A^ grows larger 
than a certain scale l/A^. In our numerical calculations we used the dipole formfactor 

F(A^) = I 1 I (5.33) 

with A2 = 1 GeV-2. 



5.8 Final results for the naive vertex 

Now, with the off-forward gluon structure function properly defined, the only thing left to 
be computed is integrands P{'y — > V). For convenience, their calculation is also given in 
Appendix B. It turns out that the results can be written in the same way for all four diagrams 
via ki given by (5.23), i.e. all quantities: 



1-z 



j{a) ^ j{b) ^ j{c) ^ 



1-Z 



can be written in a similar way 
T 



{eV*){m^ + kki) + {V*k){eki){l - 2zf - {ek){V*ki) 
L^L -Az^{l - zfQM 
T^L 2z{l- z)M{eki){l-2z) 
L^T -2z{l - z)Q{l - 2z){V*k) 

— * 

Therefore, we can cast amplitude (5.24) in a compact form with aid of functions $i and 



$9 



and 



$1 



+ 



(f+/?)2 + g' (f-/?)2 + g' (f+A/2)2 + g' (f-A/2)2 + Q' 



(5.34) 



— ♦ I — ♦ 

r + K, 



r — K 



+ 



f+ A/2 



+ 



f- A/2 



(f+«)2 + g' (f_^)2 + g2 (^+^/2)2 + q' (f-A/2)2 + g' 



(5.35) 
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With these functions, for the naive qqV vertex the whole expression in square brackets in 
(5.24) wit sign minus (which we denote here as Ix^\y) has the form: 



III 
Itt 
Itl 
Ilt 



-AQMz^{l-z)''^2; 

(eV*)[m'$2 + (^^i)] + (1 - 2zf{kV*){Si) - {ek){V*$i) ; 
2Mz{l- z){l-2z)(e$i); 
-2Qz(l - z)(l - 2z)(Vk)^2 ■ 



(5.36) 



Replacing V{i-i) by the unintegrated gluon density according to (2.39) and using approxi- 
mation for the off-forward gluon density (5.30), we obtain the final expressions for a general 
amplitude of reaction 7^ — > Vxy in the naive vertex: 



A{x,Q\A) = -is 



47r2 

'^asJ^ix, K, A) ■ /(7* V) , 



(5.37) 



5.9 Final results for S and D wave amplitudes 

Now we can use projector technique to obtain results for S/D wave states. 

^P,Pu _ rS_r , 2(Vp) 



m{M + 2m) 



1 + 



m{M + 2m) 



Pn ^ In ; 



(5.38) 



T/ T/T) -T) 2 , {M + m)p^p^ 



m 



tD t 2 (M + m)(Vp) 
m 



Note that (Vp) is 3D scalar product. While contracting, we encounter terms proportional to 
Pii ® In which should be understood as 



Pn®ln^ IvAy ^ P} + IvA'^ = V,^p, = -{2z - 1)M}. 



The result of this substitution reads: 
for ct' 



It-.t{V ^p} + It-.l{1 = 14 ^ p. = -i2z - 1)M} 



m 



(e7c)$2- (e$i)(l-2z)' 



for cq: 



Il^t{V ^p} + Il^l{1 = V,^p, = -{2z- 1)M} 
= -2Qz{l - z){2z - l)m2$2 

So, the resulting integrands for 5" wave type mesons are 



(5.39) 



(5.40) 



(5.41) 
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-AQMz'^il - zf 



1 + 



(1 - 2zf 2m 
Az{l - z)M + 2m 



$2; 



= (eV*)[m'$2 + (k^i)] + (1 - 2zf{kV*){e^i) 



M 



M + 2m 



2m 



M + 2m 



2Mz{l-z){l-2z){e<^i) 



1 + 



(1 - 2zy 2m 
4z{l - z)M+2m 



Mm 
M + 2m 



-2Qz{l- z){l-2z){V*k) 



M 



M + 2m 



2 ■ 



'l-2z){ek)<^2 



(5.42) 



and for D wave type mesons are 



= -QMz{l-z)(k^-^pi)^2; 

= {eV*)p^[m^^2 + (k^i)] + (1 - 2^)'(p' + m^ + Mm){kV*){Si) 
-p'^{ek){V*$i) - m(M + m){ke){kV*)^2 ; 



-M(l - 2z) 



(e$i) ( k'^ - —pI ) + m{M + m){ek)^2 



-2Qz{l -z){l- 2z){V*k){p^ + m^ + Mm)$2 



(5.43) 



Equations (5.42), (5.43) together with expression (5.37) constitute the ultimate sets of all 
helicity amplitudes. They give explicit answers for the vector meson production amplitudes 
within leading-log-approximation. 
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Chapter 6 

Analysis for heavy quarkonia 



The general answers (5.42), (5.43) are of course incomprehensible at a quick glance. Therefore, 
a further analysis is needed to grasp the most vivid features of the results and to disentangle 
s-channel helicity conserving and double helicity flip amplitudes. 

Since in the heavy vector mesons quarks can be treated non-relativistically, further simpli- 
fications in analytical formulas (5.42), (5.43) are possible due to the presence of an additional 
small parameter p^/m^. 

In what follows we will first perform the twist expansion and then relate simplified am- 
plitudes to the decay constants (4.37), (4.38). We will then analyze twist hierarchy of the 
amplitudes and compare results for S vs. D wave states. Though we perform this analy- 
sis for heavy mesons, we wish to stress that all qualitative features {S vs. D difference, Q'^ 
dependence etc.) will hold for light quarkonia as well. 



6.1 Twist expansion 

Here we are going to expand the amplitudes (or to be more exact, the quantities $1 and $2 

2 

(5.34), (5.35)) in inverse powers of the hard scale Q and then perform azimuthal angular 
averaging over 0^;. 

Expanding $2 in twists in the main logarithmic region 

/x^A^^^</^2<Q^ (e.i) 

one observes that twist-1 terms cancel, so one has to retain twist-2 and twist-3 terms pro- 
portional K^: 

$2 = ^ - 6.2 

WW 

— * 

The analogous decomposition for $1 reads 

^1^—4 ^ =r^ 6.3 

Q W Q 

Note that the last term does not contain k^. However, one must track it because it will be 
important in double helicity flip amplitudes. 
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6.2 Twist expansion for S wave type mesons 



With the aid of this decomposition one obtains: 
for amphtude L ^ L 

Q 



1 + 



for amphtude T^T 



(e-y*) 



. Q Q 



+ 



2m 



+ 



;i - 2zf 



M 



M + 2m 



- 1 



M + 2m 2 



(1 - 2zf 2m 
4:z{l-z)M + 2m 



Q 



(6.4) 



(6.5) 



^2/^-T*^2'^' AN 6/?2(l-2z) 



g 2g 

This amphtude is naturaUy spht into s-channel hehcity conserving and double hehcity flip 
parts 



2k' 



m' + 2k\z' + {1 - zY) + 



m 



M + 2m 



k''{l-2{l-2zf) 



/^ ^ N N / ^ A N /t7* A N 6/?(l-2z)2\r (1-2^)^ 2m 

/|^t(A7 = -AW = 4z(l-^)(eA)(l^*A)-=j 1 + 1 + ^ , . 



2g 



4z(l - ^) M + 2m 
(6.6) 



Finally, single spin flip amplitudes are 



2k' 



-2Mz{l - z){l - 2zy^{eA) 

Q 



^ (1 - 2zy 2m 
^ Az{l- z) M + 2m 



2k^-. r.2A;2 M 



IUt = -2Qz{l-z){l-2zY^{V*A)^, 



(6.7) 



' M + 2m ■ 

6.3 Twist expansion for D— type vector mesons 

Here we will need to track higher-twist terms. It will turn out later that leading contributions 
vanish, so twist-3 terms will be crucial for our results. 
For amphtude L — > L one has 



4m n\ 2k 



-QMz{l- z) {k' - —Pi) ■ ^ \1- ^ 



4k' 



Q 



Q 



(6.8) 



For T ^ T amplitude, one obtains 



(e-V*)p2 



2 








(-#) 


+ ^^2k' 


Q 




Q 



+ 



2k' 



-Az{l - z)p' + (1 - 2^)'m(M + m) { -^k\eV*) - -=^{eA){y* A) 



Q 



2g 



6^2(1 -2z) 



1 + 



2 k 



^^2 



4k' 



-m{M + m)-k\eV*)^ U " =2" 



Q 



Q 



(6.9) 
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Note that we kept track of all terms oc Again, one can separate out s-channel helicity 
conserving and double helicity flip parts: 



7|LT(A, = Ay) = (e-y*)=j 



2p2 (m^ + 2P - ) - m(M + m)P ( 1 - ^ 



Q 



Q 



-2 4m ^2 



(6.10) 



Finally, single helicity flipping amplitudes are 



-M(l - 2z) 



-2{1 - 2z){eK)K' 



Am 



M' 



k - —pi + m(M + m)^^(l - 2z)(eA) 



(l_2^)2M(eA) 



r2 4m 2 

k — -p. — m(M + ml — TT 



= -8g^(l-^)(l-2^)2(p2 + m2 + mM)^P(l^*A). 

Q 



(6.11) 



6.4 Final results for S wave mesons 

In order to grasp the major features of various S and D wave amplitudes, further simplifications 
can be achieved if one neglects spherically non-symmetric arguments of ag and gluon density. 
First we rewrite general expression (5.37) in the more convenient form 

A{x, g^ A) = -i /^^ff- / rf-V^^(p^) / • /(7* - V) . (6.12) 

In this expression everything except for integrands 1(7* — > \/) is spherically symmetric, thus 
making it possible to perform angular averaging over flp in these integrands. 



6.4.1 S wave: averaging 

Here all the calculations are rather straightforward. In the non-relativistic case one can ev- 
erywhere put 2; — > 1/2 ; M = 2m = my. The resultant integrands are: 
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/^(T ^T;A^ = Ay) = 



SQM ^2 



(g2 + MY 

8M2 



K, 



I%T ^ T; = -A 



V 



+ M2)2 

16 (eA)(f*A) 
"3"(Q2 + M2)2 



1 + 



96 



-»2 2 



5 M2(g2 + M2)_ 



7^(r ^ L) 

7^(L ^ T) 



64 M(eA) p2 ^2 



3 (g2 + M2)2 M2 
512 Q(y*A) 



15 (g2 + M2)2 M2(g2 + M2) 



(6.13) 



Note several things: since the accurate IS wave differs from the naive 7^ spinorial structure 
only by relativistic corrections, one would obtain the same results in the case of naive qqV 
vertex. The only difference would be only extra factor 2 for the L ^ T amplitude (which is 
higher-twist amplitude, anyway). 

Thus the only thing left is |p| integration. Note that these amplitudes are naturally 
expressed in terms of decay constants. Indeed, in the extremely non-relativistic case expression 
(4.37) turns into 



J 6mv 



(6.14) 



6.4.2 S wave: answers for L — ^ L up to differential cross section 

Here we would like to digress and for the sake of logical completeness show how one obtains 
the final result for the differential cross section with the example of L — > L amplitude. If 
needed, the same can be done for the other amplitudes, so we will do it just once. 
One has: 



/ 



-SQmv 



(g^ + mf.)2 



Q di^dG{x,K^) 2^ / 1r X2^ 



(Q2^^^^2 y G{x,Ql)as{Ql)exp{-^B3jpA'^) j d^p 

327rg 1 X2\r^t 7t2\ /7t2\ ^VT^ 

.^^^^-^exp(--53ipA )G(x,go)a.(go)^/, 



V 



With this result, (6.12) becomes 



is- 



is- 



47r2 
1 67r2 g 
3 my 



64:71^ Q G-a, 
3 mv (g^ + ni 



2 \2 

v) 



cvfv • V^nae 



-G-a,exp(-i53ipA2) 



(g2 + m^)2 
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The expression for the differential cross section reads {t = A^): 



da 
Itt 



1 



9 



G'ajexpj-Bsjpt) 

+ ^2^)4 



Finally, one can express this cross section through r(V^ — > e+e ) (see (4.39)): 



da 
'dt 



gW-r(v"^e+e-) 



G^Q;^exp(-i?3ipt) 



(g2 + m^)4 



(6.15) 



6.4.3 S wave: the other amphtudes 

T — > T, — Ay 

. In the non-relativistic case, this amplitude is readily obtained from the above formulas after 
Q my replacement in the numerator of the amplitude (see (6.13)). This means in particular 
that in this limit 

Ris)^f^\' (6.16) 
\AttJ ml 



T ^T,X^ = -X 



V 



. This amplitude is very interesting because of the competition of two very different terms — 
soft and hard scale contribution. Indeed, integration over gluon loop gives 

G(x,/x^) 96 G{x,q')p' 

We see that the soft contribution turns out to be of leading twist, while the pQCD contribution 
is of higher twist. This observation was first made in [6] for the naive type of qqV vertex; here 
we see that it also holds for accurate S and D wave vector mesons. 



T and L^T 

. In the case of heavy quarkonia these single spin flip amplitudes have suppressing non- 
relativistic factors. Besides, the amplitude L ^ T is of twist 3, which is another source of 
suppression. Their ratios to A{T ^ T) = A{T ^ T; = A^) read 



A{T ^L) 8 (eA) 



A{L T) 64 g(y*A) 



A{T ^T) 3 my A{T ^ T) 15 g^ + ml 

The model dependent quantities W2 and are defined via 

1 Jd^^pp^V'S 1 /rf=^Pp"V'5 



W2 = 



W4 = 



my I d'^p ips 
Within the oscillator ansatz (4.42) their values are 

3 1 



ml Jd^pi/Js 



W2= 7: 



2 (myi?)2 ' 



W4 ~ 



15 1 

T {mvRY ■ 



U'4 . 



(6.18) 



(6.19) 



(6.20) 
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6.5 Final results for D wave 



This case is much more tricky. It turns out that the leading terms in integrands pro- 
portional to m^|pP cancel out after angular averaging, so that many new terms, including 
higher twist terms come into play. This cancellation is in fact quite understandable. Indeed, 
in the very beginning we showed that vertex u'jf^u contains both 5* and D waves, with D wave 
probability being suppressed for heavy quarks due to non-rclativistic motion. This means in 
particular that the photon couples to qq pairs sitting either in 5" or D wave state. However, 
at the other end of the quark loop we have a vector meson in pure D wave state. Therefore, 
the largest items in (75+d|---|Vd) cancel out due to S-D orthogonality. 

6.5.1 D wave: Q,p averaging for L ^ L amplitude 

If we limited ourselves only to the leading p^/m^ terms, we would get 

/dn,(P-frf)^4..(?p^-2.1p^)^0, 

which is the manifestation of S-D orthogonality. Thus, we see that p^/m^ terms vanish after 
angular averaging. Therefore, one has to be extremely careful now and must take into account 
all possible sources of p^/m^ terms. To do so, one has to perform the following averaging 

1 /rs 4m 2 



Q 




Az{l-z).-,.[k^-^pi).\l-^\) (6.21) 



Before performing the averaging, let us generate a list of useful formulas : 

m = l,p' m) = ^,p' {pIpI) = Y5P' (6-22) 

Finally, remember that pi = |(1 — 22;)^M^. 
One has to perform the following avegaring 

(4.(l-.).i,.(?-^p;).(l-f 1) (6.23, 

Note that all factors should be carefully examined; all four do contribute to the final answer. 
Decomposing Q as 

Q'^m' + z{l - z)Q' ^m' + ^Q'-^{l- 2zfQ' ^qI- , (6.24) 

with Qq = m"^ + (5^/4, one gets 
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So, omitting Qq ^, one has 



1 - 



M2 



1 - 



4P 



(6.26) 



With the aid of (6.12), one obtains 



2 2 4m 2 



3^ 3M' 



2 



P - 



4 .(2 6^^^^p, 



rP 



4p2 



Vl5 15/ M2go 
16p^ 8g2p4 16p4 



15 ~ 15 



4 
4 

5^ 



* 15 157 



3"M + 2m 15M2 i^M^Q^ I5g^ 
4p4 



(6.27) 



3M2 \ 

4p4 
15M2 



4 8 g2 16 M2 \ 

^ ^ 5 " 5g2 + M2 ~ Yg2 + M2 J 
I g2 + M2 i 



to give 



/ d^p p^ V'D ^ /^^) 
6471^ g G • ^3 • exp(-iS3ipA2) 



One can now again express the integral over quark loop through the decay constant (see 
(4.38)): 

(6.28) 
(6.29) 



15 fJly yyq/ -p iiuy 

Comparison with L ^ L amplitude reveals that 



(g2 + mp2 J 1 




(6.30) 



6.5.2 D wave: the other amphtudes 

For the helicity conserving amplitude one has to repeat the same averaging procedure. The 
calculation proceeds as follows: 



2p2 ( + 2P - 4P^ ) - m(M + m)P ( 1 - ^ ) - 2p (P - 



4P 



M' 



- 2m V - 3m(M + m)p^ + 2^^ (-M^- - 3M^-) + -p^ - -p^ + p^- 



4 , " 4 , 16 M2 8 g2 4 

' " ' rP + 



3^ 15^ ^3^ ^15M2 + g2^ ' 15M2 + g2 



2pM 1 + 



4 M2 
15 M2 + g2 



(6.31) 
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The result can be written as 



m 



,2 

v_ 



m 



TT ]^ _|_ 4 f^v 



In the case of double helicity flip amplitude we again have contributions from soft and hard 
scales with the same hierarchy of twists. Namely, 



G{x,^x') 96 G(x,g')p^ 

7M2(Q2 + M2) ^^-^^^ 

so we again see the soft domination in the double helicity flip amplitude. 

In the case of single spin flip amplitudes, no dangerous calcellations among leading terms 
arise. Before giving a list of amplitudes, we wish to emphasize that in the case of D wave 
mesons there is no non-relativistic suppressing factors like ^2 and defined in (6.19). This 
means that for moderate momentum transfers helicity non-conserving amplitudes are abso- 
lutely important for that case of D wave mesons. 



6.6 S wave vs. D wave comparison 

We would like to present our final results in the form which stresses the remarkable differences 
between S wave and D wave amplitudes. Below we give a table of the ratios 

for helicity conserving and single spin flipping amplitudes. Double spin flip amplitudes are 
not given due to the presence of incalculable non-perturbative contributions. 



Pll = 
Ptt — 




m 



V 



15 g2 + 

2 

mf 



Ptl = -- — 1 + 3 



PLT = J— (6.35) 
8 

Thus, we note several things. First, the abnormally large higher twist contributions to D 
wave amplitudes are seen here as terms oc m'y/iQ'^ + rriy). They even force the opposite sign 
of L — > L amplitude in the moderate domain. Second, we see highly non-trivial and even 
non-monotonous dependence of {All/AttY ratio, which will lead to the presence of a 
dip in experimentally measured oj^joT for D wave meson production. Finally, we must stress 
that in the case of D wave mesons there is no non-relativistic suppression for single spin flip 
amplitudes as it was in S wave mesons. This leads us to a conclusion that s-channel helicity 
is strongly violated in the case of D wave meson production. 



63 



Part III 
Numerical analysis 
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Chapter 7 

Determination of the unintegrated 
gluon structure function of the proton: 
DGD2000 analysis 



The familiar objects from Gribov-Lipatov-Dokshitzer-Altarelli-Parisi (DGLAP) evolution de- 
scription of deep inelastic scattering (DIS) are quark, antiquark and gluon distribution func- 
tions qi{x,Q'^),qi{x,Q^), g{x,Q^) (hereafter x,Q^ are the standard DIS variables). At small 
X they describe the integral flux of partons with the lightcone momentum x in units of the 
target momentum and transverse momentum squared < and form the basis of the highly 
sophisticated description of hard scattering processes in terms of coUinear partons [18]. On 
the other hand, at very small x the object of the Balitskii-Fadin-Kuraev-Lipatov evolution 
equation is the differential gluon structure function (DGSF) of the target [19, 30] 

91ogQ2 ^'-'^ 

with G{x,Q'^) = xg{x,Q'^) (evidently the related unintegrated distributions can be defined 
also for charged partons). For instance, it is precisely DGSF of the target proton that emerges 
in the familiar color dipole picture of inclusive DIS at small x [31, 32, 33] and diffractive DIS 
into dijets [34]. Another familiar example is the QCD calculation of helicity amplitudes of 
diffractive DIS into continuum [35, 36] and production of vector mesons [37, 38]. DGSF's 
are custom-tailored for QCD treatment of hard processes, when one needs to keep track 
of the transverse momentum of gluons neglected in the standard coUinear approximation 
[39, 40, 41, 42]. 

In the past two decades DGLAP phenomenology of DIS has become a big industry and 
several groups — notably GRV [43], CTEQ [44] & MRS [45] and some other [46] — keep 
continuously incorporating new experimental data and providing the high energy community 
with updates of the parton distribution functions supplemented with the interpolation routines 
facihtating practical applications. On the other hand, there are several pertinent issues — the 
onset of the purely perturbative QCD treatment of DIS and the impact of soft mechanisms 
of photoabsorption on the proton structure function in the region of large being top ones 
on the list — which can not be answered within the DGLAP approach itself because DGLAP 
evolution is obviously hampered at moderate to small Q^. The related issue is to what extent 
the soft mechanisms of photoabsorption can bias the dependence of the proton structure 
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function and, consequently, determination of the gluon density from scaling violations. We 
recall here the recent dispute [47] on the applicability of the DGLAP analysis at <^ 2-4 
GeV^ triggered by the so-called Caldwell's plot [48] . Arguably the K-factorization formahsm of 
DGSF, in which the interesting observables are expanded in interactions of gluons of transverse 
momentum k changing from soft to hard is better suited to look into the issue of soft-hard in- 
terface. Last but not least, neglecting the transverse momentum k of gluons is a questionable 
approximation in evaluation of production cross sections of jets or hadrons with large trans- 
verse momentum. It is distressing, then, that convenient, ready-to-use, parameterizations of 
DGSF are not yet available in the hterature. 

Here wc perform a simple phenomenological determination of the DGSF at small x based 
on the 1978 Baltskii-Lipatov (BL) scheme [30], in which the DGSF is directly related to the 
physical observable — the proton structure function F2p{x, Q^). In early 90's the BL scheme 
has been extended to other observables and dubbed K-factorization [49]; it is also closely 
related to the color dipole factorization in the color dipole BFKL approach [31, 32, 33]. Our 
interest is in producing the ready-to-use Ansatz for JF(x,k^), so that wc take advantage of 
large body of the early work on color dipole BFKL factorization [32, 50, 51] and follow a very 
pragmatic strategy first applied in [35, 36]: (i) for hard gluons with large k, we make as much 
use as possible of the existing DGLAP parameterizations of G{x, k^), (ii) for the extrapolation 
of hard gluon densities to small we use an Ansatz [34] which correctly describes the color 
gauge invariance constraints on radiation of soft gluons by colour singlet targets, (iii) as 
suggested by color dipole phenomenology, we supplement the density of gluons with small 
by nonperturbative soft component, (iv) as suggested by the soft-hard diffusion inherent to 
the BFKL evolution, we allow for propagation of the predominantly hard-interaction driven 
small-x rise of DGSF into the soft region invoking plausible soft-to-hard interpolations. The 
last two components of DGSF are parameterized following the modern wisdom on the infrared 
freezing of the QCD coupling and short propagation radius of perturbative gluons. Having 
specified the infrared rcgularization, we can apply the resulting jF(a;, /?^) to evaluation of the 
photoabsorption cross section in the whole range of small to large Q^. 



7.1 The Ansatz for differential gluon structure function 

The major insight into parameterization of DGSF comes from early experience with color 
dipole phenomenology of small-x DIS. In color dipole approach, which is closely related to 
K-factorization, the principal quantity is the total cross section of interaction of the qq color 
dipole r with the proton target [32, 52, 53] 

a{x, r) = ^j ^ ^J^, ^^ as (max{^^^ -} j T{x, k') , (7.2) 
which for very small color dipoles can be approximated by 

a(x,r)^^as(^)G(x,^) , (7.3) 



where A «i 10 comes from properties of the Besscl function Jo (2;). The phenomenological 
properties of the dipole cross section are well understood, for extraction of g{x^ r) from the 
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experimental data see [20, 54]. The known dipole size dependence of a{x,r) serves as a 
constrain on the possible K^-dependence of J-'{x, k^). 

As we argued in section 3.2, DGLAP fits are likely to overestimate J^hard{x, k^) at moderate 
if. Still, approximation (7.3) does a good job when the hardness A/r^ is very large, and at 
large we can arguably approximate the DGSF by the direct differentiation of available fits 
(GRV, CTEQ, MRS, ...) to the integrated gluon structure function Gpt{x,Q'^): 



Hereafter the subscript pt serves as a reminder that these gluon distributions were obtained 
from the pQCD evolution analyses of the proton structure function and cross sections of related 
hard processes. 

The available DGLAP fits are only applicable at if > Q^, sec table 1 for the values of 
Ql. in the extrapolation to soft region if < Q\ wc arc bound to educated guess. To this end 
recall that perturbative gluons are confined and do not propagate to large distances; recent 
fits [55] to the lattice QGD data suggest Yukawa-Debye screening of perturbative color fields 
with propagation/screening radius i?c ~ 0.27/m. Incidentally, precisely this value of Rc for 
Yukawa screened colour fields has been used since 1994 in the very successful color dipole 
phenomenology of small-x DIS [50, 51]. Furthermore, important finding of [51] is a good 
quantitative description of the rising component of the proton structure function starting 
with the Yukawa-screened perturbative two-gluon exchange as a boundary condition for the 
color dipole BFKL evolution. 

The above suggests that dependence of perturbative hard J-'hard{x, if) in the soft region 
if < Q\ is similar to the Yukawa-screened fiux of photons in the positron, cf. eq. (2.33), with 
ttem replaced by the running strong coupling of quarks Cpasi^^) and with factor N^. instead 
of 2 leptons in the positronium, for the early discussion see [34], 

T^\k-) = CpN^^^ f ^#^1 V^(^) , (7.5) 

Here \Xpt — — ^■'^^ GeV is the inverse Yukawa screening radius and must not be interpreted 
as a gluon mass; more sophisticated forms of screening can well be considered. Following 
[50, 51, 56, 22] we impose also the infrared freezing of strong coupling: ^^(k^) < 0.82; recently 
the concept of freezing coupling has become very popular, for the review see [57]. 

The vertex function Vf^{K) describes the decoupling of soft gluons, k <C from color 
neutral proton and has the same structure as in eq. (2.34). In the nonrelativistic oscillator 
model for the nucleon one can relate the two-quark form factor of the nucleon to the single- 
quark form factor, 

F2(^, -K) = F, {^^/) ■ (7.6) 

To the extent that i?^ <C Pf^ the detailed functional form of -F2(k, —if) is not crucial, the simple 
relation (7.6) will be used also for a more realistic dipole approximation 



^,(^^)^^. (7.7) 
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The gluon probed radius of the proton and the charge radius of the proton can be somewhat 
different and A ~ 1 GeV must be regarded as a free parameter. Anticipating the forthcoming 
discussion of the diffraction slope in vector meson production we put A = 1 GeV. 

As discussed above, the hard-to-soft diffusion makes the DGSF rising at small-x even in 
the soft region. We model this hard-to-soft diffusion by matching the dependence (7.5) to 
the DGLAP fit J-'pt{x,Ql) at the soft-hard interface and assigning to Thard{^'-i^^) in the 
region of < Ql the dependence of the Born term (7.5) and the x-dependence as shown 
by the DGLAP fit J^pt{x, Ql), i.e., 

j'Har.ix, i^) = 4f\'^" fsil^;! (^iQi - + ^vt{^. ^'m' - Ql) ■ (7.8) 

^pt '{Qi) 

Because the accepted propagation radius Rc ~ 0.3 fm for perturbative gluons is short 
compared to a typical range of strong interaction, the dipole cross section (7.2) evaluated 
with the DGSF (7.8) would miss an interaction strength in the soft region, for large color 
dipoles. In Ref. [50, 51] this missing strength for large dipoles has been modeled by the non- 
perturbative, soft mechanism with energy-independent dipole cross section, whose specific 
form [50, 38] has been driven by early analysis [56] of the nonperturbative two-gluon exchange 
and tested against the diffractive vector meson production data [38]. More recently several 
closely related models for (Tsoft('") have appeared in the literature, see for instance models for 
dipole-dipole scattering via polarization of non-perturbative QCD vacuum [58] and the model 
of soft-hard two-component pomeron [59]. In the spirit of eq. (7.4) one can parameterize 
interaction of large color dipoles in terms of the genuinely soft, nonperturbative component of 
DGSF. The principal point about this non-perturbative component of DGSF is that it must 
not be subjected to pQCD evolution. Thus the arguments about the hard-to-soft diffusion 
driven rise of perturbative DGSF even at small do not apply to the non-perturbative DGSF 
and we take it energy- independent, 

•^Sk^, = asoftC^Nj^^ ( ^f, ] V^in) , (7.9) 

TT \K -i-fJ'soft/ 

where /i^^j^ <^ /x^^. Furthermore, it is natural that the soft component of DGSF decreases in 
the perturbative domain of ^ /ip^ faster than the perturbative Born term (7.5), which is 
achieved by the extrapolation of the form suggested in [35, 36] 



T{X, K^) = J^softi^^ 2 + ^hard{x, /^^) ^ 2 , (7.10) 



with Kg ~ fJ'pt- 

The above described Ansatz for DGSF must be regarded as a poor man's approximation. 
The separation of small-/?^ DGSF into the genuine nonperturbative component and small-/?^ 
tail of the hard perturbative DGSF is not unique. Specifically, we attributed to the latter 
the same small-x rise as in the DGLAP fits at Q^, though one can not exclude that the hard 
DGSF has a small x-independent component. The issues of soft-hard separation and whether 
the non-perturbative component of DGSF enters different observables in a universal manner 
must be addressed in dynamical models for infrared regularization of perturbative QCD and 
non-perturbative QCD vacuum and only can be answered confronting such models to the 



68 



experiment. We recall that in the conventional DGLAP analysis the effect of soft gluons is 
reabsorbed into the input gluon distributions. 

The K-factorization formulas (2.40) and (2.41) correspond to the full-phase space extension 
of the LO DGLAP approach at small x. For this reason our Ansatze for !Fhard{x-, Q"^) will be 
based on LO DGLAP fits to the gluon structure function of the proton Gpt{x, Q'^). Wc consider 
the GRV98LO [43], CTEQ4L, version 4.6 [44] and MRS LO 1998 [45] parameterizations. We 
take the liberty of referring to our Anzatze for DGSF based on those LO DGLAP input as 
D-GRV, D-CTEQ and D-MRS parameterizations, respectively. 

Our formulas (2.40), (2.41) describe the sea component of the proton structure function. 
Arguably these LL^ formulas are applicable at x < Xq = 1-^3- 10~^. At large the 
experimentally attainable values of x are not so small. In order to give a crude idea on finite- 
energy effects at moderately small x, we stretch our fits to x ^ xq multiplying the above 
Ansatz for DGSF by the purely phenomenological factor {1 — x)^ motivated by the familiar 
large-x behaviour of DGLAP parameterizations of the gluon structure function of the proton. 
We also add to the sea components (2.40), (2.41) the contribution from DIS on valence quarks 
borrowing the parameterizations from the respective GRV, CTEQ and MRS fits. The latter 
are only available for > Q^. At a; < 10^^ this valence contribution is small and fades away 
rapidly with decreasing x, for instance see [51]. 



7.2 The parameters of DGSF for different DGLAP in- 



Our goal is a determination of small-x DGSF in the whole range of by adjusting the relevant 
parameters to the experimental data on small-x F2p{x, Q^) in the whole available region of 

as well as the real photoabsorption cross section. The theoretical calculation of these 
observables is based on Eqs. (2.40), (2.41), (7.10). 

The parameters which we did not try adjusting but borrowed from early work in the color 
dipole picture are Rc — 0.27 fm, i.e., jipt — 0.75 GeV and the frozen value of the LO QCD 
coupling with Aqcd = 0.2 GeV: 



We recaU that the GRV, MRS and CTEQ fits to GSF start the DGLAP evolution at quite a 
different soft-to-hard interface Ql and diverge quite a lot, especially at moderate and small 
i^. The value of Ql is borrowed from these fits and is not a free parameter. 

The adjustable parameters are fj.soft, cisoft, i^u,d, ^1 and k\ (for the heavier quark masses 
we take = rriu^ + 0.15GeV and = 1.5 GeV). The both rriu^d and ^soft have clear 
physical meaning and we have certain insight into their variation range form the early work 
on color dipole phenomenology of DIS. The role of these parameters is as follows. The quark 
mass defines the transverse size of the qq — uu, dd Fock state of the real photon, whose 
natural scale is the size of the p-mcson. Evidently, roughly equal values of F2p{x, Q^) can be 
obtained for somewhat smaller JF(,t, Q^) at the expense of taking smaller rriu^d, i-e. larger size 
of the photon, and vise versa. Therefore, though the quark mass does not explicitly enter the 



puts 




(7.11) 
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parameterization for jF(a;, k^), the preferred value of rriuj could have been correlated with the 
DGLAP input. We find that it is sufficient to take the universal rriu^d = 0.22 GeV. 

The Hggjf defines the soft scale in which the non-perturbative glue is confined, and controls 
the r-dependence of, and in conjunction with asoft sets the scale for, the dipole cross section for 
large size qq dipoles in the photon. We find that it is sufficient to take the Tinivcrsal a,^oft = 2 
and = 3 GeV^ for the parameter of suppression of the hard tale of non-perturturbative soft 
glue. 

The magnitude of the dipole cross section at large and moderately small dipole size de- 
pends also on the soft-to-hard interpolation of DGSF, which is sensitive to DGLAP inputs 

for pcrturbative component Gpt{x,Q'^). This difference of DGLAP inputs can be corrected 
for by adjusting fJ^l^fi and the hard-to-soft interface parameter The slight rise of k.^ helps 
to suppress somewhat too strong x-dependence of the soft tale of the perturbative glue. The 
specific parameterizations for depend on the DGLAP input and are presented in table 1. 
Only kI and Hsoft varied from one DGLAP input to another. The soft components of the 
D-GRV and D-CTEQ parameterizations turn out identical. The eye-ball fits are sufficient for 
the purposes of the present exploratory study. The parameters found are similar to those used 
in [35, 36] where the focus has been on the description of diffractive DIS. 



Table 1. The parameters of differential gluon structure function for different DGLAP inputs. 





D-GRV 


D-MRS 


D-CTEQ 


LO DGLAP input 
Ql GeY' 
kI GeV^ 
Hsoft: GeV 


GRV98LO [43] 
0.895 

(l + 0.0018 Vi)'^' 
0.1 


MRS-LO-1998 [45] 
1.37 

(l + 0.038 log^i)'^' 
0.07 


CTEQ4L(v.4.6) [44] 
3.26 

(l + 0.047 V i)'^' 
0.1 



One minor problem encountered in numerical differentiation of all three parameterizations 
for Gpt{x, Q"^) was the seesaw K^-behavior of the resulting DGSF (7.4), which was an artifact 
of the grid interpolation routines. Although this seesaw behavior of DGSF would be smoothed 
out in integral observables like G{x,Q'^) or F2p{x,Q'^), we still preferred to remove the un- 
physical seesaw cusps and have smooth DGSF. This was achieved by calculating DGSF from 
(7.4) at the center of each interval of the Q^-grid and further interpolating the results between 
these points. By integration of the so-smoothed J-'pt{x, Q"^) one recovers the input Gpt{x, Q^). 
The values of Ql cited in Table 1 corresponds to centers of the first bin of the corresponding 
Q'-grid. 

7.3 The description of the proton structure function 

We focus on the sea dominated leading log^ region of a; < 10~^. The practical calculation 
of the proton structure function involves the two running arguments of DGSF: Xg and k^. 
We recall that in the standard collincar DGLAP approximation one has <^ <^ and 
Xg ~ 2x, see eq. (3.1). Within the ^-factorization one finds that the dominant contribution to 
F2p{x, Q'^) comes from M^ ~ Q'^ with httle contribution from M'^ ^ Q^. Because at small Xg 
the Xg dependence of J-{xg, Q^) is rather steep, we take into account the Xg — x relationship 
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Figure 7.1: The K-factorization description of the experimental data on F2p{x,Q'^) in the low 
region; black circles are ZEUS BPC data [62], open triangles denote HI shifted vertex (SV) 
data [64], open squares are E665 data [65]. Solid line represents K-factorization results for the 
D-GRV parameterization of the differential gluon structure function T{x, k^). 

(2.44). Anticipating the results on effective intercepts to be reported in section 7, we notice 
that for all practical purposes one can neglect the impact of k on the relationship (2.44), which 
simplifies greatly the numerical analysis. Indeed, the Xg dependence of J^{xg, n^) is important 
only at large which contribute to F2p(x, Q"^) only at large Q^; but the larger Q^, the better 
holds the DGLAP ordering < A;^, Q^. Although at small to moderate the DGLAP the 
ordering breaks down, the Xg dependence of J-{xg, k^) is weak here. 

Obviously, achieving a good agreement from small to moderate to large is a highly non- 
trivial task, because strong modification of the soft contribution to J-{x, k^) unavoidably echos 
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□ E665 



NMC 



Figure 7.2: T/ie K-factorization description of the experimental data on F2p{x, Q^) in the 
moderate and high region; black circles and triangles are ZEUS data [60], [61], open circles 
and triangles show HI data [63], [64], open squares are E665 data [65], stars refer to NMC 
results [66]. Solid line represents K-factorization results for the D-GRV parameterization of 
the differential gluon structure function J-{x, k^) . 



in the integrated gluon SF throughout the whole range of and shall affect the calculated 
structure function from small to moderate to large Q^. 

The quality of achieved description of the experimental data on the small-x proton structure 
function is illustrated by figs. 7.1, 7.2. The data shown include recent HERA data (ZEUS 
[60], ZEUS shifted vertex [61], ZEUS BPC [62], HI [63], HI shifted vertex [64]), FNAL E665 
experiment [65] and CERN NMC experiment [66]. When plotting the E665 and NMC data, 
we took the liberty of shifting the data points from the reported values of to the closest 
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Figure 7.3: A comparison of the K-factorization description of the experimental data on 
F2p{x,Q'^) for several values of with the D-GRV, D-CTEQ and D-MRS parameterizations 
of the differential gluon structure function T(x, k^). The contribution to F2p{x, Q"^) from DIS 
off valence quarks is shown separately for larger Q^. 



boxes for which the HERA data were available. For Q'^ < Q'^ — 0.9 GeV^ the parameterizations 
for valence distributions are not available and our curves show only the sea component of 
F2p{x,Q'^), at larger the valence component is included. 

At X < 10^^ the accuracy of our D-GRV description of the proton structure function is 
commensurate to that of the accuracy of standard LO GRV fits. In order not to cram the 
figures with nearly overlapping curves, we show the results for D-GRV parameterization. The 
situation with D-GTEQ and D-MRS is very similar, which is seen in Fig. 7.3, where wc show on 
a larger scale simultaneously the results based on the D-GRV, D-CTEQ and D-MRS DGSFs 
for several selected values of Q^. Here at large we show separately the contribution from 
valence quarks. The difference between the results for F2p{x, Q"^) for different DGLAP inputs 
is marginal for all practical purposes, see also a comparison of the results for cr'"'^ for different 
DGLAP inputs in Fig. 7.4. 
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Figure 7.4: A comparison of the K-factorization description of the experimental data on real 
photoah sorption cross section based on the D-GRV, D-CTEQ and D-MRS parameterizations 

of the differential gluon structure function Tix^ k^). The squares show the experimental data 
from 1992-93 direct measurements, the bullets are the results of extrapolation of virtual pho- 
toabsorption to = ([62] and references therein). The soft component of photoabsorption 
cross section is shown separately. 

7.4 Real photoabsorption cross section a^^ 

In the limiting case of = the relevant observable is the real photoabsorption cross section 
a"'^. Although the Bjorken variable is meaningless at very small Q^, the gluon variable Xg 
remains well defined at = 0, see eq. (2.44). In Fig. 7.4 we present our results alongside the 
results of the direct measurements of cr'''^ and the results of extrapolation of virtual photoab- 
sorption cross sections to — 0, for the summary of the experimental data see [62]. The soft 
contribution to the cross section is shown separately. We recall that our parameterizations 
for T{x,K^) give identical soft cross sections for the GRV and CTEQ inputs (sec table 1). 
The barely visible decrease of o^1lf^ towards small 1^ is a manifestation of oc (1 — x)^ large-a; 
behaviour of gluon densities. The extension to lower energies requires introduction of the 
secondary reggeon exchanges which goes beyond the subject of this study. 

We emphasize that we reproduce well the observed magnitude and pattern of the energy 
dependence of a'^^ in an approach with manifestly energy-independent soft contribution to 
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the total cross section. In this scenario the energy dependence of (7"<^ is entirely due to the 
Xp-dependent hard component J-'hardi^g, Q"^) and as such this rise of the total cross section for 
soft reaction can be regarded as driven entirely by very substantial hard-to-soft diffusion. Such 
a scenario has repeatedly been discussed earher [50, 51, 67]. Time and time again we shall 
see similar effects of hard-to-soft diffusion and vise versa. Notice that hard-to-soft diffusion 
is a straightforward consequence of full phase space calculation of partonic cross sections and 
we do not see any possibility for decoupling of hard gluon contribution from the total cross 
sections of any soft interaction, whose generic example is the real photoabsorption. 
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Chapter 8 

Properties of differential gluon 
structure function 

8.1 DGSF in the momentum space 
8.1.1 Soft/hard decomposition of DGSF 



Now we focus on the x and behavior of the so-determined DGSF starting for the reference 
with the D-GRV parameterization. The same pattern holds for DGSF based on CTEQ and 
MRS DGLAP inputs, see below. In figs. 8.1 and 8.2 we plot the differential gluon density 
T{xg, Q"^), while in Fig. 8.3 we show the integrated gluon density 



Here the subscript D is a reminder that the integrated Goix^Q"^) is derived from DGSF. As 
such, it must not be confused with the DGLAP parameterizations Gpt{x, Q"^) supplied with 
the subscript pt. 

Figs. 8.1 and 8.2 illustrate the interplay at various x of the nonperturbative soft component 
of DGSF and perturbative hard contribution supplemented with the above described contin- 
uation into < Ql. The soft and hard contributions are shown by dashed and dotted lines 
respectively; their sum is given by solid line. 

Apart from the large-x suppression factor {1 — x)^ our non-perturbative soft component 
does not depend on x. At not so small x — 10~^ it dominates the soft region of ^ 1 2 
GeV^, the hard component takes over at higher k^. The soft-hard crossover point is close to 
/Xpj, but because of the hard-to-soft diffusion it moves with decreasing x to a gradually smaller 

-»2 

In this determination of DGSF we focus on the ready-to-use parameterizations, the dynam- 
ical evolution properties of the so-found DGSF will be addressed elsewhere. In what concerns 
the relationship between DGSF and the observable proton structure function, the early work 
by Kwiecinski et al. [68] is close in spirit to ours, the difference being in a treatment of the 
nonperturbative soft component and subjecting DGSF to unified BFKL/DGLAP evolution. 
In Fig. 8.2 we present DGSF read off the plots in [68]; the agreement with our results is good, 
which indicates a consistency of our purely phenomenological parameterizations of DGSF with 
general expectations from the BFKL dynamics. 




(8.1) 
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Figure 8.1: D-GRV differential gluon structure function T{x^ k^) as a function of at sev- 
eral values of X. Dashed and dotted lines represent the soft and hard components; the total 
unintegrated gluon density is shown by the solid line 

8.1.2 Soft/hard decomposition of the integrated gluon structure 
function 

The role of the soft component if further illustrated by Fig. 8.3, where we show the inte- 
grated gluon density (8.1) and its soft and hard components Gsoft{x,Q'^) and Ghard{x,Q'^), 
respectively. The soft contribution Gsoftix,Q'^) is a dominant feature of the integrated gluon 
density Gd{x,Q^) for < 1 GcV^. It builds up rapidly with and receives the major 
contribution from the region ~ 0.3 -r- 0.5 GeV^. Our Ansatz for J^softi^, k^) is such that it 
starts decreasing already at ~ 0.2 GeV^ and vanishes rapidly beyond > i^-ioft^ ^gs- 
10,11. Still the residual rise of the soft gluon density beyond ~ 0.5 GeV^ is substantial: 
Gsoftix,Q^) rises by about the factor two before it flattens at large Q^. We emphasize that 
GsoftiQ'^) being finite at large is quite natural — a decrease of Gsojt{Q^) at large only is 
possible if J-'soft{Q^) becomes negative valued at large Q^, which does not seem to be a viable 
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Figure 8.2: The same as in Fig. 8.1 hut overlaid onto one graph for illustration of the x- 
dependence of J-'{x, k,"^) . The dashed lines shows the soft component J- so fti^, k,"^) and its slight 
variation with x due to the finite-x factor (l — x)^. The dot-dashed curves show the Kwiecinski 
et al. [68] results for DGSF from a ^-factorization phenomenology of F2p{x, Q^) based on the 
solution of the unified BFKL/DGLAP equation. 

option. 

At moderately small x — 10~^ the scaling violations are still weak and the soft contribution 

Gsoft{x,Q'^) remains a substantial part, about one half, of integrated GSF Gd{x,Q'^) at all 
Q^. At very small x < 10~^ the scaling violations in the gluon structure function are strong 
and Ghardix.Q"^) > Gsoftix.Q"^) starting from Q'^ ~ 1-2 GeV^. 

8.1.3 Soft/hard decomposition of the proton structure function 

Eqs. (2.40), (2.41) define the soft/hard decomposition of the proton structure function. In 
Fig. 8.4 we show F2p^'^{x, Q^) and F2p^^{x, Q^) as a function of for two representative values 
of X. Notice how significance of soft component as a function of Q'^ rises from fully differential 
J^{x,Q'^) to integrated Gr,{x,Q'^) to doubly integrated F2p^*{x,Q'^). At a moderately small 
X ~ 10~^, the soft contribution is a dominant part of F2p{x,Q'^), although the rapidly rising 
hard component F2p'^'^{x, Q^) gradually takes over at smaller x. 

Notice that not only does F2p^^{x, Q^) not vanish at large Q^, but also it rises slowly with 
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Figure 8.3: The same as Fig. 8.1 hut for integrated gluon structure function Gd{x,Q^) as 
given by the D-GRV parameterization of the differential gluon structure function T{x, k^), for 
the discussion see Section 6.2. 



Q"^ as 



^2p {x,Q)^ }_^ef — log-——. 



(8.2) 



Again, the decrease of F2p'^^{x, Q'^) with would only be possible at the expense of unphysical 
negative valued GgoftiQ"^) at large Q"^. 



8.2 DGSF in the x-space: effective intercepts and hard- 
to-soft diffusion 

It is instructive to look at the change of the x-dependence from the differential gluon structure 
function T{x,n^) to integrated gluon structure function Gd^x^Q"^) and further to proton 
structure function F2p{x,Q'^). It is customary to parameterize the x dependence of various 
structure functions by the effective intercept. For instance, the effective intercept Tg// for 
differential gluon structure function is defined by the parameterization 



■l\^e//('t^) 

J-{x, K^) OC ( — 



(8.3) 



One can define the related intercepts Thard for the hard component J^hard{x, k^)- Notice, that 
in our Ansatz Tgoft = 0. 
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Figure 8.4: The soft-hard decomposition of K-factorization results for the proton structure func- 
tion F2p{x,Q'^) evaluated with the D-GRV parameterization of the differential gluon structure 
function J^{x, k^) . 



The power law (8.3) is only a crude approximation to the actual x dependence of DGSF 
and the effective intercept Tgff will evidently depend on the range of fitted x. To be more 
definitive, for the purposes of the present discussion we define the effective intercept as 



^,,l0g[J^{x2,K')/T{x,,K')] 



^'"^^ ' = — — 

taking X2 = 10~^ and Xi = 10~^. The effective intercept Thardi^'^) is defined by (8.4) in terms 
of ^hard{x,i^). 

One can define the related intercepts Ag//, Xhard for the integrated gluon structure function 
GD{x,Q^y. 

Gd{x,Q^) oc^-j . (8.5) 

In the case of F2p{x, Q^) we define the intercept A{Q^) in terms of the variable x defined 

as 

^-^F^-^.' (8-6) 

where My is the mass of the ground state vector meson in the considered flavor channel. Such 
a replacement allows one to treat on equal footing <l GeV^, where the formally defined 
Bjorken variable x can no longer be interpreted as a lightcone momentum carried by charged 
partons. For the purposes of the direct comparison with t((5^), X{Q'^) and in order to avoid 
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Figure 8.5: Effective intercepts for total, and hard components of, (a) the differential gluon 

structure function J-'{x,Q'^); (b) integrated gluon structure function Gd{x,Q'^) and (c) pro- 
ton structure function F2p{x,Q'^) evaluated with the D-GRV parameterization of the differ- 
ential gluon structure function T{x,k^). In the box (d) we compare the effective intercepts 
Teff{Q'^),Xeff{Q'^) and Aeff{Q'^) for J^{x,Q'^), Gd{x,Q'^) and F2p{x,Q'^), respectively. 

biases caused by the valence structure function, here we focus on intercepts ^effi^hard for 
the sea component of the proton structure function F2p°'{x, Q'^): 

F,7(x,g2)ocy . (8.7) 

The results for the effective intercepts arc shown in figs. 8.5, 8.6 and 8.7. 

In our simplified hard-to-soft extrapolation of J-'hard{x, Q^) we attribute to J-'hard{x, Q^) 
at < Ql the same x-dependence as at = Ql modulo to shght modifications for the 
x-dependence of /?|. This gives the cusp in ThardiQ'^) at = Ql, i.e., the first derivative of 
ThardiQ"^) is discontinuous at = Ql- 

A comparison of Fig. 8.2 with Fig. 8.3 and further with Fig. 8.4 shows clearly that only 
in DGSF J-'{x, Q^) the effect of the soft component is concentrated at small Q^. In integrated 
Gd{x,Q^) and especially in the proton structure function F2p{x,Q^) the impact of the soft 
component extends to much larger Q'^. The larger the soft contribution, the stronger is the 
reduction of Te// from Thard and so forth, the pattern which is evident from Fig. 8.5a to 8.5b 
to 8.5c, see also figs. 8.6 and 8.7. 

The change of effective intercepts from differential J-'{x,Q'^) to integrated Gd{x,Q^) is 
straightforward, the principal effect is that XhardiQ'^) < ThardiQ'^) and XeffiQ"^) < Teff{Q'^) 
which reflects the growing importance of soft component in Gd{x, Q"^). The change of effective 
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Figure 8.6: Effective intercepts for total, and hard components of, (a)the differential gluon 

structure function J-'{x,Q'^); (b) integrated gluon structure function GniXyQ"^) and (c) pro- 
ton structure function F2p{x, Q"^) evaluated with the D-MRS parameterization of the differ- 
ential gluon structure function J-'{x,k,'^). In the box (d) we compare the effective intercepts 
Teff{Q'^),Xeff{Q'^) and Aeff{Q'^) for J^{x,Q'^), Gd{x,Q'^) and F2p{x,Q'^), respectively. 

intercepts from T{x, Q^) and Gd{x, Q"^) to F2p{x, Q"^) is less trivial and exhibits two dramatic 
consequences of the hard-to-soft and soft-to-hard diffusion. If the standard DGLAP contri- 
bution (3.2) were all, then the change from the intercept A((5^) for integrated gluon density 
to the intercept A((5^) for the proton structure function F2p{x,Q'^) would be similar to the 
change from t{Q^) to X{Q^), i.e., the effective intercept Aeff{Q^) would have been close to 
zero for <1 GeV^. However, by virtue of the hard-to-soft diffusion phenomenon inherent 
to the K-factorization, F2p{x, Q^) receives a contribution from gluons with > Q^, which 
enhances substantially A^ardiQ^) and Ae//(Q^). The net result is that at small to moderately 
large we find AhardiQ"^) > XhardiQ"^) and Aeff{Q'^) > XeffiQ"^)- As we emphasized above 
in sections 5.3, the rise of real photoabsorption cross section is precisely of the same origin. 

The second effect is a dramatic flattening of effective hard intercept, AhardiQ"^)-, over the 
whole range of Q^. For all three DGLAP inputs A^ardiQ'^) flattens at approximately the same 

^hard ~ 0.4. 

The whole set of figs. 8.5-8.7 also shows that the systematics of intercepts in the hard 
region of > Ql is nearly identical for all the three DGLAP inputs. In the soft region 
we have a slight inequality Thard{l<-'^)\D-MRS ^ '^hard{i^'^)\D-GRV^ which can be readily at- 
tributed to a slight inequality Ql{MRS) > Ql{GRV). In the case of CTEQ4L(v.4.6) in- 
put the value of Ql{CTEQ) is substantially larger than Ql{M RS) , Ql{G RV) . In the range 
Ql{MRS),Ql{GRV) < R'^ < Ql{CTEQ) the effective intercept rhard{K^) rises steeply with 
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Figure 8.7: Effective intercepts for total, and hard components of, (a) the differential gluon 

structure function J-'{x,Q'^); (h) integrated gluon structure function Gd^x^Q"^) and (c) pro- 
ton structure function F2p{x,Q'^) evaluated with the D-CTEQ parameterization of the differ- 
ential gluon structure function J-'{x,k,'^). In the box (d) we compare the effective intercepts 
Teff{Q'^),Xeff{Q'^) and Aeff{Q'^) for J^{x,Q'^), Gd{x,Q'^) and F2p{x,Q'^), respectively. 

K^. This explains a why in the soft region Thard{f^'^)\cTEQ significantly larger than for the 
D-GRV and D-MRS parameterizations. The difference among intercepts for the three pa- 
ramctcrizations decreases gradually from differential J-{x,k^) to integrated Gd{x,Q^) gluon 
density to the proton structure function F2p{x,Q^). 

Finally, in Fig. 8.8 we compare our results for Ae//((5^) with the recent experimental data 
from ZEUS collaboration [61]. Since in the experimental fit the range of x — [xmax,Xmin] 
varies from point to point, we mimicked the experimental procedure in our evaluation of 
Ag// from eq. (8.11) by taking X2 = Xmax and xi = Xmin- This explains the somewhat 
irregular dependence. The experimental data include both sea and valence components. 
At > Q^{GRV) = 0.9 GeV^ we included the valence component of the structure function 
taking the GRV98LO parameterization. For CTEQ4L(v.4.6) and MRS-LO-1998 the values of 

are substantially larger. However, the valence component is a small correction and we took 
a liberty of extracting the valence contribution F2p\x, Q"^) from GRV fits for Ql{GRV) < < 
Ql{MRS),Ql{CTEQ). The overall agreement with experiment is good. Difference among 
the three parameterization is marginal and can of course be traced back to figs. 8.5-8.7. 
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Figure 8.8: Effective intercepts A((5^) of the proton structure function F2p{x, Q^) in the HERA 
domain evaluated for the D-GRV, D-MRS and D-CTEQ parameterizations for the differential 
gluon structure function J-{x, k^); the experimental data points are from ZEUS [61] 

8.3 How the gluon densities of factorization differ from 
DGLAP gluon densities 

It is instructive also to compare our results for integrated GSF (8.1) with the conventional 
DGLAP fit Gpt{x,Q'^). In Fig. 8.9 we present such a comparison between our integrated D- 
GRV distribution (the solid curves) and the GRV98LO distribution (the dashed curves). As 
was anticipated in section 3.2, at very large the two gluon distributions converge. We also 
anticipated that at small x and moderate the DGLAP gluon structure functions Gpt{x, Q"^) 
are substantially larger than the result of integration of DGSF, see eq. (8.1). At a; = 10~^ 
they differ by as much as the factor two-three over a broad range of Q"^ ^100 GeV^. The 
difference between integrated DGSF and the DGLAP fit decreases gradually at large x, and 
is only marginal at x = 10~^. 

Recall the substantial divergence of the GRV, MRS and CTEQ gluons structure functions 
of DGLAP approximation Gpt{x, Q"^) at small and moderate Q^. Contrary to that, the k- 
factorization D-GRV, D-CTEQ and D-MRS gluon structure functions Gd{x,Q'^) are nearly 
identical. We demonstrate this property in Fig. 8.10 where we show integrated Goix, Q"^) and 
their DGLAP counterparts Gpt{x, Q^) for the three parameterizations at two typical values of 
X. Because of an essentially unified treatment of the region of < Ql and strong constraint 
on DGSF in this region from the experimental data at small Q^, such a convergence of D-GRV, 
D-CTEQ and D-MRS DGSFs is not unexpected. 

One can also compare the effective intercepts for our integrated GSF Gd{x,Q'^) with 
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Figure 8.9: Comparison of our results for integrated gluon density Gd{x,Q^) evaluated with 
the D-GRV parameterization of the differential gluon structure function J^{x,k^) with the 
GRV98L0 DGLAP input parameterization Gpt{x, Q'^). 
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Figure 8.10: A comparison of the divergence of GRV98L0, CTEQ4L(v.4.6) and MRS-LO-1998 
gluon structure functions Gpt{x, Q"^) in the left box with the divergence of our integrated gluon 
structure functions Gd{x, Q^) evaluated for the D-GRV, D-CTEQ and D-MRS parameteriza- 
tions for differential gluon structure function T{x^ Q^) at two typical values of x 
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Figure 8.11: A comparison of the intercept A^*j((5^) of the x-dependence of the GRV98L0, 
CTEQ4L(v.4-6) and MRS-LO-1998 gluon structure functions Gpt{x.Q'^) with their counter- 
part Aejj(Q^) for integrated Gd{x,Q'^) evaluated with D-GRV, D-GTEQ and D-MRS param- 
eterizations for differential gluon structure function J-'{x, Q"^). 

those obtained from DGLAP gluon distributions Gptix^Q"^). Fig. 8. 11 shows large scattering 
of \^!}]{Q^) from one DGLAP input to another. At the same time, this divergence of differ- 
ent DGLAP input parameterizations is washed out to a large extent in the /^-factorization 
description of physical observables (see also (8.8)). 



8.4 How different observables probe the DGSF 

The issue we address in this section is how different observables map the dependence 
of T{xg^ K^). We expand on the qualitative discussion in section 3.2 and corroborate it with 
numerical analysis following the discussion in [52]. We start with the two closely related quan- 
tities — longitudinal structure function Fl{x,Q'^) and scaling violations dF2{x,Q'^)/dlogQ'^ 
— and proceed to F2p{x, Q^) and the charm structure function of the proton F2p{x, Q^). This 

mapping is best studied if in (2.40) and (2.41) we integrate first over k and z. In order to 
focus on the dependence we prefer presenting different observables in terms of T{2x, k^) 
and GdC^x, k^) 

Fl(x, Q') = E 4 / ^')^(2^, , (8-8) 



dF2{x,Q^) ^ as{Q^ )^^2 f dn 

aiogg2 ~ 



j:e}J^ei"~\Q^K'M2x,K'). (8.9) 
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Figure 8.12: The weight function Ql for mapping of the differential gluon structure function 
K^) as a function of for several values of Q^. We show separately the results for light 
flavours, u, d, and charm. 

In the numerical calculation of Fl{x,Q'^) starting from eq. (2.41) we have Xg and as 
the two running arguments of J^{xg, k^). As discussed above, the mean value of Xg is close to 
2x, but the exact relationship depends on k^. The A;, z integration amounts to averaging of 
J^{xg, N?) over certain range of Xg. The result of this averaging is for the most part controlled 
by the effective intercept Teffin^): 

{J^{xg, K^)) = h^{2x, ti^) i-^j j = r{ie)T{2x, te) . (8.10) 

Because the derivative of Teffin^) changes rapidly around — Ql, the rescaling factor r{K^) 
also has a rapid variation of the derivative at — Ql, which in the due turn generates 

the rapid change of derivatives of ^LpiQ'^ ^ i^'^) around — Ql- As far as the mapping of 
differential J-{2x, k^) is concerned, this is an entirely marginal effect. However, if we look at 
the mapping of integrated gluon structure function Gd{x,Q'^), which is derived from (8.8), 
(8.9) by integration by parts: 

(8.11, 

aiogg^ ^ — 3^ -W diog^ ^ ^ ' ^^-^^^ 



87 




Figure 8.13: The weight function 02 for mapping of the differential gluons structure function 
K^) as a function of for several values of Q"^ . We show separately the results for light 
flavours, u, d, and charm. 

then the weight functions dQ2^[\Q'^,K^)/d\ogK^ will exhibit a shghtly irregular behaviour 

around = Q"^. Evidently, such an irregularity appears in any region of fast variation of 
Tf.ff^K^)] in our simplified model it is somewhat amplified by the cusp-like dependence of 

Finally, starting from (8.12) one obtains a useful representation for how the proton struc- 
ture function F2p{x, Q^) maps the integrated gluon structure function: 

= ^E4/ ^wi"\Q',K')as{K')Gn{2x,K') (8.13) 

In figs. 8.12 and 8.13 we show the weight functions Ql and G2. Evidently, for light flavours 
and very large they can be approximated by step-functions 

e^i\Q',k')^9{CL,2Q'-K'), (8.14) 

where the scale factors C^, ~ | and C2 ~ 2 can be readily read from figures, for the related 
discussion sec [52]. Note that the value C2 ~ 2 corresponds to (72 ~ 8 introduced in Section 
3.2. Recall that the development of the plateau-like behaviour of 0^ and ©2 which extends 
to ~ signals the onset of the leading logQ^ approximation. For large in the approxi- 
mation (8.14) the K,^ integration can be carried out exphcitly and Fl{x, Q"^) oc Gd{2x, ClQ"^)- 
Similarly, dF2{x,Q'^)/d\ogQ'^ oc G'd(2x, CaQ^), cf. eq. (3.6). 
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Figure 8.14: The same as Fig. 8.12 but for mapping of the integrated gluon structure function 
Gd{x, k^) as a function of for several values ofQ^. We show separately the results for light 
flavours and charm. 
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Figure 8.15: The weight function W2 for mapping of the integrated gluons structure function 
Gd{x, k^) as a function of for several values of Q^. We show separately the results for light 
flavours and charm 

Still better idea on how and scaling violations map the integrated GSF is given by 
figs. 8.14,8.15, where we show results for —dQ{ff)j^/d\ogK^ and W2'^'^\ The first quantity 
is sharply peaked at ~ ClQ^- The second quantity visibly develops a plateau at large 
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Figure 8.16: A comparison of the experimental data from ZEUS [69] on the charm struc- 
ture function of the proton with K-factorization results for F2'^{x,Q^) based on the D-GRV 
parameterization of the differential gluon structure function J^{x, Q"^). 

Q"^ . As can be easily seen, scaling violations do receive a substantial contribution from the 
beyond-DGLAP region of > Q"^. 

Because of the heavy mass, the case of the charm structure function F|p(a;, Q"^) is somewhat 
special. Figs. 8.14 and 8.15 show weak sensitivity of F2p{x, Q'^) to the soft component of 
T{x,K^), which has an obvious origin: long- wavelength soft gluons with k <j rUc decouple 
from the color neutral cc Fock state of the photon, which has a small transverse size <j 
Our results for F2p{x, Q^) are shown in Fig. 8.16, the agreement with the recent precision 
experimental data from ZEUS [69] is good. 
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Chapter 9 



Improved determination of the 
differential glue in proton: DGD2002 
analysis 

During the last two years new high-accuracy measurements of structure function in the 
expanded parameter space were presented by ZEUS [70], [71] and HI [72]. In the hght of 
our analysis, the most important improvement was further plunging into \ow-x region (down 
to X = 6.7 ■ 10~^, [70]) and 2 ^ 3% accurate determination of in the small-Q^ (0.25 
GeV^ < < 0.65 GeV^) and relatively small-a; (10~^ < x < 10~^) region. Even a brief 
comparison of /cffactorization predictions for in this region based on the differential glue 
from the previous section showed a systematic several-sigma deviation from experimental 
points. The second problem with old differential glue was that we did not quite match the 
experimentally measured exponent of rise towards high 1/x. As our analysis showed (see 
Fig. 8.8), the predicted intercept was 1 -j- 2(T lower at > 10 GeV^. 

Thus, not being able to claim that we reproduce with old differential glue the F2p data 
well enough, we re-extracted the differential gluon structure function. This time, the eye-ball 
fits were not acceptable, and therefore a ^^-minimization procedure was carried out. This 
communication presents the results of this reevaluation. 

9.1 Fitting procedure and parameters of DGSF 

Our goal is a determination of small-x DGSF in the whole range of by adjusting the relevant 
parameters to the experimental data on small-x F2p{x, Q^) in the whole available region of Q^. 
The /ct-factorization predictions for F2p were calculated at iV = 191 points in {x, Q^) space 
and compared to the experimentally values measured at HERA. These points include all data 
available at the moment for low-Q^ region (Q^ < 1 GcV^) as well as all data points for several 
higher values of Q\ namely, at = 1.5 2.5 2.7 3.5 4.5 5.0 6.5 10 35 GeV^.The was 
calculated according to 

N ( zptheor. rpexp.\2 

x- = E ' ' . (9-1) 

j=l ^ stat ' ^syst 

The parameters which we did not try adjusting but borrowed from early work in the 
color dipole picture are Rc — 0.27 fm, i.e., Upt — 0.75 GeV and the frozen value of the 
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LO QCD coupling with Aqqd = 0.2 GeV (7.11). Our previous analysis showed that the 
parametrizations D-GRV, D-MRS, and D-CTEQ (that is, DGSF based on GRV, MRS, and 
CTEQ fits to integrated gluons) run very close to one another throughout the whole x, 
space. This can be verbalized as no matter from what particular DGLAP parametrization we 
start, we will arrive at virtually the same shape of the differential gluon structure function. 
To this end, we used only GRV98LO fits, since they are available for the widest range of x 
and values. 

The adjustable parameters are Hsoft, C'soft, ^u,d, i^l and (for the heavier quark masses 
we take rUg — mu,d + O.lSGeV and rUc — 1.5 GeV). The both m^^d and yLgoft have clear 
physical meaning and we have certain insight into their variation range form the early work 
on color dipole phenomenology of DIS. The role of these parameters is as follows. The quark 
mass mu,d defines the transverse size of the qq = uu, dd Fock state of the real photon, whose 
natural scale is the size of the p-meson. Evidently, roughly equal values of F2p{x, Q^) can be 
obtained for somewhat smaller T{x, Q"^) at the expense of taking smaller mu,d and vise versa. 
Therefore, though the quark mass does not explicitly enter the parameterization for jF(a;, k^), 
the preferred value of rriu^ could have been correlated with the DGLAP input. We find that 
it is sufficient to take the universal mu,d = 0.21 GeV, which is slightly lower than 0.22 GeV, 
used before. 

Parameter nj^ff defines the soft scale in which the non-perturbative glue is confined, and 

controls the r-dependence of, and in conjunction with agoft sets the scale for, the dipole cross 
section for large size qq dipoles in the photon. We find that it is sufficient to take the universal 
Hsoft = 0.1 GeV and = 1 GcV^. 

The magnitude of the dipole cross section at large and moderately small dipole size de- 
pends also on the soft-to-hard interpolation of DGSF, which is sensitive to DGLAP inputs for 
perturbative component Gpt{x,Q'^). This difference of DGLAP inputs can be corrected for 
by adjusting the hard-to-soft interface parameter k,\. The slight rise of k\ helps to suppress 
somewhat too strong ,T-dcpcndcncc of the soft talc of the perturbative glue. The specific pa- 
rameterizations for depend on the DGLAP input and are presented in Table 2. Only k\ 
and Hsoft varied from one DGLAP input to another. 

In order to be able to assess the uncertainty in determination of differential glue, we 
performed several minimization procedures using slightly different sets of free parameters. 
The resulting parameters of the fits are shown in Table 2; below we comment on each fit in 
detail. 

Fit 1 

When obtaining Fit 1, we adjusted values of agoft, transition point in (7.8), functional form 
of 

i^h — '^hi^)j which we took as a first order polynomial a + b\og{l/x), and the power ( of 
the high-x suppressing factor (1 — x)''. 

The resulting value of ( turned our uncomfortably large ( ^ 11. Numerically, such a strong 
suppression serves as a remedy against somewhat too slow x behavior of as we approach 
X ~ 10^^ from the low-x side. Although rather artificial, this suppressing factor does not 
invalidate our approach since the exact behavior of the glue in the limit x — > 1 lies far beyond 
the scope of the present approach. Still, we would like to note the alarming tendency that 
even with this factor large part of the overall comes precisely from the region we wanted to 
correct (10~^ < x < 10~^). This might be an indication that our understanding of this region 
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is far from perfect. More analysis is needed to settle up this issue. 



Fit 2 

In Fit 2 we opened up another degree of freedom, namely, we allowed for shifting in the 
second argument of G{x, Q^). As early analyses showed, the best fcj-factorization vs. DGLAP 
correspondence would be 



,2, dG{x,Q') 



(9.2) 



The early analysis gave C ~ 2, but in our approach we treated C as free parameter. 

We started from Fit 1 and let the slide to its minimum as we freed parameter C. We 
expected C > 1 at minimum point. Indeed, our predictions with Fit 1 for high-Q^ region 
suffered from too high values at high x and too low values at very small x. As early analysis 
indicated, when we shift the scale according to (9.2), we make the unintegrated gluon 
density increase a,t x < 10~^ and decrease at higher values of x. We expect this tendency to 
survive the multiple integration procedure and to echo in the proton structure F2p. 

The minimization procedure gave C 1.1 with slight adjustment of other parameters, see 
Table 1 for details. 



Fit 3 



In the region of very hard gluons both Fit 1 and Fit 2 rely on the leading order DGLAP 
parametrizations of G{x,Q^). Although it is desirable that integrated gluon structure func- 
tion Gd{x, Q^) based on our parametrizations approaches in the double logarithmic limit the 

conventional gluon density obtained from DGLAP evolution, there is certainly no requirement 
that our fits be built on these DGLAP fits. One should only make sure that at log(l/a;) ^ 1 
and log{Q^ / Aqqjj) ^ 1 behavior of our fits is compatible with the corresponding behavior of 
the DGLAP fits. 

The properties of the DGLAP-evolved gluon density in the limit is well understood. Since 
the anomalous dimension of gluons is higher than that of the sea quarks, the secondary gluons 
in this limit tend to be radiated off gluons as well. The evolution of the integrated gluon 
density G{x, Q"^) = xg{x, Q^) separates out and is governed by 



dG{x,Q^) _ as{Q^ 
d log (5^ 



with splitting function 



Pgg{z) = 27Ve 
The Regge-type behavior 



27r 



z{l -z) + 



J^dzG(^^,Q'^PGG{z), 



(1-z) 



+ 



1-z 



+ 



+ ^7Ve5(l - Z) 



G{x,Q')=fiQ')-Q' 



(9.3) 



(9.4) 



(9.5) 



with constant A is compatible with the DGLAP equations and leads to 

z 1 — z 



ds Jx 



z(l-z) + 



(1 - 



+ 
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= s. 



1 

A 



1 



1 



A + 1 A + 2 A + 3 12 



+ + C + V(A + 2) 



(9.6) 



Here 



log 



log (gVA 



2 

QCD 



log f g^/A^cD 



(9.7) 



and C and ^ are the Euler constant and the digamma function respectively. One immediately 
gets the solution for the integrated gluon structure function 



log 



(9.8) 



which leads to the differential glue of the form 



J-'{x, K^) OC 



log 



K 



S-1 



X 



(9.9) 



This analysis inspires us search for a parametrization of Tpt that would be power-like both in 

log (k^/^qcd) log(l/^) ill the double logarithmic regime. 
The functional form of J^pt Fit 3 is based on is 



:^pt(^,«')|«« = 0.245 



log 



A, 



QCD, 



I) 



0.4 



(9.10) 



Table 2. The parameters of DGD2002 fits to differential gluon densities. 





Fit 1 


Fit 2 


Fit 3 


hard input 


GRV98LO 


GRV98LO 


Eq.(9.10) 


/?^-shift 




C = 1.1 




Ql GeV^ 


1.45 


1.45 


1.4 


dsoft 


2.66 


2.63 


2.6 


kI GeV^ 


0.4 + 0.245 log i 


0.38 + 0.245 logi 


0.31 log i 


total 


257 


245 


226 



9.2 The properties of the gluon structure function 

In Fig. 9.1 we plotted the three fits to the differential gluon structure functions obtained above 
vs. at several values of x. One observes here a number of interesting features. First, the 
three curves display rather similar behavior at very small x, but at x as high as 10~^ the 
difference among them throughout the region > 1 GeV^ is quite sizable. Still, all of these 
parametrizations of DGSF do provide a reasonably accurate description of F2p. Thus, we 
conclude that modern experimental data on do not place severe restrictions on the shape 
of differential glue at not very small x (say, x > 10~^). 

The second feature of the curves presented is their salient two-peak shape. The technical 
origin of this clear separation of the soft and hard exchange mechanisms is, of course, the 



94 



F(x,/c'), 2002 




-2-1 2 3-2-1 2 3 

10 10 1 10 10 10 10 10 1 10 10 10 

(GeV') 1^ (GeV') 
Fit 1 Fit 2 Fit 3 

Figure 9.1: Differential gluon structure function as function of k^: DGD2002 analysis 

very abrupt extrapolation of soft gluons into hard region and vica versa, generated by high 
powers C- This might seem artificial, but as we described, such abrupt extrapolation seemed 
necessary in order to obtain the correct behavior of effective intercept of and therefore it 
was quite essential for getting better x^- Thus, we inclined to think that such clear separation 
of the soft and hard mechanisms is indeed preferred by experiment. 

The integrated gluon structure function is shown in Figs. 9.2. A similar observation, 
though in a subdued form, can be made. At x — 10~^ 10~^ one can notice some departure 
among the curves, which die out as we shift towards lower and lower x. As we go to smaller 
X, wc observe plateau around ~ 1 -j- 3 GeV^, which originates from the two-peak shape of 
J-{x, K^), become more and more prominent. 

Note also that at small x and very large all three Gd{x,Q^) curves, including Fit 3 
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Gd(x,Q'), 2002 
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Fit 1 



Fit 2 



Fit 3 



Figure 9.2: Integrated gluon structure function as function of : DGD2002 analysis 



with simple-formula parameterization of J-'hard, approach each other. 



9.3 The observables 



9.3.1 Structure function and its derivatives 

Since the overall behavior of structure function F2p{x, Q"^) was the quantity we tried to fit, 
one can expect a very good description of the data. Indeed, as Figs. 9.3 and 9.4 show, our 
calculations for based on all three DGSF fits go almost directly through the experimental 
points. This trend is somewhat spoiled beyond the region of fitting, namely, at a; > 0.01, 
but still we do not run into any severe discrepancy even here. Note also that throughout the 



96 



2p 



0.3 
0.25 

0.2 
0.15 

0.1 
0.05 















0^= 0.045 GeV^ 















































10 



0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 



10 



10 



10 



low Q 

0.3 
0.25 

0.2 
0.15 

0.1 
0.05 





Q^= 0.11 GeV^ 






































10 



10 





















"""■"■^ 






= Q^= 0,40 GeV" 











0.7 
0.6 
0.5 
0.4 
0.3 
0.2 
0.1 




10 



10 



10 

X 



10 10 10 10 

kj— factorization description 

DGD2002, Fit 1 

DGD2002, Fit 2 

DGD2002, Fit 3 



10 



10 









































= 0^= o.e 


5 GeV^ 
















10 10 10 10 

# ZEUS BPC 1995 

A ZEUS BPT 1997 

A HI SV 

□ E665 



10 

X 



Figure 9.3: The kt-factorization results for the structure function F2p{x, Q'^) in the small 
region 
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Figure 9.4: The kt-factorization results for the structure function F2p{x, Q^) in the moderate 
and high region 
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fitting region all three curves difi^er less than 5%. 




Figure 9.5: The improvement of the kt-factorization predictions for the structure function 
based on DGD2000 (left pane) and DGD2002 (right pane) parametrizations of the differential 
gluon structure function 

Fig. 9.5 illustrates the improvement in the description of the structure function by the 
/cj-factorization calculations as we switch from old DGD2000 parametrizations of DGSF to the 
new DGD2002 fits. The curves now go directly through very constraining new data points 
and therefore have less tendency to depart from each other outside the fitting range. 

Fig. 9.6 shows the effective intercept of structure function in the moderate domain 
together with recent HI data [73]. The intercepts were calculated according to 

l0g(x2/Xi) 

For each differential gluon density fit we plotted here two curves: the upper one corresponds 
to the effective intercept taken between xi — 10~^ and X2 — 10~^, while the lower one is for 
xi — 10""^ and X2 — 10"^. Significant deviation between the two curves indicates the fact that 
the powerlike law 

F2p{x,Q^)cx (-j (9.12) 



X 



is only a very crude approximation. 

The agreement with the data is reasonable, especially when one takes into account that 
at lower one should compare the data with the upper curves and at higher the data 
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Figure 9.6: The kt-factorization predictions for the effective intercepts of the structure function 
in the moderate region confronted with HI data (solid circles) 



should be compared with the lower curves. This is due to the experimental procedure used 
by HI to determine the intercepts: at smaller the value of the intercept comes from data 
points in the range x e (10~^, 10~^), while at higher only x ~ 10~^ 10~^ were available. 

9.3.2 Structure function Fl 

In 1999 HI pubhshed improved data on determination of structure function Fl{x,Q'^). We 
show them in Fig. 9.7 together with our predictions. The general agreement can be seen, 
however, since the data are not very constraining, little futher information can be extracted 
from this plot. 

9.3.3 Real photoabsorption cross section 

Finally, in Fig. 9.8 we show experimental data for real photoabsorption cross section a^p 
together with our predictions. 
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Figure 9.7: The longitudinal structure function Fl of the proton in the moderate region 
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Figure 9.8: The predictions of the real photoabsorption cross section as function of total energy 
based on DGD2002 parametrization 
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Chapter 10 



Numerical analysis of the vector 
meson producton 

10.1 IS states: a brief look 

10.1.1 Absolute values of cross sections and scaling phenomenon 

An overview of the experimental results on vector meson production together with our predic- 
tions is given by Fig. (10.1). Here we plotted the experimentally measured values of p meson 
cross sections for photoproduction (HI [74] and ZEUS [75]) and electroproduction (HI [76] 
and ZEUS [77]) cases, mesons electroproduction cross sections (HI [78] and ZEUS [79], [80]), 
J/t/j photoproduction (HI [81]) and electroproduction cross sections (HI [82] and ZEUS [77], 
[83]), and T(IjS') photoproduction cross sections (HI [81] and ZEUS [84]). All experimental 
data points are either taken at 1^ = 75 GeV or are extrapolated to this energy. Whenever 
possible, we used the energy dependence measured experimentally in the corresponding pa- 
pers. In the case of T mesons photoproduction no reliable data on energy are available, and 
we used power law cr{T) oc with 5 — 1.75, which comes from our calculations. 

In the case of 0, J/ijj, and T mesons the cross sections were nmltiplied by appropriate 
flavor factors in order to remove trivial sensitivity of the cross sections to the mean-square 
charge (ef ) of the quark content in a vector meson. 

As suggested by the twist expansion analysis, the total production cross sections should 
exhibit an approximate scaling in variable Q . Indeed, as Fig. (10.1) shows, the experimental 
data for various vector mesons do possess such a scaling property: data points for p, (j), and 
J ftp taken at the same values of Q2 almost coincide. It is worth noting that the scaling 

— 2 

phenomenon takes place even at small Q . Note that the /c^-factorization predictions also 
exhibit approximate scaling phenomenon. 

10.1.2 The energy and |t|-dependence 

2 

The remarkable scaling in variable Q is observed not only in the magnitude of the production 
cross sections, but also in the patterns of energy dependence and |t|-dcpcndcnce. 

The energy dependence of the vector meson production cross sections is compatible within 
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cr(7P ^Vp), nb 




Figure 10.1: The total cross section of diffraction electroproduction of various vector mesons 
against scaling variable Q^+rriy. Results are scaled according to flavor content to equilibrate the 
electric charges of different flavours. The kf factorization predictions based on oscillator (solid 
lines) and suppressed Coulomb (dashed lines) are also shown. All calculation are performed 
for W — Ih Ge V; the experimental points are either taken in this energy range or consistently 
extrapolated to this energy. 

experimental errors with the power-hke Regge-type growth: 

a(xW^; 5 = 4[q;]p-1]. (10.1) 

The latter equality reflects the assumption that the energy behavior comes from the gluon 
content of the proton, which is usually linked to the Pomcron intercept ajp at \t\ = 0. 

The t-dependence of the differential cross section da/d\t\ can be approximated at \t\ ^0.5 
GeV^ by a simple exponential law 
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The magnitude of the slope parameter B shows how "fragile" the proton and the produced 
meson are. 
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Figure 10.2: The effective intercepts (upper pane) and effective slopes (lower pane) of the vector 
meson production cross sections vs. scaling variable + rriy. The data points represent p 
meson and J/%1) meson results. The curves are the ktO- factorization calculations for the p 
meson with the oscillator wave function. 

Both quantities — the intercept and the slope — do exhibit dependence, but again via 
scaling variable Q . Fig. 10.2 shows a plot from [16], where the fct-factorization predictions 
for these quantities are compared with the experimental results then available. Although the 
agreement is not perfect, the tendency is caught by the /cj-factorization calculations. Below 
we will give more detailed investigation of both quantities. 

10.1.3 The level of accuracy 

Before proceeding further, we should have a clear idea of what level of accuracy one can expect 
from the /cj-factorization predictions. 

There is a number of quantities that can be a source of uncertainty in the final answers. 
Below we discuss them. 
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The gluon density 

The gluon content of the proton is an intrinsically non-perturbative quantity, and therefore 
it is not calculable within the /crfactorization approach. Our extensive analysis of the structure 
function lead us to an accurate determination of the differential gluon density T{x, k^) in 
the proton. The /ct-factorization results for various physical observables based on distinct fits 
to T{x, K^) differ from each other at the level of several percent. 

As we discussed above, the vector meson production amplitudes are related not to the 
diagonal quantity jF(a;, k^), but to the off-diagonal (skewed) gluon distribution JF(a;i, X2-, A). 
The situation can be partially cured by linking the skewed gluon distribution to the diagonal 
one by means of Shuvaev formula (5.28) or its simplified version (5.30). Since such a linking 
is strictly valid only when xi or X2 vanishes, and therefore it is a good approximation in the 
cases Q"^ 3> my and ^ ruy. 

If ~ rriy, this linking becomes not well justified, and the whole procedure introduces an 
uncertainty. In order to test the magnitude of the uncertainty associated with the skewness 
of the gluon density, we performed an additional check, namely we calculated the p meson 
cross sections at low using strictly forward and simphfied off-forward (5.30) Ansatze for 
the gluon structure function. We found that at = using forward instead of non-forward 
(the later is default option for all calculations here) reduced the cross section by a factor of 
1.07. Obviously, smallness of this factor originates from low average value of energy growth 
exponent A. 

Thus, the inaccuracy introduced at low by forward/off- forward Ansatz for gluon density 
is no more than 10%. 

The wave function 

With the gluon density being brought under reasonable control, the only major uncalculable 
piece of the pie is vector meson wave function. As described above, we focus on the soft 
wave function, in particular, we used the oscillator and "suppressed Coulomb" wave function 
Ansatze. Being virtually the two limiting cases of how the radial part of the wave function can 
look like, they represent fairly well the region of uncertainty introduced by a specific choice of 
the wave function. 

As can be seen from the above Figures, the calculations based on the oscillator wave 
function Ansatz are roughtly twice higher than those obtained with the suppressed Coulomb 
WF. 

The width of vector meson 

In our calculations we treated the produced vector mesons as particles with negligible width. 
This is not the case for the p meson, whose width is about 1/5 of its mass. Usually, the 
incorporation of a finite width of a produced particle is conducted via effective 'smearing' of 
the results (which depend on the mass of the particle produced) over a certain mass interval. 
If the cross sections calculated for a given mass happen to have convex dependence on mass 
(which is precisely the case in the p meson production), then such a smearing will lower 
the values of the cross sections. Thus, from a very general arguments we can expect such a 
smearing in our case as well. Roughly, it should amount to decreasing of the cross sections by 
factor of {1 + Tp/rup) fa 1.2. 
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A more accurate calculation of this effect is a non-trivial task. The problem is that in our 
treatment of the vector meson production we never refer to the vector meson mass. We deal 
only with the effective invariant mass of the qq pair. Therefore accurate calculation of the 
smearing effect requires solution of conceptually non-trivial problems. 

The limits of /cj-factorization approach 

Finally, the very approach we use has limited applicability domain. In particular, it would 

be fallacious to extend our calculations to high region. A rough criterion to border of the 
applicability domain can be given by Qmai ~ y/^'^l-^lofu which is about 50 GeV^ for the HERA 
energy range. Above this values, the logarithms \og{Q'^ / ^1^^^ will be more important than 
log(l/a;), and one can expect that /cj-factorization will underestimate observed cross sections. 

10.2 The p meson production 

In this section we provide a detailed description of the /ct-factorization predictions for the p 

meson production. Note that throughout this chapter we treat the physical p meson as pure 
IS state. Whether this is indeed realized in Nature and what changes if S/D wave mixing 
occurs, will be discussed in Chapter 11. 

Whenever experimental data available, we compare them with our results. 

10.2.1 Q"^ dependence 

The Q^-dependence of the p meson production cross section is shown in Fig. 10.3. One sees 
that as we slide to higher values of Q^, the cross section drops sharply. Although for the major 
part of the region shown the experimental data points fall roughly between the oscillator 
and Coulomb wave function predictions, two separate places of discrepancy are easily visible. 

Low region 

The first problematic point is \ow-Q^ region. As we come from moderate Q"^ down to = 1 
GeV^ and below, our predictions, if compared to the experimental points, tend to climb too 
high. 

The first thought would be to suspect that our predictions rise too steeply as — > 0. An 
accurate analysis shows, however, that this suspicion is off the point. In fact, the ^ 
behavior of our predictions is perfectly compatible with experimentally observed tendencies. 
This ather unexpected fact is illustrated by Figs. (10.4) and (10.5). 

In Fig. (10.4) we zoom in on the region < 1 GeV^, where the available experimental 
data include ZEUS 95 BPC [77], ZEUS 94 photoproducton [75] and HI photoproduction 
[74] measurements. The /ct-factorization prediction based on the oscillator wave function 
(which seems to be a more reasonable choice for the p mesons than the Coulomb WF) and on 
DGD2002, Fit 3 are shown with solid line. The kt predictions overshoot the data throught the 
whole region shown here. However, when simply divided by 2, the predictions pass exactly 
through all the data points (dotted curve in Fig.(10.4)), including the photoproduction point. 
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Figure 10.3: Total cross section of the diffractive p meson production as a function ofQ^. The 
kt-factorization predictions based on oscillator (solid lines) and suppressed Coulomb (dashed 
lines) are also shown. All calculation are performed for W — 75 GeV using DGD2002, Fit 1. 



Alternatively, one can try parametrizing the cross sections at low by a simple formula: 



^(Q^T.^^^; m')^'^, (10.3) 



which is equivalent to 



In order to find parameters Tn^ff, n that reproduce the low-Q behavior of our predictions, we 
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Figure 10.4: Total cross section of the diffractive p meson production as a function of in 
the small-Q^ region. The solid curve represents the kt-factorization predictions; dashed line 
shows the same prediction divided by two. 



plotted in Fig. (10.5) the quantity 

- iog(i + gv</) ^ ^ ^ 

against { = 1 + Q'^/mlff. If Eq.(10.4) holds, this quantity should be equal to n for all values 
of ^. Wc sec that when rnlj^j: ~ 0.4 GeV^, TZ has the most fiat shape vs. ^ and is ~ 2. 

This result is in perfect agreement with low-Q^ analysis of the ZEUS BPC data [77]: when 
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Figure 10.5: Effective + Tn^ff exponent of the transverse p meson production cross section 
in the low-Q^ region as a function of ^ = 1 + Q"^ /mljj at various values of nieff 



fitted by the formula 



the data yield 

rueffiexp) = 0.66 ± 0.11 GeV , 

in perfect agreement with our rrie// ~ 0.6 -r- 0.65 GeV. Note that although in our analysis 
power n was a free parameter, it turned our close to 2. 
Three conclusions can be drawn from this analysis. 

First, we showed that although there is a sizable (by factor of 2) departure between the 
magnitude of the experimentally measured cross sections and the kf factorization predictions. 
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the shape of low-Q^ behavior is perfectly reproduced by the calculations based on the oscillator 
wave function. 

Second, this analysis proves that the value mg// ~ 0.65 GeV is dynamically generated. 
Indeed, the light quark mass used in our approach was pre-fixed by the gluon structure function 
analysis at the level = 0.21 GeV. Still irilff > 'im^, which means that approximately one 
half of fnlj^j comes from quark momentum. The exact proportion depends of course on the 
particular choice of the vector meson wave function. It is still to be investigated whether the 
oscillator wave function was an unexpectedly lucky guess, or whether the perfect reproduction 
of low-Q^ behavior is accidental. 

Third, the above means that the actual region where our predictions and the experimental 
data really mismatch is not the low-Q^, but rather moderate region (Q^ ~ 1 GeV^). 



High Q region 

The second region where our predictions tend to depart from the data is high-Q^ region, 
^ 5 -j- 10 GeV^. If the cross section is (locally) fitted by the powerlike fall-off 

then this discrepancy can be expressed numerically in terms of effective exponent n. The 
experimental determination of this exponent resulted in the following values. The ZEUS 95 
data [77] with > 5 GeV^ are consistent with fit (10.7) with energy independent rif^^p, whose 
avegare value is found to be 2.32 ± 0.10. More copious HI data sample [76] taken at W = 75 
GeV results in Uexp = 2.24 ± 0.09, which is in agreement with ZEUS fit. 

The kt factorization predictions for this exponent is shown in Fig. (10.6). Here we plotted 
local analogue of the exponent riexp, i-e. quantity 

n (q^ ^ ^) ^ I . (10.8) 

V " J ^og[{Qi + ml)/{Qi + ml)] 

Again, if fit (10.7) holds, n should be independent of Q^. 

Graphs in Fig. 10.6 show that this is not the case. At intermediate n starts already 
from about 2.5, then grows s increases, and at > 20 GeV^ it is even higher than 3.0. 

Although being in constrast with fits to experimental data, such a Q'^ growth of n{Q'^) is 
still firmly grounded theoretically. Qualitatively, this can be understood from the analysis of 
the leading logQ^ result. At fixed W'^ and high enough Q^, the dependence of the cross 
section comes from 



a(Q') oc 



l + R{Q 



2^ 



2 



(10.9) 



(Q2 + m2)4 

where c ^ OAl comes from the approximate representation of the off- forward gluon distribu- 
tion (5.30). The dependence of slope B and of running coupling is inessential for our 
point. 

The non-trivial behavior arises from the integrated gluon structure function. In Fig. 10.7, 
left pane, we show its behavior at = 75 GeV. The origin of the peaked shape is obvious: 
at moderate the integrated glue grows due to the sharp explicit dependence (that is. 
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Figure 10.6: Effective + w^e// exponent of the total p meson production cross section in 
the high-Q^ region as a function ofQ^. The solid and dashed curves represent calculations for 
oscillator and suppressed Coulomb wave functions respectively 

due to large values of unintegrated gluon density F{x,k^)), while for larger the effect of 
decreasing Xg// overpowers and leads to decreasing of G as grows further. 

On the right pane of Fig. 10.7 we show the local exponent uq of the gluon density (the 
gluon density contribution to the local exponent n is equal to ^nc) 



One sees that at moderate Q^, when G is stil rising, it tampers the Q'^ fall, but when Q"^ > 10 
GeV^, gluon density starts decreasing on its own. This is precisely the reason why at higher 
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Figure 10.7: (Right pane) The Q'^ behavior of the integrated gluon density G{0A1Q^ /W"^ , (5^/4) 
at fixed value of W = 75 GeV obtained by integration of the DGD2002 Fit 1; (left pane) The 
effective exponent of this integrated gluon density 

the exponent n (Fig. 10.6) boosts up. 

The arguments that justify such a behavior seem to be universal, and to this end, it is 
surprising why the experimental data do not possess such a behavior. 

10.2.2 gl — ot decomposition 

An important insight of the behavior of the p meson production cross sections comes from 
the separate analysis of (TLiQ"^) and (7r((5^), that is proportion of the p meson production 
rates cause by transverse and longitudinal photons. 

Fig. 10.8 represents the results for these cross sections within the /ct-factorization approach 
confronted with experimental data [85]. 

One sees that at high we do provide a reasonably good description of the ctl, but our ctt 
curves sink significantly deeper as grows. Thus, it is mostly a way too steep Q^-behavior 
of (Tt that causes departure of our curves from the data. 

If analyzed in terms of power-like fits 

(JT(g') oc (Q' + m^)--- ; GLiQ'') oc (Q^ + m^""- , (10.11) 
the experimental data yield [85] 

nT{exp) = 2.47 ± 0.03 . (10.12) 

We found no direct expermental results for ul, but clearly it should be even less than ut- 

The /cj-factorization predictions for the local values of and ut, defined similarly to 
(10.8), are shown in Fig. 10.9. At higher Q'^, ul and ut grow up to 3 and 4 respectively, the 
latter being in stark contast to the data. 
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Figure 10.8: Experimental data on transverse ax and longitudinal a l cross sections of p meson 
production compared with the kt-factorization approach. 



Fig. 10.10 depicts the ratio 

(10.13) 

c7t(7 P PP) 

experimental points taken from [76]. The glaring diagrccmcnt at higher is, of course, caused 
by too much suppressed (Tt out calculations predict. Evidently, if we found a way to increase 
(Tt, the ratio R{Q^) would be authomatically cured. 



10.2.3 Energy dependence 

The growth of the vector meson production cross sections is a well-established fact. It is 
linked basically to the energy growth of the Pomeron exchange, and therefore fitting the cross 
sections to the energy power law 

a{W) oc (10.14) 

seems a natural way to quantitize the energy growth. 

Fig. 10.11 shows the experimental data on S from ZEUS [77] and HI [76] together with the 
fcf-factorization predictions based on oscillator and Coulomb wave function. Since the true W 
dependence of the cross sections can deviate from simple power law (10.14), the exponent 6 
can depend on W as well. At lower the upper pair of curves corresponds to 6 calculated 
between 50 and 75 GeV, while the upper pair corresponds to energy range between 110 and 
150 GeV. 
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Figure 10.9: The effective exponents shown separately for ai and cjt as functions of Q"^ . 
Solid and dashed lines correspond to oscillator and suppressed Coulomb wave functions respec- 
tively 

One sees that the agreement is rather good, although a tendency that our curves go shghtly 
higher than the (ZEUS) data is noticeable. However due to still significant experimental errors, 
it is too early to draw any more definite conclusions. 

10.2.4 i-dependence 

The analysis of t-dependence of the differential cross sections within a perturbative frame- 
work has an ambiguous status. On the one side, if we deal with proton in final state, this 
dependence is governed largely by the intrinsically nonperturbative (multiparticle) formfactor 
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Figure 10.10: Ratio R = a^/ as a function of for p meson as a function of . Solid 
and dashed lines correspond to oscillator and suppressed Coulomb wave functions respectively 

of the proton. Therefore, in order to have a plausible ^-dependence, we have to introduce a 
certain "educated guess". On the other hand, our gluon density analysis already uses this 
formfactor, and therefore places some contraints on it. Finally, as the analysis suggests, the 
at low t, within diffraction cone, several other mechanisms modify the t dependence of the 
proton formfactor. If one parametrizes the differential cross sections at low t {t < 0.5 GeV^) 
by an exponential law with slope parameter b, 

^ (X e-^l*l , (10.15) 
d\t\ ^ ' 

one finds that various sources of t-dependence can be treated in terms of cotributions to the 
overall slope b approximately additively. 

In principle, there can be three sources of the non-zero contributions to the slope: the 
proton transition, the 7* — > y transition, and the exchange (the Pomeron propagation): 

b = bp^p + be^ch. + b^*^v ■ (10.16) 
First term appears in our calculations explicitly as an effective slope of the energy- and k^- 
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Figure 10.11: Effective exponents 6 of energy growth for the p meson production cross sec- 
tions. Solid and dashed lines correspond to oscillator and suppressed Coulomb wave functions 
respectively. The upper pair of curves ( at low Q^) corresponds to 5 calculated between 50 and 
75 GeV, while the lower pair is the result for 5 calculated between 110 and 150 GeV. 



independent dipole formfactor 

F(A2) = -2 ; A = 1 GeV ; ^ bp^p = 4 GeV-^ . (10.17) 

(l + AVA2) 

The second term is also introduced explicitly. It is responsible for the shrinkage of the diffrac- 
tive cone with energy growth. The third term appears from the accurate QCD treatment of 
the 7* — > y transition. It possesses a characteristic l/((5^ + my) shape, which leads to the 
scaling phenomenon mentioned earlier. 
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Due to the above blend of non-calculable soft and perturbative hard contributions, the 
/cj-factorization approach should not be expected to yield first-principles predictions for the 
absolute value of the slope parameter. However, several kinematical dependencies observed in 
the experiment must be confronted with the predictions. 

In principle, one can invoke several definitions of the effective slope parameter. However 
there is no siginficant difference among them. Below in Table 3 we compare results for three 
definitions of the effective slope: 



6(def.l) 



d\og{da/d\t\ 



d\t\ 



1 da 



■ 6(def.2) ^, , 
|t|=o '^^l^l 



6(def.3) 



1 



! da 



1*1=0 



calculated at several points. 



\t\) !\t\da' 

(10.18) 



Table 3. Various possible definitions of the effective slope and their values obtained from 
/ct- factorization calculations at three characteristic values of Q^. 



GeV^ 


6((lef.l) 


6((lef.2) 


6((lef.3) 





14.0 


13.1 


11.4 


2.2 


10.2 


9.6 


8.7 


27 


6.7 


6.3 


6.3 



In Fig. 10.12 we present our results for the slope parameter of the p meson production 
cross sections. The results for our calculations give somewhat too high values of the slope 
than the experimentally measured numbers (ZEUS [77] and HI [76]), but nevertheless they 
exhibit the right dependence. 

It is interesting to note that the slopes for aL and ar slightly differ, see Fig. 10.13. Note 
also that at small the calculations based on the Coulomb wave function give higher results 
for slopes. 

Fig. 10.14 shows a typical pattern of the |t|-depcndcnce of the differential cross sections 
in the region of small to moderate t (0 < \t\ < 1.5 GeV^). Oscillator wave function was used 
everywhere in this Figure. One sees that initial approximately exponentially decreasing of the 
differential cross section flattens at higher |i| as the process leaves the diffractive peak. 



10.2.5 Helicity amplitudes 

Our analysis exphcitly takes into account all possible helicity amplitudes 7*(A-y) V{\v)-, 
with A^, Ay = 0, ±1. Since the Pomeron exchange does not distiguish left from right, only flve 
independent helicity amplitudes survive. 

Fig. 10.15 shows the absolute values of the five helicity amplitudes against the momentum 
transfer squared \t\. Within diffraction cone one sees the charachterictic behavior of all the 
amplitudes. In the region of moderate \t\ one can observe diffractive dips, whose location 
changes from one amplitude to the other. Certainly, since we did not focus on large we 
cannot be sure that the dips are located exactly where we predicted. Still this picture shows 
a general pattern of the t behavior of the helicity amplitudes. 

It is clear that the presence of helicity flip amplitudes leads to the breaking of the s-channel 
helicity conservation. It is therefore interesting to check what is the magnitude of the helicity 
flip amplitudes. 
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Figure 10.12: Effective slopes b of the p meson differential cross sections within diffraction cone 
as functions of . Solid and dashed lines correspond to oscillator and suppressed Coulomb 
wave functions respectively. 

Fig. 10.16 gives the answer to this question. Here we show ratios of the hehcity flip to 
hehcity non-flip amplitudes 

\A{^t^V^t)\ \A{^t^Vl)\ \A{^o^Vt)\ , . 

\A{^t^Vt)\' |A(7o^Fi)|' 1^(70 ^ ■ ^ 

and whenever possible compare /c^-factorization predictions to the available experimental data, 
taken from [85]. 

Finally, we made predictions to the full number of spin density matrix elements for the p 
meson production and compared them with HI [76] and ZEUS [86] data. Results are shown 
in Fig. 10.17. 
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Figure 10.13: Effective slopes bi and bt of the p meson longitudinal and transverse differential 
cross sections within diffraction cone as functions of . Solid and dashed lines correspond to 
oscillator and suppressed Coulomb wave functions respectively. 

10.3 (j) mesons 

Production of (f) mesons has much similarity with p meson production. Therefore we will not 
provide as detailed discussion of the predictions as we did for p meson, but will rather show 
our direct predictions for the quantities that have been measured experimentally for the 
mesons. 

Fig. 10.18 shows the total cross sections of diffractive meson production as a function of 
data taken from HI [78] and ZEUS [79], [80]. Remembering the scaling phenomenon of the 
vector meson production cross sections, we should expect the picture similar to what happens 
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Figure 10.14: The transverse momentum squared dependence of the differential p meson pro- 
duction cross sections at two characteristic values ofQ^. 

in the p meson case. Indeed, we see at this figure a reasonably good description of the data 
for moderate Q"^ and a shght overshooting of our predictions as we shift towards small Q^. 

It is interesting to directly compare p meson and meson production cross sections taken 
at equal Q^. Fig. 10.19 shows our predictions for the ratio cr['-f*p <pp)/<y{l*P pp) toether 
with experimental data. If the scaling phenomenon holds, at higher the ratio is expected 
to approach constant value of 2/9, a tendency to do so is indeed visible at the figure. However, 
at smaller the ratio goes down. This behavior should be expected, for ta smaller Q'^ the 
ratio can be approximated by 



with rupeff < m^eff- 

One sees that our predictions for this ration agree with the data very well. Moreover, note 
a remarkable coincidence of the results based on the oscillator and Coulomb wave functions. 
The reason for that is of course the fact that we study here not absolute values of cross sections, 
but their ratios. This removes a significant part of ambiguity present in this or that specific 
choice of the wave function and reveals the features of the kt factorization approach in its pure 
form. 

Finally, Fig. 10.20 shows the experimental results for the meson density matrix measure- 
ment, published by HI [78]. Our predictions agree with the data well. 

10.4 J/tp and T mesons 

Fig. 10.21 shows the /cffactorization predictions for the total cross sections of the J/ip meson 
electroproduction as a function of compared with available data from HI [81], [82] and 
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Figure 10.15: The transverse momentum squared dependence of the five helicity amplitudes of 
p meson production at — 5.2 GeV^ . 

ZEUS [77], [83]. A reasonable agreement throughout the whole range is seen. 

The photoproduction cross sections are shown at Fig. 10.22 versus total energy W of the 
7P collision. One can see again that the /c^-factorization calculations give prediction with 
accuracy roughly of factor 2, uncertainty coming from the meson wave function. Still, the 
experimental data fit between the two curves. 

Fig. 10.23 shows ul — (Tt decomposition of the J/V' meson production cross sections, 
taken from [85]. An agreement at a similar level of accuracy is observed as well. The ratio 
R{Q'^) = o-i/ax is shown at Fig. 10.24. The data are taken from ZEUS and HI. 

The energy dependence of the T(ls) meson photoproduction is shown at Fig. 10.25. One 
observes that the experimental points soar about 10 times higher than the both curves. The 
origin of this discrepancy is not clear. At the one hand, the kt factorization predictions 
for the T production cross sections follow the scaling behavior in Fig. 10.1, while the T 
photoproduction data do not comply with this tendency. It should be noted also that a 
similar behavior is observed for other processes: that it, the production rates for T seems to 
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Figure 10.16: The dependence of the spin-flip to non-spin-flip amplitudes in the p meson 
production compared with combined HERA data. 

be higher than expected. 

10.5 Production of excited states 

In this section wc will give some of the most prominent features in the reaction of diffractive 
production of 2S and D-wave vector mesons. 

In Fig. 10.26 we show the ratios of the excited-to-ground state production cross sections 



r{2S/lS) 



ajjp ^ Vi2S)p) 
a{jp ^ V{lS)p) ' 



r{D/lS) 



ajjp -> ViD)p) 
a{jp ^ V{lS)p) • 



;io.2i) 



for p system and for charmonium. 

A brief look reveals that although the production rate both of 25* and of D-wave states 
are suppressedin respect to ground state, a radically different behavior of the suppress- 
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Figure 10.17: Experimental data on spin density matrix for the p meson production as function 
of Cf" compared with kt-factorization predictions based on oscillator wave function. 



ing factors r{2S/lS) and r{D/lS) is observed. This difference is due to distinct nature of 
suppression in these two cases. 

The suppression for 25* state production comes from much discussed node effect of 25* state 
wave function. Indeed, if one looks at the 25" production amphtude in the impact parameter 
space, then one product of two factors: the photon wave function multiphed by dipole cross 
section ■ a^ip and the radial wave function of the produced meson. At low Q^, ■ a dip 
is a wide function, peaked at a hadronic scale, that is precisely where the node of the vector 
meson wave function is located. This results in significant calcellation between contributions 
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Figure 10.18: Total cross section of the diff'ractive meson production as a function ofQ^. The 
kt-factorization predictions based on oscillator (solid lines) and suppressed Coulomb (dashed 
lines) are also shown. All calculation are performed for W — 75 GeV using DGD2002, Fit 1. 

from impact parameter regions smaller and higher than the node position. At higher Q^, the 
photon wave function shrinks, and the peak of \l/^ ■ adip shifts toward smaller r, that is away 
from the node. Thus, the effect of calcellation vanished. 

As for the D wave meson suppression, it arises from the angular part of the integrals. In 
fact, if the initial photon were built only of the S wave qq pair, then, due to ortogonality 
between pure S and D waves, there would be no D wave vector mesons at all. However, the 
spinorial structure of the photon coupling to the quark line does not correspond to the pure 
S wave, but contains an admixture of D wave as well. This D wave part leads to the D wave 
meson producton, and since it is rather small, the cross sections turn out suppressed as well. 
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Figure 10.19: Ratio of (j) meson to p meson total production cross sections as function . Solid 
and dashed lines correspond to oscillator and suppressed Coulomb wave functions respectively. 

The situation is basically the same both in the p system and in the charmonium. The major 
difference is the energy scale of the node effect suppression, which can be directly related to 
the mass of the corresponding meson. 

In the case of charmonium the fc^-factorization predictions are compared with the available 
HI data on ^^(25) production [87], [88]. A good agreement is seen. 

As vividly illustrated by Fig. 10.26, a mere measurement of the production cross section 
for high- lying states in the p spectrum (such as p'(1450) and p'(1700)) at high might be 
enough to get some insight into the spin-angular structure of these states. However, even 
more dramatic distinction between S and D wave states is offered by the Ry — crL{V)/aT{V) 
measurements. 
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Spin density matrix for (|) meson 




Figure 10.20: Experimental data on spin density matrix for the meson production as function 
of compared with kt-factorization predictions based on oscillator wave function. 

Fig. 10.27 shows rations R{IS), R{2S), and R{D) for the p system. We see that at higher 
Q"^ the following hierarchy takes place 

R{1S) > R{2S) » R{D) , (10.22) 

with about one order of magnitude difference between the 15" and D-wave states production 
rates. At smaller Q^, R{2S) exhibits very characteristic spectacular wiggles. Starting from 
zero at = 0, it springs to a local maximum of at ^ GeV^, then rapidly fall down, 
and then rises again. This behavior is the manifestation of the node effect as well. Indeed, 
since the structure of the transverse and longitudinal photons is not identical, the strongest 
cancellation takes places at different values in these two cases. In particular, the curve 
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Figure 10.21: Total cross section of the diffractive J/ip meson production as a function of 
. The kt-factorization predictions based on oscillator (solid lines) and suppressed Coulomb 
(dashed lines) are also shown. All calculation are performed for W — 75 GeV using DGD2002, 
Fit 1. 

suggests that as start from — Q and we slide along scale, the strongest cancellation 
occurs in ax earlier than in o"l. 

A word of caution shoud be said now. The exact position of the node in the 25* state wave 
function cannot be predicted accurately and depends on the particular Ansatz for the wave 
function chosen. This means that the wiggly shape of R{2S) in Fig. 10.27 — with the 
positions and the values of local extrema — should not be regarded as an accurate prediction. 
In fact, one cannot even be sure that the non- monotonous behavior of R{2S) actually 
takes place. This curve simply illustrates what kind of effect can be expected. The specific 
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Figure 10.22: The energy dependence of the total J/ip meson photoproduction cross sec- 
tion. The kt-factorization predictions based on oscillator (solid lines) and suppressed Coulomb 
(dashed lines) are also shown. All calculation are performed for W — 75 GeV using DGD2002, 
Fit 1. 

profile of this curve is not of course predictable, as long as we rely on simple vector meson 
wave functions. 

An immediate conclusion from here is that we found an observable that is extremely 
sensitive to the minute details of the vector meson wave function — a powerful tool that 
would help distunguish among various models of the vector meson structure. 

The t-dependence of the 25* state production also shows remarkable features originating 
from the node effect. Fig. 10.28 shows how the t-profilc of the differential cross section changes 
over a narrow region. The strongest cancellation due to the node effect takes place at 
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Figure 10.23: Decomposition of the J/ip production cross sections in (Tl and (Tt 

= , |t| = 0. When we increase \t\ or Q^, the strength of the cancellation diminishes, and 
the differential cross section grows. 

Finally, in Fig. 10.29 we show the spin-angular density matrix for the 15, 25, and D states 
in the p system. A dramatic difference among the profiles both SCHC and SCHNC matrix 
elements is seen. Note that some of the matrix elements are even of opposite sign ofr S and 
D wave states. 
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Figure 10.24: Ratio R = ai/crT as a function of for J/ip meson as a function ofQ^. Solid 
and dashed lines correspond to oscillator and suppressed Coulomb wave functions respectively 
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Figure 10.25: The energy dependence of the total T(IS') meson photoproduction cross sec- 
tion. The kt-factorization predictions based on oscillator (solid lines) and suppressed Coulomb 
(dashed lines) are also shown. All calculation are performed for W — 75 GeV using DGD2002, 
Fit 1. 
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Figure 10.26: The behavior of the ratio of excited to ground state vector meson production 
cross sections in the case of p system and charmonium. Dashed and dotted lines represent 
a{2S)/a{lS) and a{D) / a{lS) ratios respectively. In all calculations, oscillator wave function 
was used. 
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Figure 10.27: The behavior of the (Tl/(Tt ratio for IS (solid line), 2S (dashed line) and 
D wave ( dotted line ) states in the p system. In all calculations, oscillator wave function was 
used. 
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dcr(7'*V^P(2S))p/dt, nb*GeV'' 




Figure 10.28: Strong forward dip for 2S state in p system: the change in t profile of the 
differential cross section wihtin narrow region. 
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Figure 10.29: The behavior of the spin ensity matrix for IS, 2S, and D wave states in the 
p system. 
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Chapter 11 
The a^p puzzle 



The detailed comparison of the numerical results for p meson production obtained in the 
fcf-factorization approach with the experimental data shows that the scheme used fails to 
reproduce the correct behavior of the transverse cross section ut, which was dubbed by 
us as (Jt puzzle. In this chapter we explore two possible causes of this mismatch. Namely, we 
explore the effect of color Coulomb interaction of qq pair and show that it leads to ctt increase 
in respect to ai. Then we analyze the issue of possible S/D wave mixing in the p system and 
try to obtain a better description of the ai/cTT- 

As we will see, the first method might turn our the remedy for the ctt puzzle, but even an 
accurate formulation of this method forces us to go beyond the lowest Fock state, thus making 
only rough estimates of the effect possible. The second mechanism will be shown indeed to 
decrease the (Tl/o't ratio at higher down to an acceptable values, but this happens at the 
expense of too low values of ai, not increased values of ctt- Therefore, this mechanisms does 
not provide the solution to the ax puzzle, at least in its pure form. 

11.1 The Coulomb tail of the wave function 

As we discussed before, the major source of ambiguity in our approach is the vector meson 
wave function. Without accurate knowledge of the radial wave function, one is bound only to 

an educated guess. 

In our treatment we stuck everywhere to the soft wave function Ansatz, that is, we assumed 
that any integrated involving the wave function would be saturated by small-]? region ~ 
1/Rv)- This implies that at large-p the wave function must vanish fast enough, or to be 
precise, faster than any power present in numerators of the above formulas. Since the 
highest momentum power is P, which appears, for instance, in the transverse amplitudes, we 
conclude that the "softness" of the wave function means that 



The Gaussian wave function Ansatz used everywhere above certainly satisfies this requirement. 

It turns out, however, that the above requirement is violated when the true short-range 
Coulomb-like interquark potential is taken into account. As well known, if the color Coulomb 
potential were the only source of the quark-antiquark interaction, the resulting Coulomb-type 




(11.1) 
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wave function of the vector meson would look like 



In reality, the interquark forces are much more complicated. However, without any need to 
know the precise form of interquark forces, one can assert that at short distances there must 
be some resemblance of the Coulomb-like wave high-p^ tail. 



11.1.1 The strategy 

Certainly, there is a multitude of approaches that claim to account for the color Coulomb 
interaction at short distances. One of them would consist in obtaining an accurate numeric 
solution of a given potential model. If one prefers to used a fully QFT-based approach, one 
must deal then with a set of diagrams. Ideally, one should start with free quarks and then, by 
taking their interaction into account, arrive at the physical vector meson with a (presumably) 
uniquely defined qq wave function. 

Realizing such a program is still unresolved task and attempts to accomplish it would 
certainly go far beyond this study. Moreover, even if we were capable of doing it, we would 
still run into interpretational difficulties, at least in the framework of our scheme. Indeed, when 
calculating the diagrams, we will see the higher Fock states intrude into our qq description of 
the vector meson. So, even if wc still wanted to follow only quark-antiquark distribution in a 
meson, we will be forced to switch from the wave function to the density matrix description. 
Another problem would be the presence of higher angular momenta of the qq pair due to 
gluon radiation. This will leave no room for our simple S wave/Dwave description of the 
vector meson. Finally, the qq wave function (or density matrix) will have no unique, process 
independent definition. Indeed, in order to preserve the gauge invariance at any given order of 
perturbation theory, we will have to include corrections to the rest of a diagram (the kernel) 
as well as corrections that entangle the kernel with the wave function. As a result, the two 
particle irreducibihty of the process will be lost. 

The conclusion is that we cannot expect a reasonable answer to the question "how to 
account for short range Coulomb interaction" and still stay within the framework of our 
approach. 

Nevertheless, the mere formulation of this question does not force us to go beyond our 
lowest Fock state, two-particle irreducible approach. So, it can be reiterated as "The Coulomb 
tail of the wave function must be there; given that, can we develop a reasonable understanding 
of its impact on dynamics of the vector meson production?" In this formulation, we now ask 
for a QFT-inspired model that would produce an estimate of the hard tail without forcing us 
to run into problems just mentioned. 

We suggest the following procedure that would satisfy this need. We start with the soft 
wave function Ansatz and perform a sort of evolution procedure, that will generate the hard 
tail. This evolution will basically consist in allowing for gluon exchange between the quark 
and the antiquark. The result will be interpreted as the hard piece of the vertex function F, 
and eventually as the hard piece of the wave function itself. 
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11.1.2 The quantum mechanics of the Coulomb tail 

Suppose that our potential is a sum of an oscillator and Coulomb potential, the latter being 
a "perturbation" and bearing an intrinsic small parameter: 



Voscip) = ^S"ip) , v^ouiip) = --r—T^s{p')CF . (11.3) 



v{p) = Voscip) + Vcouiip) ; 

P P ^ H'G 

Here coupling constant is supposed to be the small parameter. Effective parameter /iq accounts 
for confinement. 

The total wave function is represented as sum of the soft part ips = c-exp(— A;^a^/2), which 
is the solution of Schrodinger equation with oscillator potential only, and the hard part iph, 
determined via: 

- Mp) = J dkVcoui{k)MP - k) (11.4) 

Let's for a moment neglect running of the coupling constant. Then, 

V'.(p)=g^-^|^-/(P), (11.5) 



SO that 



Here 



m - exp (^-^ j + ^ . . m . (11.6) 



/(p) = / £^ cxp = 27r r^dk^log 

^ ' J [p-kf + fil \ 2 J Jo 2pk ^ 

/■°° k'^dk e 

2ti 



jp + kf + 1^1% 
{p - kf + 



k^a^' 



exp -■ 



/■°° k^dk 6 2—62 

Jo k^ + i^il k^^ ■ ^ ' 

This function cannot be evaluated exactly. One can, however ask for its the asymptotic large- 
p behavior, which can be evaluated directly from (11.4) by replacing k ^ p — k and taking 
Vcouiip - k) ^ Vcouiip) out of the integral: 

Mp ^ oo) = 47rC^a.^ • M = 0) . (11.8) 

Note, however, that one should not use this simple analytic form because it spoils the large 
distance behavior of the wave function. The wave function (11.8), even when regularized as 
p2 ^ p2 /ig, leads to exponential decrease of the wave function as oc exp(— r/xc), while the 
honest wave function exhibits gaussian decrease oc exp(— r^a^/2). 



11.1.3 Derivation of Aip, 



Coul 



The diagrammatical representation of this procedure is show in Fig. 11.1. Instead of imple- 
menting soft vertex function T^ofv properly normalized 

r?o*a^ = ^(r^o/t + An , (11.9) 
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Figure 11.1: Diagrmmatic representation of the effect of Coulomb tail at small distances 



where 1/yN factor accounts for the proper normahzation. The Coulomb correction AF" 
reads then 

^ ' y (27r)4 ^ {ll - w? + ie){ll - w? + ie)[{ki - hf - 11% + is] ^ ' 

The Sudakov decomposition of the momenta of the initial and intermediate particles' reads 

1^2 = {l-y')p'^ + il-z')q"^-l>l; 
{h-hr = {y-y')p"' + {z-z')q"^+{k^-l^r. (11.11) 

We already know from the above analysis that one of the initial particles is on-mass shell. 
Therefore, one of the two conditions 

P + P + 

y = o— or y = 1 



zruy (1 — z)my 

is fulfilled. 

The numerator constitutes the most challenging part. Certainly, there is absolutely no 
problem to perform gluon index summation: 

-7^(-f2 + m)7"(/i + m)7^ = 2mV - 2/i7"^^2 - ^{h - hT ■ (11-12) 

The problem however is that this spinorial structure cannot be reduced to the scalar multi- 
plicative renormalization of the original structure 7". Moreover, after integration the spinorial 
structure cannot be expressed in terms of S and D waves only. Finally, the resulting spinorial 
structure will lead to not gauge invariant correction to the wave function. 

In order to go as far as we can in trying to obtain a reasonable estimate for the Coulomb tail 
and avoid trouble, we propose the following procedure. First, we look at the transverse vector 
meson production. Since the vector meson polarization vector possesses only transverse 
components {V"), only transverse terms 

2m?T + ^l^'^l- 2[z\l - y')p'TQ' + ^'(1 - z')q'TP'] - sf" . (11.13) 

will survive from the whole expression (11.12). Then, in order to simplify life, we will proceed 
for the p meson as if it were a heavy quarkonium. That is, we will neglect P in comparison 
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with and put = |. Finally, the last stroke will be to put niy = 4m^ in numerator, so 
that y' can be replaced by 1/2 as well. The result of this procedure will give 

-7^(-/2 + m)T{h + m)7^ ^ 4m V, • (11-14) 

Finally, we postulate that the same scalar renormalization holds for longitudinal vector mesons 
as well. 

The denominator is calculated in the same way as before. We note that the numerator 
in our calculation does not depend on y', and therefore, the pole structure on the y' plane 
comes from the propagators only. The position of poles, however, changes with z': ai z' < 
or z' > 1 all poles lie in the same half-plane, and thus by closing the contour we can nullify 
the contribution of these z' regions to the integral. Then, the interval < z' < 1 naturally 
breaks into two parts: z' < z and z' > z intervals. In the former region we prefer to close the 
contour from below and take the residue at the position of pole 1. The two other propagators 
become 

(l-z')(m''y-M")+is; 



pole 2 
pole 3 



{z - z') yrriy - 



J' 



:il.l5) 



In the second z' region we close the contour from above and take residue at = 1 — (P + 
m^)/[(l — z')my], so that the other poles become 

pole 1 z'{my - M'^) + ie ; 



pole 3 



-(z'-z){M''-m'y)-iz'-z)[ym'y 



m 



Note that in both cases the familiar M'^ — my appears in the denominator, which fuses with 
TsoftiM'^) to produce ipsoftiM'"^). Thus, the result of the y' integration is 



with 



Ai 

A2 



dy' 



1 . . 1 



e{z - z') 



-{z - z') 



Ai + (A; - ly + //^ 



+ e{z' - Z) 



{z' - z){M''^ - ml) 



:il.l7) 



ki on mass shell; 



Ai 
A2 



,z — z 



P- 



1-z 



1-2' 



+ {z' - z){M'^ - ml) 



k2 on mass shell. 



The analysis of the positive definiteness of the gluon propagator reveals that in one case 
the positiveness indeed holds due to inequality 

2 

/ / ^/ \ 

\2 



■^-^k'^ + ^-^t' + ik-iy 

z z' 



>0, 
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while in the other case the gluon pole can in principle arise in the allowed kinematical region. 
Such a pole would correspond to the situation when both quark and gluon are simultaneously 
on mass shell. However, in the case of p mesons such configurations are avoided due to small 
enough mass of the vector meson [my < 2ij,q). 

— * 

So, the expression for Ar(2;, k) takes form 



X 



47r2 



dz' 



z'{l - z') 

1 



Ai + (A; - /)2 + /^^ 



+ - z) 



. (11.18) 



With Ar{z. k) calculated according to (11.10), we can now construct the hard part of 



the wave function ip{z,k). However, the straightforward answer Sip = Ar/(M^ — my] 
not be a satisfactory option due to unphysical pole at = my. 



will 

Wc remember that when 
constructing the soft wave function, we forced the corresponding vertex function Tgoft to have 
zero at = my in order to cancel the unphysical pole. Here we do not have such a freedom 
in manipulation with AF, which is, by the way, always positive. However, we can again tune 
the soft vertex function Tgoft so that the entire expression Tgoft + AF does have zero at the 
required point. This can be achieved by shifting the node position of the soft vertex function 



Tsoft{z,k^)^i^soft{M^) 



,,2 2 Ar(^,A;^) 
-ml- ^ ' °^ 



;il.l9) 



where hpIL — ui"^ xhe total wave function will then turn into 

z(l—z) y 



i'total{z,k'^) 



- ml - 



Ar(^,fc-^) 



V ipsoftirriy) 



+ Ar{z,k^ 



M2 - ml 



V 



= lljsoft{M^) + 



1 



M2 - 



Ar(.,P)-Ar(..«sf^, 

Wsoft[TTly)_ 



(11.20) 



The wave function constructed in this manner is regular at — my and does have the 
expected large behavior. 



11.1.4 The large asymptotics of ctt' analytical result 

The honest integration in (11.18) cannot be done analytically, but in order to grasp the 
asymptotics, a simple estimate can be performed. 

For k"^ ^ P, fiQ, m^ the denominators turn into k'^z' /z and k!^{l — z')/{l — z) for on 
mass shell initial quarks 1 or 2 respectively. If we then suppose that the soft wave function is 
given by the Gaussian Ansatz with Rym^ 3> 1, we can put z' — > 1/2 and do the integration 
completely. The answer reads 

AT{z, P) = C,'^ ■m'-'^-j dHiPsoftiM'^) = ^ C,'^ ■ Ama^ ' i^] ■ (11-21) 

7r2 mk^ J f.2 \RvJ 
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We see that the hard tail of the wave function AF falls off as at large P, in accordance 
with our expectations. 

The Coulomb tail of the vector meson wave function being finally obtained, we are now 
ready to estimate the expected change in the (Tl/cft ratio for high-Q^ production of p mesons. 
Before we take a look at the numerical values, let us make an analytical estimate of the Q"^ 
asymptotics. We plug the hard tail of the wave function A'^(2;,A;^) into the leading twist 
expression for At^Ti which we write as 



ImAn 



dz 



z{l-z) 



d^kiptotai{z,k^ 



2k2 



{m' + 2[z'' + {l-zf\k^} 



47r2 



The integral 



•JtotaliQ'^) 
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zy]k^ 



z{l — z) 



[z{l - z)Q^ + K^ + k 
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is naturally split into JsoftiQ'^) and JhardiQ^)- The former is calculated in the heavy- 
quarkonium approximation with Gaussian wave function and yields 



/■ 4 

JsoftiQ^) = / d^kj/jsoft 

J my 



4m^ 



V 



16mv 
— ^(^^ 



27^ 



The hard piece is calculated as 



JhardiQ ) 



Jq z{1 — z) J/j. 



= 2mya^ • 
The ratio Ahard/Aoft is then 



as{k'^ 



2n_ 



'2n_Y Cf16 



C dzl 
Jo 



z^ + a- z?] 



(11.23) 



q;,(A;2) (11.24) 



AhardiQ'') 4C, 



AsoftiQ'') 37r 



fQV^dk^ ,-^0, 64, 



Q^s(/^^) 



:il.25) 



This expression represents the asymptotic large-Q^ behavior of the Ahard/Aso/t ratio. As it 

could be expected, it is governed by log ""i^^, , which rises with growth, but extremely 

slowly. Within the experimentally studied region of (not higher than 100 GeV^), this 
logarithm, which is supposed to be large, is still less than unity (basically 0.5 0.7). This 
indicates that the asymptotic regime is not reached at HERA, and the estimate just derived 
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seems to be irrelevant to the real magnitude of the hard tail effects. Nevertheless, this estimate 
tells us that no abnormally suppressing factors appear in our problem. 

The numerical analysis of the Coulomb tail impact on the production amplitudes will be 
given elsewhere. 



11.2 The S/D wave mixing 

The presence of tensor forces in two-fermion bound states is always a natural consequence of 
relativistic corrections. They lead mixing of states with definite angular momentum, which in 
the case of a vector particle translates into S wave - D wave mixing. 

The most famous example is provided by the Breit potential in a hydrogen atom. Such 
forces are present in a deuteron, where they lead to about 5% admixture of D wave state, and 
they can be present in vector meosns as well. Review [89] provides further examples. 

How are we going to mix S and D wave? When calculating the normalization or decay 
constant for a vector meson, we observed a two-fold dependence on the vector meson state: 
in the spinorial structure and in the particular form of the radial wave function. In all cases, 
the rest was absolutely insensitive to the vector meson state. Therefore, the two quantities we 
should mix, are S^i/jsip"^) and P^'0£)(p^): 

|r^z^(p2)) = \S^Mp')) COS0 + |P^Vd(p') sin^) . (11.26) 

If mixing angle 6 is constant, then this vector meson state satisfies automatically the normal- 
ization condition: 

{r^ip\r,ip) = {SM^^^i^s) cos" e + (p^v^dIi^.V^d) sin^ e^g^,. (11.27) 

There can be several types of S/D mixing. First, we can mix D state with IS", 25* or even 
higher S states. It is not clear a priori which mixing should be stronger. Therefore in our 
analysis below we accounted for possible IS/ D and 2S/D mixing on the equal foot. Second, 
the mixing angle 6 can be an explicit function of p^. This case is a bit more complicated but 
nothing seems to rule out such possibility. 



11.2.1 Cooking up 7^ vertex 

The first issue we want to elaborate is how to construct the 7^^ vertex from accurate 5^ and 
Vj^ spinorial structures. This question arises when one tries to check how the naive u'jf^u 
treatment of the quark- antiquark- vector meson vertex differs from the accurate one. To do 
so, in the expression 



^Mp') = (7^ - M^+2m ) "^^^P'^ cos^(p2) + (p^^^ + (M + m)p,) Mp') sme{p') . 

(11.28) 

all terms with must be banished by properly adjusting mixing angle 9 — ^(p^). Assuming 
for simplicity, that we use oscillator ansatz wave functions with all equal, one obtains for 
IS/D mixing 



2\ ^ v30 

tan6/(p ) = - + 2m)(M + m) ^ " ~ (M + 2m) (M + m)a^ ' ^^^'^^^ 
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A very rough estimate for p system gives tan^ ~ 0.2. For 2S/D mixing one has 



V45(l-2pV/3 

ts^neip^) = -ji- — ; ,,,, — 11.30 

' {M + 2m){M + m)a^ ^ ' 

This time the mixing angle even swings with growth. After terms are canceled 
in (11.28), the residual expression is proportional only to 7^ and has slightly modified 
dependence via explicitly p^ dependent mixing angle. However, this modification is completely 
inessential. Moreover, for 2S/D mixing this modification does not shift the position of the 
node. Therefore, it seems that the impact of this specific mixing on the radial wave function 
is minimal. The magnitude of this mixing can be also described with only one number — the 
integrated mixing angle: 

= (11-31) 

(cos^6'(p^))5 

where {■■■)s,d mean the normalization integral for S/D state. For \S/D and 2S/D types of 
mixing, ia.n^ 6* is equal to 0.027 {6* = 9.3°) and 0.062 {9* = 14°) respectively. 




Figure 11.2: The impact of IS/ D mixing in the p system on the oj^Igt ratio for the lowest 
energy state production. 



145 



11.2.2 Impact of S/D wave mixing on p meson production 

For the simplicity, we will restrict ourselves to the constant-angle variant of S/D wave mix- 
ing. Since our primary motivation here is to check how mixing alters the (Tl/(Tt ratio in p 
production, we will focus only on p meson. 

Since a generic p meson production amplitude is linear in the meson wave function (un- 
derstood in its complete sense, i.e. spinorial structure times the scalar wave function), one 
immediately has 

COS 9 + Ad ^^T^O . (11.32) 

With this expression in mind, we can now predict the effect of S/D wave mixing by simply 
looking at the pictures for pure S and D wave states production. 



Oi vs, aj decomposition 




1 10 1 10 

Q\ GeV' Q\ GeV' 



Figure 11.3: Changes in profiles of aT and aL cross sections for p system ground states 
caused by IS/D mixing. 

Fig. 11.2 and Fig. 11.3 show the changes in the longitudinal and transverse p meson cross 
sections caused by the IS/D mixing. One sees that at (pmix > grows, the ratio ai/crx 
decreases. Fig. 11.2 implies that (p ~ 0.7 would do the best job in describing the data points. 
Note however that the shape of oi^/ut ratio remains the same and does not significantly 
flatten, which would be needed for a better description of the experimental points. Even more 
problems come from separate analysis of 0"^ and cross sections. One sees that decrease of 
(Jl/i^t arises not from making higher, but at the expense of making ai significantly lower 
than the data. Thus, although the ratio gl/^^t can indeed be corrected in a simple mixing 
scenario, the cross sections themselves still deviate at large from the measured values. 
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This analysis lead us to a conclusion that the gt puzzle in p production still persists in our 
approach. 

Although the S/D wave mixing failed to completely resolve the ut puzzle, it is still an 
interesting issue on its own. In constrast to all previous calculations of the vector meson 
production amplitudes, our approach allows for a clear formulation and detailed analysis of 
this phenomenon. 
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Chapter 12 
Conclusion 



In the present work we formulated the fcr factorization approach to the calculation of the 
diffractive production of vector mesons in DIS. Since this approach is organically linked to 
the BFKL dynamics, no requirements on the value of + my was placed, which allowed us 
to investigate the production cross sections from the photoproduction limit up to highest 
values attainable in the experiment. 

When describing the vector meson, we limited ourselves to the lowest Fock state only, 
that is we treated a vector meson as a bound state of quark- ant iquark pair only. Being of 
course an approximation, this limitation allowed for the strict construction of the pure 5'-wave 
and D-wave states of the qq pair sitting inside the meson, which was then applied to the 
corresponding production amplitudes. 

On the other side of the reaction, we related the production amplitudes to the unintegrated 
gluon structure function of the proton. We undertook an extensive study of the differential 
gluon density, which included its first ever extraction from the experiment and detailed inves- 
tigation of its properties. This was then used in the vector meson production calculations and 
yielded reliable numerical results. 

Here we give a detailed list of main results of the present work. 

1 . Differential gluon structure function of the proton 

• First ever extraction of the unintegrated gluon structure function and casting it 
into a form of simple and ready-to-use parametrizations; 

• detailed analysis of the soft-hard decomposition of various observables and the 
impact of soft-to-hard diffusion phenomenon on them; 

• observation of self-generated two-reggeon-like structure of F^^ predicted by kt- 
factorization; 

• showing by means of analysis that the high-precision experimental data on 
suggests rather strong separation of soft and hard parts of the unintegrated gluon 
density. 

2. Diffractive vector meson production: analytical study 

• Developing for the first time a theory for vector meson spin-angular coupling and 
using it consistently in the derivation of the meson production amplitudes; 
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• observing remarkably different Q^-dependencc of 5* and D wave type amplitudes, 
providing thus a way to discern S and D wave states that are indistinguishable at 
e'^e~ colliders; 

• at large Q^, observing a dramatic role of higher twist contributions to the D wave 
vector mesons, which even forced sign change for L — > L amplitude and led to 
non- monotonous (Tl/(7t ratio; 

• observation of very large helicity violating effects for D wave vector states, which 
do not get suppressed even in the case of heavy quarkonia; 

• confirmation of the soft dominance of the double spin flip amplitude in the case of 
both S and D wave states; 

• establishing the borders of our approach during discussing the hard Coulomb tail 
of the wave function. 

3. Diffractive vector meson production: numerical study 

• Showing that /c^-factorization approach leads to a good overall agreement with 
availble experimental data on all types of vector mesons. Namely, we showed that 
the overall shape of dependence, energy growth, ^-dependence, picture of s- 
channel helicity violation observed in p, 0, and J/-^ mesons production are in good 
agreement with /cj-factorization prediction. 

• Recognizing that two particular issues — the transverse cross section aj- for p 
mesons at higher and the magnitude of T(ls) state — still remain unresolved. 
The causes of this discrepancy and the resolution possibilities are discussed. 

• Predicting many previously unknown features of the excited vector meson produc- 
tion reactions, including opposite signs for the largest spin flip amplitude T — > L 
for S and D wave vector mesons and dramatically different (Tl/ut ratio; 

• investigating S/ mixing induced phenomena; indicating that large S/ D wave mix- 
ing in p system can be the origin of ctt puzzle. 

Since the work contains a large number of predictions on the observables that have not yet 
been investigated in experiment, I hope that the thesis will serve as a guide to the directions of 
future experimental research. The dramatic differences between S and D wave vector meson 
production predicted in this work demand confirmation and offer a novel way to study the 
structure of hadrons. 

Finally, in this work I intended not merely to list the results, but also give a detailed 
and pedagogical presentation of all steps. In particular. Part I, where I introduced the kf 
factorization scheme. Chapter 4, where I construct the theory of the vector meson spin-angular 
coupling, and Appendices contain very detailed intermediate calculations, so that everyone can 
follow the entire line of derivation. To this end, I hope that the present text can be used as a 
means of learning as well. 
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Appendix A 

Denominator evaluation: details 



Below we give a detailed denominator calculus that is used in all /cj-factorizatoin processes. 
The major guideline in this derivation will be the analysis of pole positions and setting some 
of the propagators on mass shell by taking appropriate residues. 



A.l Forward Compton scattering 
A. 1.1 The s-channel diagram: all details 

We first start with the case of forward virtual Compton scattering amplitude. The integral to 
be calculated is 



J dy dz da df3 



[(g — kY — rm? + ie\ [A;^ — rn? + ie] [A;^ — rm? + it] [[k + kY — rn? + ie] 
^ (A.l) 



[k^ — jj? -\- ie]^[(p — kY — w? + ie] 
With Sudakov's decomposition 

S S 

one can rewrite all 7 propagators as: 

k'^ + m^) + ie , 



' 2 



(1) 




kf - 


+ ie - 


(2), (4) 




+ ie - 




(3) 


{k + 


\2 2 

K) — m 


+ ie 


(5), (6) 


k'- 


1^ -\-ie — 


: aj3s 


(7) 


{p- 


\2 2 

K) — m 


+ ie ■■ 
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Let us analyze the position of propagator poles (namely, the sign of ie) in complex y plane as 
functions of z and (3. One gets 

1 y^-^-- - + - ; 

s 1 — z)s 1 — z 



(2), (4) y 



k!^ + ie 
zs z ' 



(3) y = -a + ^ ,^ , —s . (A.3) 



For the poles on a plane we get 

(3) a = -y + 
(5), (6) a 



{(5 + z)s z + p 

{k + Ky + m'^ ie 

{(5 + z)s ^ + 

-*2 1 2 



/3s /3 ' 

7 a = 1 + ^ + — . A.4 

Since the function we integrate (namely, the product of propagators) is an analytic function, 
which vanishes at large \y\ and large \a\ fast enough, we can switch from (— oo, +oo) integration 
to the complex plane contour integration both on y and a planes. This is done by adding 
half-circle of infinite radius either on the upper or the lower half planes. What we now have 
to done is just to trace how many poles we have inside the contours. 

As clearly seen from (A.3) and (A.4), the position of the poles on y and a planes depends 
on the values of z and /3. If values of of z and (3 are such that no poles happen to fall inside 
the contours, the integral turns zero, and corresponding {z, [3) region does not contribute to 
the total integral. So, our task now transforms into searching for such (z, P) regions that both 
y and a contour integrations result in non-zero values. 

A convenient way to perform this analysis is to do it graphically. Figure illustrates the 
sign of ie in 1/ = . . . ± ie (left pane) and a = . . . ± ie (right pane) for all values of (2;, 13). Each 
shaded line here corresponds to a propagator from (A.3) and (A.4), the shaded side indicating 
the half plane with positive value of ie. For example, propagator (1) from (A.3) has positive 
ie when z <1. 

If we now take a closer look at each of possible regions (-2, /3), we see that the only (z, /3) 
pairs that result in a non-zero expression lie inside a triangle. For convenience, we break up 
the whole area into three sub-regions A, B, and C. For each of these regions, we get only one 
pole inside y and one pole inside a contours: 

pole (1) for y; pole (3) for a. 
pole (1) for y] pole (7) for a. 
pole (3) for y; pole (7) for a. (A. 5) 
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Thus, the integrals over y and a turns into taking the residues at certain poles. Let us 
write the result of this procedure for each of the three regions. 



Region A. 

Here we get 



y = 



^ _ Q + ^ {k + kY + 

s (1 — z)s ' s (1 — z)s {z - 



We know that /3 < 0, so we explicitly switch for more clear notation 
propagators turn now into 



(A.6) 
The other five 



(2), (4) 
(5), (6) 

(7) 



1 



1-z 



[k^ + m^-z{l-z)Q'']+i4; 



->2 2 
-K, — H 



->2 2 I 

—K — m + 



F + {k + kY + m? 



1-z 



+ 



k +171 {k + k) + m 



1-z 



z-\P\ 



m 



+ +ie.(A.7) 



Here symbol ^ means that the rule for pole passing is not needed since the propagator has 
definite sign for the physical values of the parameters. 
Thus, the honest result of the integration equals 



dy dz da d(5 



1 



X 



propagators 
1 



27ri\2 /-i 



] I dz r d\(3\ — 

J Jo Jo z — 



1-z 



1 



[F + m2 + ^(l -^)Q2]2 



2 _i_ fc2+m2 _j_ (fc+K)2+m^ 

z-m 



g2 + ^i±^ + 
1 



+ mP' — {m? + \l3\s) 



s 1-z z-\B\ 



(A.8) 



le 



Now comes the last step. The result (A.8) has both real and imaginary part. Since we are 
hunting for the imaginary part only, we extract it by using 

V.P. f 1) + i7^5{X) 



X-ie 



X 



and kiUing the /3 integral with the means of S function. The value of P is 

\(3\ «^ - « 1 

and therefore it can be neglected everywhere else in (A.8). The result reads: 



Im 



dy dz da dp) _ 47r^ ^ 1 - z 1 1 

propagators } ~ Jo ^ z [k^ + m'^ + z{l - z)Q'^Y + 



(A.9) 



(A. 10) 
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Of course, the z integration limits in (A. 10) should not be understood literally as and 1. 
In fact, when z is close enough to or 1 (by a typical value of ~ w? /s), (3 can no longer be 
neglected in comparison with z. Therefore, the integrand in (A. 10) does not have the correct 
small z and small 1 — z asymptotics. Nevertheless, since the physical quantities will have 
regular z ^ Q and z \ behavior, the impact of this difference will be suppressed by ~ m^/ s 
simply due to small integration measure. 
Regions B and C. 

Although region B has much smaller area than region A, one cannot guarantee a priori that 
the third pole, which produces the imaginary part of the integral, does not happen to lie inside 
it. However, one can check that it is not the case. 

One extracts first the values of y and a according to (A. 5) 

y = Ts T' — ~2 — IT- ^-^^^ 

s (1 — z)s — ps 

We remind that < /3 < m? / s, so that no singularity arises here. The remaining propagators 
are 

(2), (4) ^ -J—[k^ + m^-z{l-z)Q^]+ijt- 
1 — z 

(5), (6) - 

(3) ^ -[(k + ny + rn']-(Q' + ^:^]{z + P)-^^{z + P)s+zAA.12) 

\ 1 — z I — ps 

We see that throughout the whole phase space of the remaining kinematic parameters all the 
propagators keep definite sign, which makes the answer purely real. 

The similar picture occurs in region C as well. Therefore, regions B and C do not contribute 
to the imaginary part of our integral. 

A. 1.2 The other three s-channel diagrams 

The evaluation scheme just described can be applied to the other three diagrams. There will 
be slight modifications in quark propagators, but the strategy is unchanged. 

For example, consider diagram (b) in Fig. 5.2. The only difference is the expression for 
propagator (4): 

(4) = {q-k-Kf-m'^ + ie = {l-z-(3)(^-^-y-o}j-[{k + K)'^ + m'^] + -e; 
_ Q"^ {k + K,y + 171^ ie 

^ y_-a-- ______ + ____; 



If 



^ a = -y-^-^- +--, n- (A-13) 

s [1 — z — p)s 1 — z — p 

The graphic representation of the pole position on y and a planes is shown in Figure. Again, 
one can check that the main nonzero contribution comes from the same region A and the same 
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residues. So, if one again calculates the imaginary part of the integral, one finds 
dy dz da d(5 \ 



Im 



(/ 



propagators 
47r2 



/•II 1 
- Q / d.^— =^ ^(A.14) 

S3 Jo + r7l2 + ^(1 - ^)Q2] [(^ + + ^2 + ^(1 _ ^)g2] [«;2 + ^2]2 V 

Again, one can make sure that the small regions B and C do not contribute to the imaginary 
part as well, although the proof might not appear as clean as before. Namely, one can find 
that regions B or C will contain a propagator (namely, propagator (4)) that can change its 
sign, and therefore will contribute to the real part. However, this can happen at abnormally 
small 



?2 



m 



s 



Therefore the contribution of this configuration is 1/s-suppressed, and we cannot take it into 
account at our level of accuracy. 

Finally, the answers for the remaining two diagrams are 

47r2 /-i 1 1 

I dz— =i -r=;z rr- diagram C; 

Jo [(k + K)^ + m^ + z{l- z)g2] [A;2 + ru^ + ^{1 - z)Q^] + A*']' 

/ dz^i rrr diagram D. 

Jo p + + ^2 + ^(1 _ ^)Q2]2 [«2 + ^2]2 

A. 1.3 The M-channel diagrams 

We claimed before that the w-channel diagrams do not contribute to the imaginary part. Here 
we show that it is indeed so. 

The only difference between diagrams t-channel and w-channel diagrams is that in the 
li-channel case we switch the direction of particle p. Propagator (7) will be now 

rrfi + 13 s 13 + w? / s 

Following the same line as before, we recognize three regions A, B, and C. In region A (—1 < 
(3 < —rri^/s; —{3 < z <1) the last propagator turns into 



(7) ^ -K^ -m^ - 



s + <9 + + 



z z-\(3\ 



2 ^ 

m 



+ i^, (A. 16) 



which means we do not get any contribution to the imaginary part. One can check that no 
contribution to the imaginary part arises from the other regions on [z^ f3) plane (except for 
abnormal cases when ~ m*/s). The conclusion is that the li-channel diagram gives no 
leading - contributions to the imaginary part of the process. 
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A. 1.4 The same integral in the ^-representation technique 

Here, for the purpose of completeness, we show an alternative way to do the integrations over 
y and a. Sure enough, the underlying meaning of all manipulations will be absolutely the 
same as before. However, in a certain sense this way might appears simpler, since it does not 
require one to perform any graphic analysis or to think what propagator should be put on 
mass shell. Everything is done automatically here. It seems that this method is similar to the 
so-called a representation of the loop integrals. 

We start with expression (A.l) and use the following representation for each of the seven 
propagators 

1 i 

— = = i / dtexp[-te - itX] . (A. 17) 

X-te e + tX Jo ^ ^ ^ ^ 

The integral (A.l) transforms into 

/roo _^ 
dy dz da dj3 J dti- ■ -dtr exp[— tie — iti{ys + (5^)(1 — z) — iti{k'^ + m^)] 

X exp[— + it2yzs — it2{k^ + m^) — tse + iH{y + a){z + (3)s — it\[(k-\- + '^\[ 

X exp[— + it^yzs — it2{k^ + m^) — t^e + it^a(5s — it^{K^ + /i^)] 

X exp[— tge + iteaPs — ite{K^ + f/) — t-je + it'j{m'^ — /3s)(l — a) — ^^{k^ + m^)](A.18) 

Since the total integral is convergent and thanks to the factorization of y and a dependencies, 
we can do the y and a integration first and obtain 

j dy e^2/*[-*i(l-^)+*22:+t3(^+/3)+t4z] j gias[t3(z+;S)+t5^+t6,8+t7(;S-mVs)] ^ ^ 27ri ^ ^ 

X 5 [-ti(l -z)+ t2Z + t3{z + Uz] S [tsiz + /?) + t5(3 + te(3 + t7{(3 - mVs)jA.19) 

Since all parameters ti are positive, the two delta-functions can be simultaneously non-zero 
only for certain {z, (3) pairs. As should be expected, these pairs lie precisely within the allowed 
regions. 

Now we can take two of . . . integrations to kill the two delta-functions. Let them be 
^1 and is (of course, one can take other pairs as well). We get 

3 771^ 
ts = ^{t5+t6 + h)+t7- 



z + p' ' ' ^ ' ' '{z + (3)s' 

z B ni^ 

h = {t2 + U)- T^it5 + t6 + t7) + tr- -. (A.20) 

1 — z I — z (1 — z)s 

Since the integrations in ti remain factorized, we easily obtain 



JO 
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K^ + H^-ie- ^[(32(1 -z)+k^ + m^ + {k + /?)2 + m2 



r 

Jo 



he 



-I 



N? + jj?-ie- ^[Q\l -z)+P + w? + {k + Kf + w?] ' 

Jo 



—I 



If one now brings all pieces together, one will arrive at already familiar expression. 



(A.21) 



A. 2 Vector meson production: the fully off- forward case 

Since this case is done similarly to the Compton scattering, we do all calculations in less details 
and pay special attention only to distinctions. 

Strictly speaking, the (virtual) photoproduction of vector mesons is always an off-forward 
process, even though the final state might have no transverse momentum. The reason is that 
when the initial and final states have different masses, the exchange Pomeron must carry a 
non-zero longitudinal momentum. In the language of two-gluon exchange it means that the 
momenta of the gluons are not identical, in contract to what we had in previous section. 

The integral we deal with is 



/ dy dz da dp j-. — — 



r(M2 



[(g — ky — m? + ie][A;2 — w? + ie\[{k + A)^ — w? + ie][{k + — m? + ie] 
^ (A.22) 



[ki — ij? -\- ie\ [^2 — + ie] [{p — i^iY — 'rn? + ie\ 

With Sudakov decomposition and notation ni = n + A/2 ; = k — A/2, we have following 
seven propagators: 

^ -y - ^ j (1 - z)s - {P + m^) + ie , 

{y + a){z + 5)s - [(fc + A)^ + m^] + ie , 

1 1 ^ 

{y + a + -(T){z + (3 + -8)s-[{k + K + A/2)2 + m^] + ie , 

'J? + m^) + ie , 
1 , r..^ , 1an2 , ..21 , 



(1) 


(g — kY — + ie = 


(2> 


{k + A)^ - + ie = 


(3) 


{k + ACi)^ — + ie - 


(4> 


k"^ — + ie = yzs - 


(5) 


kI — /i^ + ie — {a + - 


(5) 


— jj? -\- ie — {a — 



2 /vr- 2 ' LV 2 



1 r..^ ^ , ,.21 
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(7) (p - Kif - + ie = (1 - a - ^(7) ^ « - [(« + ^ A)^ + m^] + ie . 

One can write now poles on the y and a plane, draw the regions of positivity of the correspond- 
ing ie's on the {z, (3) plane, and select the regions that leads to non-zero contributions (there 
will be five of them this time). One can again make sure that the leading 1/s contribution to 
the imaginary part comes from the main region. Another thing that helps avoid extra poles 
is the presence of /i^ (namely, /i^ > my — 4m^ for all vector mesons). 
Setting propagators (1) and (3) on mass shell, we obtain 

a 



y = 

s 



m 



[l-z)s 

The resulting quark propagators read 

1 



Q' 1 



k'^ + _^ (A; 



(4) 



(2) 



{l~z)s 

P + m'^ + z{l - z)Q'^] + ijl ; 



I- /? + A/2)2 + 
{z + (3 + S/2)s 



(A.23) 



1-z 

{z + 6)sa - {z + 5)sQ'^ 



(z + S) 



1-z 

Note that since the invariant mass of the produced qq state is equal 

^ ^ ,2 

:2 



[{k + Kf + m"] + ie . (A.24) 



^2 _ + ^ (A; + A)^ + m^ _ 



I- z ' z+S 

propagator (2) turns into {z + 6){my — M^). Together with the vertex factor r(M^), this 
propagator gives rise to the wave function. 
The gluon propagators are 



(5) = iP+ls){Q' 



P + w?' {k + K + K/2f 



+ 



m 



z + i3 + 5/2 



[(/5+^A)' + //']+i/; 



(6) 



(/3 - J) 



-m 



V 



P + m'^ {k + K + A/2f 



+ 



m 



z + (3+6/2 



P-iA)2 + /x2]+z/. 



In the case of gluon (5) the absence of pole was a trivial thing (due to negativity of both (3 
and S), in the case of the second gluon the pole can in principle arise. In our case however 
this is avoided due to large enough /x^ (namely, /x^ > my — 4m^ for all vector mesons). 



Finally, the propagator (7) acquires the form 



(7) = -(/3 + -<J) 



2 P+TV? {k + K + K/2Y + 771^ 

s 



l-z 



z + 13 + 5/2 



-[(« + ^A)2+//2]+ie. (A.25) 



This propagator, as usual, is used to extract the imaginary part of the integral, which leads 
to 



/3 



< 1. 



With all these manipulations done, one finally obtains 

r 



Im j dy dz da dp 



[all propagators] 



S 



2m \ ^ /■ dz 



-) ■ 



z{l-z) 



ipv{z,P 



[kl + m? + z{l - z)Q^] (/?2 + /x2)2 



(A.26) 
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Appendix B 

Helicity amplitude technique 



Here we give the derivation of expression for traces of the following type 

Sp |e {k^ + m) p (£3 + m) V* {k2 + m) p {ki + m)| 



(B.l) 



in full detail. Though one can calculate this trace covariantly, a particularly convenient way 
to do so is given by light cone helicity amplitude technique [25]. We emphasize that both 
ways are absolutely equivalent. In the helicity amplitude approach, we recall that all fermion 
lines in (B.l) can be taken on mass shell (see detailed derivation of LCWF normalization in 
Sect. 4. 5) and decomposed into a sum of hght cone hehcities 



{k + m) ^ ^ u\Ux for quark lines; 

— ^ vxv\ for antiquark lines. 

A=± 



{k + m)^-[{-k) 



Since the specific choice of this decomposition does not affect the final result, we are free to 
take the most convenient choice of spinors (see [25] for details), namely. 
Here are the spinors we use: 



u{p, A) 
v{p, A) 



1 



where 



1 



^^^^ ^ 71 



{\/2p'^ + (3m + ap 
(^\/2p'^ — (3m + apj 



x(T) 


A = 


+1 


xii) 


A = 


-1 


xii) 


A = 


+1 


x(T) 


A = 


-1 



(B.2) 



/1\ 



1 

V J 



^^^^ = 7=2 



( \ 

1 



V -1 / 



We stress that we use here our normal convention for + components of 4- vectors. For conve- 
nience, we also give the explicit expressions for all these spinors in the Dirac representation of 
7 matrices: 



I 
-I 



7 



—a 



a" 
k 



a 



a" 




7 



/ 

1 



(B.3) 
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The explicit expressions for all spinors W^ip), u {p), v (p), read: 



/ V2p+ + m\ 

Px + iPy 

•\/2p+ — m 
V Px + iPy J 



I -Px + iPy \ 

\/2p"*" + m 

Px - iPy 

y —\f2p'^ + m y 



m 



\ 



Px + iPy 
V Px + iPy J 



( -Px + iPy \ 

•\/2p+ — m 

Px - iPy 

y — V^p"*" — m J 



(B.4) 



Using these explicit formulas, one can do straightforward calculations and indeed prove 
that comply with orthonormality and completeness rules: 



u{p, X)u{p, A') = -v{p, \)v{p, A') = 2m6x\' ; 

u{p, X)u{p, X) = p + m; ^ v{p, X)v{p, A) = p — m. 



(B.5) 



For A = ±1 we defined 
Cross product is defined as 



a(A) = — Aoj 



2/ ■ 



— * 

[ah\ = aa;6j/ - h^a^ 



Moreover every matrix element should be multiplied by common factor ^p^q^. 
Some useful relations: 

l^vx{q) = -Xu^xiq) ; l^uxiq) = Xv^xiq) ; 
a(-A)6(A) = -db + iX[ab] ; ^a(-A)6(A) = -2ab; 



a(-A) (ba + iX[ba\) = b{-X)a'^ 



(B.6) 
(B.7) 



(B.8) 



a) 



P 




b) 




X . P_ 



d) 



Figure B.l: Four types of transitions, for which we give the amplitudes. Dot with label A 
indicates the spinor whose helicity is used as A. 
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Straight line elements 

For line in Fig. B. la, one has 

uxip)...ux{q) u-xip)---u\{q) 

7+ 2 

1 TTl 

y+H+ P+Q+ 



a ■ 7 

1 ^(^ + ^) 



p(A) g(A) 



7 7^7 ■ ^ 1^ {ap — iX[ap\) a(A) 



2m 



a-^-f^-f -^{aq + i\[a^) a(A) 

a -7 7+7-6 2{ab + i\[ab]) 

7+7^ 2A 



For line in Fig. B. lb, one has 



vx{p)...vx{q) vx{p)...v-x{q) 

7+ 2 

7~ vToT (^^ + P^- «A[pg]) — ^ (p(A) - g(A)) 



a • 7 



• - - p^q^ 

VP+ q+J 

. a ^{ap + iX[ap\) ^a(A) 



a-^-i^-f -^{aq-i\[a^) -^a(A) 

a-^-i^^-h 2{ah-i\[ah]) 

7+7^ -2A 



P+ 
2m 
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Vertex lines 

For vertex in Fig.B.lc, one has 



7+ 2 

t'-i-^-t P+q+ 



- - rx\ f 1 , l\ ap aq ( [ap\ [ag] 
a -7 ma( A) + ?A 



I p(A) _ g(A) _rn m 



P+ Q+ 



P+ ^ Q+ 



-i~-i^^-a ^«(A) —{dp-i\[dp\) 

U-^j+j-b 2(a6 + iA[a6]) 

7+7^ 2A 



For vertex in Fig.B.ld, one lias 

ux(p)...vx{q) ux(p)...v-x(q) 

7+ 2 

7~ (P(-A) + g(-A)) ^(-m'+pg + iA[pg]) 

a -7 -ma(-A) + 77- ^ iX { - — --- — - 

pj-X) g(-A) m^_m_ 

P+ Q+ 



1 



1' 



P+ Q+ 



7 7+7 -a -^a(-A) (ap-iA[ap]) 



a • 7 7+7- -|^^a(-A) — {aq + iX[aq\) 



P+ 
2 



• q^ 



a-7 7+7-6 2(a6 + iA[a6]) 

7+7^ 2A 

-A(^-^) Xm{- + - 
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Transverse photon poleirization 

As wc sec from the table, neither of gluon legs can flip quark helicity. Therefore, there are 4 
combinations of all possible helicity assignments. 

In the case both qq helicities are the same, one gets 



/ z{l — z)q+ 
2zq+p- ■ 2zq+p- ■ ^ ^ me* (A) 



z{l - z)q+ 



z{l-z)q+ 



— m)e(— A) 



l-z 



(B.9) 



In the case quark and antiquark helicities are opposite, one obtains 



2zq+p- ■ 2zq+p- ■ — ^ [{-e*k)z + {e*k){l - z) - iX{-[e*k]z - [e*k]{l - z))] 
^(1 - z)q+ 



z{l -z)q+ 



2^(1 - z)q+ 



[{ek){l -z)- {ek)z - iX{[e*k]{l - z) + [e*k]z)] 



l-z 



[{e*k){l - 2z) + iX[e*k]] [{ek){l - 2z) - iX[el 



(B.IO) 



Thus, performing summation over qq helicities, one obtains 



2s' 



1 - z 



m'(ee*) + (eA;)(e*A;)(l - 2zf + [ek][e*k] 



Finally, averaging over azymuthal angle, omitting factor 2s^ and including (—1) due to one 
antiquark propagator reveals 

l-z 



:/(«) = - 



m' + {z' + (1 - zf)k'' 



(B.ll) 



Note that the expression for Diagr.(d) can be instantly obtained from the above expressions 

— * — * 

by replacement A;— >— A;; z ^ 1 — z. In this case 3 antiquark propagators also give factor 
(— 1)^ = —1. So, the answer for this diagram is 



l-z 



^- rr^ + iz^ + il- zf)k'' 



(B.12) 



Finally, one has 

+ {z' + {I - zf){k,k + k) 
m^ + {z'' + {l-zf){k,k-K) 

for Diagrams (b) and (c) respectively. 



With all these results, Eq.(2.21) turns into 



IS 



X 



(27r) 



dz 



(Pk 



k^ + + z{l - z)Q^ J [K' + fi 



+ 



[A;2 + m2 + z{l - z)Q'^] [{k + kY + + z{l - z)Q'^] 
Note that we changed k — > —k in Diagr.(c), so that it became identical to Diagr.(b). 



(B.13) 
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Scaleir photon poleirization 

The virtual photon scalar polarization is described by 

1 1 



(B.14) 



One notes that in this case only qq states with opposite helicities do contribute. For Diagr. (a) 
an amplitude reads 



2zq+p- ■ 2zq+p- 

o 1 Z 



'^z{l-z)q+ 



Q 



q^{—m^ — k^) 
z{l - z)ql 



+ 2xp- 



q+{—m? — ]?) 
z{l - z)ql 



+ 2xp- 



— s 



m' + k'- z{l - z)Q' 



l2 



Therefore, 



l-z 



/(») = — 



Q 



^ V + P-z(i-^)g2i' 



Obviously, for Diagr. (d) one has 



l-z 



.fid) ^ t_ 



(B.15) 



(B.16) 



(B.17) 



The expressions for Diagrs. (b) and (c) can be obtained in the same way and give 
1 



1 



+ r - z{l - z)Q'\ [m' + {k + ny - z{l - z)Q' 
rn^ + k'^- z(l - z)Q'^] \m^ + (k - nf - z{l - z)Q' 



(B.18) 



Therefore Eq.(2.21) now reads 

32 i-) i 



IS 



2 



/ 



dz(fk 



[m^ + ^2 _ ^(1 _ ^)g2] ^ d?KV{K) 



k^ + + z{l - z)Q^ J (k^ + I^^) 



I 



[m^ + k^- z{l - z)Q^] ^ K + {k- k)2 - z{l - z)Q^] 



+ m2 + z{l - z)Q2] [(k - ny + m2 + z{l - z)Q^] 



(B.19) 



In fact, this expression in brackets can be greatly simplified. First, trivial algebra leads to 
[m2 + P- z{l - z)Q^\ [m^ + - - z{\ - z)Q''\ 



[F + m2 + z{\ - z)g2] p - ny + m2 + z{\ - z)Q'^\ 
2z(l-z)Q2] 2^(1 -^)g2 



-1 + ^ 



2z{\-z)Q 



[A;2 + m2 + z{l - z)g2] 
1 



+ 1- ^ 



p - k)2 + m2 + z(l - z)g2] 

1 



[A;2 + m2 + z{\ - z)Q2] p - ny + m2 + z{\ - z)Q^] 



Second, in we replace cq oc {q' + xp') hy q — {q' — xp'), we will get zero due to the gauge 
invariance of the totaf sum of all diagrams. Effectively, it means that residual [m^ + P — 
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z{l — z)Q'^] can be replaced by —2z{l — z)Q^. Of course, it does not mean that these two 
expressions are identical by themselves. Only after integration over quark momenta they give 
the same answer. Thus, the final answer for scalar photons reads 



1 



(27r) 



X 



1 



cPkV{k) 

2^2 



+ m2 + z{l - ^)Q2] [{k - + + z{l - z)Q^ 



B.l Photon vertex amplitudes 

Notation is given in Fig.5.3. We start with transverse photon. 

u'e-TV = u{—je)v . 

Equal qq helicities give 



^(1 - z)q+ 
^(1 - z)q+ 

Opposite qq helicities give 



(-m)e(-A) 



me(— A) 



^(1 - z)<l+ 
1 



{eki){l -z)- {eki)z - i\ ([^i](l -z) + [^i\z) 



{eki){l-2z) -iX[eki] 



In the case of scalar photon 

u'eov = u-^{q+j- + xp-j+)v 
the same helicities give exactly zero while opposite helicities result in 



^(1 - z)q+ 1 



^(1 - z)q+ Q 
11 

~~^z{l-z)Q 



— i—rn? — kV) + x2zil — z)p_q+ 



m^ + kf- z(l - z)Q' 



;b.2o) 



(B.21) 



(B.22) 



(B.23) 



B.2 Vector meson vertex amplitudes 

— * 

This case is more tricky due to the nonzero vector meson transverse momentum A. We start 
with the transverse vector meson polarization: 



^,1* 2(y*A) 
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The same qq helicities give again 



\/z(l — z)q+ 1 

^ -mV^*(A) = , mV*{X) (B.24) 



^(1 - z)q+ \/z{l-z) 
while opposite hehcities give 

^Jz{l - z)q+ 



^(1 - 



-[{V%)z + {V%)il -z)-iX {[V%]z - [V%]{1 - z))] 



p_2z(l - z)q+ 



{{V*k){l - 2z) + i\[Vk]} (B.25) 



Here we used definitions and properties (see also Fig. 5. 3): 

k'i^k + zK- k3 = -k+ (1- z)A (B.26) 

=^ {1 - z)k2 - zks = k ] {1- z)k2 + zk3 = {l-2z)k + 2z{l- z)A] 

{k2h) = -1? + (1 - 2z){kA) + z(l - z)A^ 

, ^9 1*9 ^9 +nn? H + P + 

+ = ^ + ^ ^ = ^ + . 

z (1 ~ -2) z[l — z) 

For the longitudinal vector mesons one has for equal quark-antiquark helicities 

^Jz{l - z)q+ 1 



-mA(A)+m^[A;2(A) + A;3(A)] 



= (B.27) 



z{l - z)q+ M 
and for opposite helicities 

^^/5iE^i^/_[(A4)^ + (a4)(1 -z)+ tX[Ak]] + ^i-m' + {k2h) + ^AM2]] 

+ — - — p-M^ - z)u } 

^ 2z(l - z)M (B.28) 



B.3 Final trace calculation 

Once we have the expressions for vertex amplitudes, the rest is done quickly. We first note 
that each gluon vertex attached to the lower or upper line gives factor 

2zq+-p- = sz- 2{l - z)q+ ■ p- = s{l - z) (B.29) 

correspondingly. 
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So. let's start with T ^ T amplitude and calculate it for Diagr.(c) at Fig.5.2. 
Equal qq helicities give 

s(l -z)-sz- ~^ mV*{\) ^ =me(-A) = -s^m^e{-X)V*{X) (B.30) 

^Jz{l-z) \lz{l-z) 

Summing over A gives 

2sW(ey*) . (B.31) 

The opposite helicities yield 

[{y*k){l - 2z) + i\{V*k^ [{eh){l - 2z) - iX[eh^ (B.32) 
Summing over helicities and making use of identity 

[ah\ [cd\ = {ac) (bd) - (ad) (be) (B.33) 

one obtains 

2s^ [{V*k){eki){l - 2zf + {eV*){kk^) - {ek){y*k^)\ . (B.34) 
Since we factored out when deriving (5.18), we finally get 

- \eV*)[m^ + kkx) + {y*k){ekx){\ - 2zf - {ek){y*kx)\ . (B.35) 



Ac) 



We included factor (—1) since in this diagram we have one antiquark propagator. 

An important observation here is that all other integrands, namely I^"'^ ■{l — z)/z, I^'^^ ■ 

— * 

z/{l — z) give absolutely the same result (with their own definitions of ki of course). The only 
thing one should not forget is that diagrams (a,d) enter with sign '— ' while diagrams (b,c) 
enter with sign '+': 

1 ~ ^_j{a) ^ j{b) ^ j{c) ^ ^_j{d) _ 



Z 1 — Z 

For L ^ L amplitude one has immediately 

/El = -^[m' + kl- z{l - z)Q'] ■ ^2z{l - z)M' (B.36) 
In fact, using simple relation 

w? + kl + z{l-z)Q'^ + kj + z{l - z)Q^ 

and noting that all unity terms will eventually cancel out in (5.24), one can rewrite (B.36) as 

I^lIl = -^QMz\l -zf. (B.38) 

For T ^ L amplitude one has 

41l = 2^(1 - z)M{eki) (1 - 2z) (B.39) 
and for L ^ T amplitude one has 

= 1[^2 ^ fc2 _ _ ^)g2](i _ 2z){V*k) . (B.40) 

The same transformation as in L — L amplitude, leads to 

I^^l^ = -2z{l - z)Q\l - 2z){V*k) . (B.41) 
Note that in the last three amplitudes only opposite qq helicities contributed. 
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